poj 2115 C Looooops 扩展欧几里得
来源:互联网 发布:黑色星期五禁曲 知乎 编辑:程序博客网 时间:2024/04/30 15:46
题意:
在while(x=a;x!=b;x+=c) statement;中,问statement会被执行多少次,计算在2^k下进行。
思路:
等价于计算同余式a+c*x=b(mod2^k)用扩展欧几里得算法。设g=gcd(a,b)在计算a*x+b*y=g过程中,x的结果可以用b/g调整,y的结果可以用a/g调整,因为a*(b/g)==b*(a/g)。
代码:
//poj 2115//sep9#include<iostream>using namespace std;void gcd(__int64 a,__int64 b,__int64 &d,__int64 &x,__int64 &y){if(b==0){d=a,x=1,y=0;}else{//a/b==(a-a%b)/bgcd(b,a%b,d,y,x);//by+(a%b)x=d ==> ax+b(y-x(a-a%b)/b))=by+(a%b)x=dy-=x*(a/b);}}int main(){__int64 a,b,c,A,B,C,k,d,x,y;while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k)==4){if(a==0&&b==0&&c==0&&k==0)break;k=1LL<<k;A=c;B=-k;C=b-a;gcd(A,B,d,x,y);if(C%d!=0)printf("FOREVER\n");else{x*=C/d;y*=C/d;if(d<0)d=-d;x=x%(k/d);if(x<0)x+=(k/d);printf("%I64d\n",x);}}return 0;}
0 0
- POJ 2115 C Looooops(扩展欧几里得)
- poj 2115 C Looooops(扩展欧几里得)
- poj-2115 C Looooops(扩展欧几里得)
- poj 2115 C Looooops (扩展欧几里得)
- POJ 2115 C Looooops 扩展欧几里得算法
- poj 2115 C Looooops 扩展欧几里得
- POJ - 2115 - C Looooops (扩展欧几里得)
- POJ 2115 C Looooops(扩展欧几里得)
- POJ 2115 C Looooops(扩展欧几里得)
- poj 2115 C Looooops 【扩展欧几里得】
- POJ 2115 C Looooops (扩展欧几里得)
- poj 2115 C Looooops (扩展欧几里得算法)
- POJ-2115-C Looooops(扩展欧几里得)
- C Looooops poj 2115 扩展欧几里得
- [poj 2115] C Looooops 扩展欧几里得
- poj 2115 C Looooops (扩展欧几里得)
- 【POJ】2115 - C Looooops(扩展欧几里得)
- poj 2115 C Looooops 扩展欧几里得算法
- Unity之Handles旋转控制器-十三
- CSS样式
- 【记录】LINK : fatal error LNK1123: failure during conversion to COFF: file invalid or corrupt
- 在CentOS6.2下安装DNS服务软件Bind并快速配置简单实例
- webpy源码阅读(2)——CherryPyWSGIServer
- poj 2115 C Looooops 扩展欧几里得
- Unity之Handles缩放控制柄-十四
- 【记录】VIM的基本操作
- java inputStream之read、available
- 矩阵链乘法(递归法)
- C++ 基础回顾
- USACO 1.2 Problem 1
- Unity之Handles绘制球体控制柄-十五
- 如何将PSD模型转化成XHTML和CSS文件