首页 > 代码库 > POJ 2115 C Looooops(扩展欧几里得应用)
POJ 2115 C Looooops(扩展欧几里得应用)
题目地址:POJ 2115
水题。。公式很好推。最直接的公式就是a+n*c==b+m*2^k.然后可以变形为模线性方程的样子,就是
n*c+m*2^k==b-a.即求n*c==(b-a)mod(2^k)的最小解。(真搞不懂为什么训练的时候好多人把青蛙的约会都给做出来了,这题却一直做不出来。。。。。这两道不都是推公式然后变形吗。。。。。)
代码如下:
#include <iostream> #include <cstdio> #include <string> #include <cstring> #include <stdlib.h> #include <math.h> #include <ctype.h> #include <queue> #include <map> #include <set> #include <algorithm> using namespace std; #define LL __int64 LL X, Y; LL exgcd(LL a,LL b) { if(b==0) { X=1; Y=0; return a; } LL r=exgcd(b,a%b); LL t=X; X=Y; Y=t-a/b*Y; return r; } int main() { LL a, b, c, k, d, L, z; while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k)!=EOF&&k) { z=pow(2,k); d=exgcd(c,z); L=b-a; if(L%d) { printf("FOREVER\n"); continue ; } else { LL ans=X*L/d; LL s=z/d; if(ans<0) { ans=ans%s+s; } else { ans=ans%s; } printf("%I64d\n",ans); } } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。