POJ 2115 C Looooops

转化成求解a+cx=b(mod 2^k)的最小整数解，还是求解一次同余式

 

注意当k=32时，int会溢出

 

代码：

#include<iostream> using namespace std; __int64 a,b,c,k,i,d,x,y,ans; __int64 extended_eculid(__int64 a,__int64 b,__int64 &x,__int64 &y) { if(b==0) { x=1; y=0; return a; } __int64 t,res; res=extended_eculid(b,a%b,x,y); t=x; x=y; y=t-a/b*y; return res; } int main() { while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&i)) { if(a==0&&b==0&&c==0&&i==0) break; k=((__int64)1)<<i; b-=a; d=extended_eculid(c,k,x,y); if(b%d!=0) { printf("FOREVER/n"); continue; } ans=x*b/d; k/=d; ans=(ans%k+k)%k; cout<<ans<<endl; } return 0; }