Poj 2115 C Looooops (模线性方程)
来源:互联网 发布:淘宝优惠券怎么发放 编辑:程序博客网 时间:2024/05/22 05:01
题意:对于循环语句for(i=A ; i!=B ;i +=C),问在无符号的k位存储系统中循环几次才会结束。若在有限次内结束,则输出循环次数,否则输出死循环。
思路:例如k=16,也就是unsigned int,那么能保存2^16个数据,即最大数为65535,当循环使得i超过65535时,则i会返回0重新开始计数。
于是可得 循环次数 x=[(B-A+2^k)%2^k] /C
即 Cx=(B-A)(mod 2^k) ,求最小正整数解。
以下参考:http://972169909-qq-com.iteye.com/blog/1104538
解方程:ax == b (mod n);【ax % n == b % n】
设线性模方程的一个解为x0
条件①:有d = gcd(a, n)
条件②:有d = ax1 + ny, 由扩展欧几里得(Egcd)得到x1的值
条件③:b % d == 0 (有解的条件)
对条件③进行解释:
原方程化为:ax + kn = b (设k为某一整数)
那么如果a与n的最大公约数为d,那么ax + kn 必然可以提取一个d的因子,也就是说b必然有d这个因子,所以如果b%d!=0,说明b没有d这因子,与前面的结论相互矛盾,所以无解
则x0 = x1*(b/d);
证明:
因为:容易求得d = gcd (a, n), 则存在一个x1、y使得d = ax1 + ny①
方程①2边同时模n得:d % n == ax1 % n②
又因为:b % d == 0, 即b是d的倍数;
所以(b/d)必为整数;
所以由②得: b % n == b*d/d%n == d*(b/d) % n == ax1*(b/d) % n == ax % n
所以很容易可以看出x = x1*(b/d)是方程的一个整数解,得证
#include <cstdio>#include <cstring>#define i64 __int64i64 Extended_Euclid (i64 a,i64 b,i64 &x,i64 &y){//扩展欧几里得算法,求ax+by=gcd(a,b)的一组解(x,y),d=gcd(a,b) i64 d; if (b==0) { x=1;y=0; return a; } d=Extended_Euclid(b,a%b,y,x); y-=a/b*x; return d;}bool MLE (i64 a,i64 b,i64 n,i64 &x0) //解模线性方程:ax == b (mod n);【ax % n == b % n】{//无解返回false,有解返回true,解存放在x0i64 d,x,y;d=Extended_Euclid (a,n,x,y);if (b % d)return false;x0=x * (b/d);i64 t=n/d;if (t<0) t=-t;//以防万一,有的题目t有可能是负数x0=(x0%t + t)%t;return true;}int main (){i64 a,b,c,k,x0;while (scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k),a||b||c||k){if (MLE(c,b-a,(i64)1<<k,x0))printf("%I64d\n",x0);elseprintf("FOREVER\n");}return 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 (单变元模线性方程)
- C Looooops--模线性方程
- [ACM] POJ 2115 C Looooops (扩展欧几里得求解模线性方程)
- pku 2115 C Looooops(解模线性方程)
- POJ2115 C Looooops 一元模线性方程
- C Looooops poj2115 (扩展欧几里得+模线性方程)
- poj 2115 C Looooops
- So cute are you python 1
- 6410的存储空间分配
- Linux 系统调用
- HTTP协议详解(真的很经典)
- 如何利用C/C++逐行读取txt文件中的字符串(可以顺便实现文本文件的复制)
- Poj 2115 C Looooops (模线性方程)
- Leetcode: Sqrt(x)
- 农夫过河问题(一个有趣的问题-位运算)
- 大话数据结构五:线性表的链式存储结构(双向链表)
- 最长公共子字符串
- MINOS(IIS) 软件 -用于线性问题的不可行解分析
- 关于b=++b和b = b += 1的问题探讨
- HTTP协议详解
- 【VC++6.0】getline需要输入2次回车才会结束的BUG修复方法