简单扩欧悬关求调
  • 板块学术版
  • 楼主_lqs_
  • 当前回复15
  • 已保存回复15
  • 发布时间2023/8/21 19:31
  • 上次更新2023/11/3 02:09:58
查看原帖
简单扩欧悬关求调
664744
_lqs_楼主2023/8/21 19:31

ABC315G AC 55,WA 46 调了一年。应该是边界问题但是没看出来(

#include<bits/stdc++.h>
using namespace std;
#define N 1000005
__int128 n,m,i,j,ans,a,b,c,x,l,r;
__int128 gcd(__int128 a,__int128 b){
	if(b==0) return a;
	return gcd(b,a%b);
}
void exgcd(__int128 a,__int128 b,__int128 c){
	if(b==0){
		l=c/a,r=0;
		return;
	}
	exgcd(b,a%b,c);
	__int128 tmp=l;
	l=r,r=tmp-a/b*r;
}
inline __int128 read()
{
    __int128 x=0,f=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
    {
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}

inline void write(__int128 x)
{
    if (x<0) putchar('-'),x=-x;
    if (x>9) write(x/10);
    putchar(x%10+'0');
}
signed main(){
	n=read(),a=read(),b=read(),c=read(),x=read();
	for(i=1;i<=n;i++){
		__int128 k=x-a*i,d=gcd(b,c);
		if(k%d!=0 || k<=0) continue;
		exgcd(b,c,k);
		__int128 L=c/d,R=b/d;
		l%=L;
		if(l==0) l+=L;
		r=(k-b*l)/c;
		if(r<=0) continue;
		if(r>n){
			if((r-n)%R==0) r=n;
			else r=r-((r-n)/R+1)*R;
			l=(k-c*r)/b;
			if(l>n || r<=0) continue;
		}
		__int128 p1=0,p2=0;
		if(r%R==0) p1=r/R;
		else p1=r/R+1;
		if((n-l)%L==0) p2=(n-l)/L+1;
		else p2=(n-l)/L+1;
		ans+=min(p1,p2);
//		printf("%d\n",ans);
	}
	write(ans);
	return 0;
}
2023/8/21 19:31
加载中...