POJ 2115（解线性同余方程）

与POJ 1061青蛙很像。

注意：for（int i=2；i！=1；i++） 循环到 i 值溢出时，会重新从0开始加。

#include<cstdio>#include<cmath>#include<algorithm>#define ll long longusing namespace std;void gcd(ll a,ll b,ll &d,ll &x,ll &y){    if(!b){d=a;x=1;y=0;}    else{ gcd(b,a%b,d,y,x);y-=x*(a/b);}}int main(){    ll a,b,c,k;    while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&k)&&(a||b||c||k)){        ll M = 1ll<<k;        ll d,x,y;        gcd(c,M,d,x,y);        ll C = b-a;        if(C%d!=0){        //如果推导的公式中的B-A不是d的倍数，则一定无解            printf("FOREVER\n");        }else{            x = (x*(C/d))%M;          //gcd中的x的最小整数解            ll ans = (x%(M/d)+M/d)%(M/d);   //通式（结果）的最小正整数解。            printf("%lld\n",ans);        }    }    return 0;}