C Looooops DescriptionA Compiler Mystery: We are given a C-language style for loop of type for (variable = A; variable != B; variable += C) statement;I.e., a loop which starts by setting variable to value A and while variable is not equal to B, repea…

poj2115 C Looooops 题意: 对于C的for(i=A ; i!=B ;i +=C)循环语句,问在k位存储系统中循环几次才会结束. 若在有限次内结束,则输出循环次数. 否则输出死循环. (k位==mod $2^{k}$) 列出方程:$A+Cx\equiv B(mode\quad 2^{k})$ 转换一下:$Cx+ky=B-A$ 用exgcd解出 $Cx+ky=gcd(C,k)$ 然后把求出的$x*(B-A)/gcd(C,k)$ 再$\% (k/gcd(C,k))$求个最小正整数解…

题目链接:https://vjudge.net/problem/POJ-2115 题意:模拟for循环for(int i=A;i!=B;i+=C),且数据范围为k位无符号数以内,即0~1<<k-1,如果能循环为有限次则输出循环次数,否则输出FOREVER. 思路:典型的扩展欧基里德题.题意即求Cx=B-A (mod 1<<k),可化为Cx+(1<<k)y=B-A (mod 1<<k).不访令a=C,b=1<<k,c=B-A,即求ax+by=c (…

欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - POJ2115 题意 对于C的for(i=A ; i!=B ;i +=C)循环语句,问在k位存储系统中循环几次才会结束.若在有限次内结束,则输出循环次数.否则输出死循环. 题解 原题题意再次缩略: A + xC Ξ B (mod 2k) 求x的最小正整数值. 我们把式子稍微变一下形: Cx + (2k)y = B-A 然后就变成了一个基础的二元一次方程求解,扩展欧几里德套套就可以了. 至于扩展欧几里德(ex…

C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 24355   Accepted: 6788 Description A Compiler Mystery: We are given a C-language style for loop of type for (variable = A; variable != B; variable += C) statement; I.e., a loop w…

无符号k位数溢出就相当于mod 2k,然后设循环x次A等于B,就可以列出方程: $$ Cx+A \equiv B \pmod {2^k} $$ $$ Cx \equiv B-A \pmod {2^k} $$ 最后就用扩展欧几里得算法求出这个线性同余方程的最小非负整数解. #include<cstdio> #include<cstring> #define mod(x,y) (((x)%(y)+(y))%(y)) #define ll long long ll exgcd(ll a,…

题意:很明显,我就不说了 分析:令n=2^k,因为A,B,C<n,所以取模以后不会变化,所以就是求(A+x*C)%n=B 转化一下就是求 C*x=B-A(%n),最小的x 令a=C,b=B-A 原式等于ax=b(mod n) 这就是标准的解模线性方程 该方程有解的充要条件是d=gcd(n,a) && d|b(可以根据这一条判断是否FOREVER) 然后参考算法导论应用扩展欧几里得求解x a*x0+n*y0=d x=x0*(b/d)(mod n) 然后应用多解条件求最小正整数解,即解的…

题目大意 求同余方程Cx≡B-A(2^k)的最小正整数解 题解 可以转化为Cx-(2^k)y=B-A,然后用扩展欧几里得解出即可... 代码: #include <iostream> using namespace std; typedef long long LL; void extended_gcd(LL a,LL b,LL &d,LL &x,LL &y) { if(!b) { d=a,x=,y=; } else { extended_gcd(b,a%b,d,y,x…

题目链接. 分析: 数论了解的还不算太多,解的时候,碰到了不小的麻烦. 设答案为x,n = (1<<k), 则 (A+C*x) % n == B 即 (A+C*x) ≡ B (mod n) 化简得 C*x ≡ (B-A) (mod n) 设 a = C, b = (B-A) 则原式变为 ax=b.解 x. 到这里,以为求出来 a 的逆, 然后 x = b*a-1. a 的 逆好求,用<训练指南>上的模板(P122) inv函数. 例如,求 2 模 66536 下的逆, inv(2,…

C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 23700   Accepted: 6550 Description A Compiler Mystery: We are given a C-language style for loop of type for (variable = A; variable != B; variable += C) statement; I.e., a loop w…