POJ 2891 Strange Way to Express Integers 中国剩余定理解法
一种不断迭代,求新的求余方程的方法运用中国剩余定理。
总的来说,假设对方程操作。和这个定理的数学思想运用的不多的话。是非常困难的。
參照了这个博客的程序写的: http://scturtle.is-programmer.com/posts/19363.html
这个博客举例说的挺好的:http://blog.csdn.net/mishifangxiangdefeng/article/details/7109217
hdu 3579 Hello Kiki 中国剩余定理(不互质的情况)
对互质的情况,处理起来比較方便,能够直接套模板
本题给出不互质的模线性方程组,求出满足方程的最小正整数解
方案:对于不互质的模线性方程组,能够进行方程组合并。求出合并后的方程的解。这样就能够非常快地推出方程的终于解。
两个方程合并的一种方法:
x = c1 (mod b1)
x = c2(mod b2)
此时b1,b2不必互质的。
显然能够得到x = k1 * b1 + c1 x = k2* b2 + c2。
两个方程合并一下就能够得到:k1 * b1 = c2 - c1 (mod b2),
这样能够设g=gcd(b1,b2),于是就有b1/g*k1-b2/g*k2=(c2-c1)/g,
显然推断(c2-c1)/g是否为整数就能推断是否存在解。
这样在经过类似的变换就能得到k1 = K (mod (b2/g)),
最后得到x = K*b1 + c1 (mod (b1 * b2/g))。
对于题目所给正整数的要求,仅仅有一种反例,就是结果输出为0的情况,
这个能够特殊考虑。仅仅须要考虑全部数的最小公倍数就可以。
各个式子各个变量的含义都须要理解才干写好这个程序;最后0MS过。这个程序竟然上榜了。
__int64 s, t, g; void extGCD(__int64 a, __int64 b)
{
if (b == 0)
{
s = 1, t = 0, g = a;
}
else
{
extGCD(b, a % b);
__int64 tmp = s;
s = t;
t = tmp - a / b * t;
}
} int main()
{
__int64 m1, m2, r1, r2, m10, m20, c;
int n; while (scanf("%d", &n) != EOF)
{
bool flag = false;
scanf("%lld %lld", &m1, &r1); for (int i = 1; i < n; i++)
{
scanf("%lld %lld", &m2, &r2);
if (flag) continue;
extGCD(m1, m2);//由于定理条件是除数互质,所以除以公约数使得其互质
c = r2 - r1;//k1*m1 == (r2 - r1) (mod m2)
if (c % g)
{
flag = true;
continue;
}
m20 = m2 / g;//这个为新的mod除数。和以下新的m1互质
c /= g;
__int64 r0 = (c * s % m20 + m20) % m20;
r1 = r0 * m1 + r1;
m1 = m1 * m20;//得到新式子的系数: m1 * x + r1 == r2 即:x = r1, r2...(mod m1, m2)
}
if (flag) puts("-1");
else printf("%lld\n", r1);
}
return 0;
}
POJ 2891 Strange Way to Express Integers 中国剩余定理解法的更多相关文章
- POJ 2891 Strange Way to Express Integers 中国剩余定理 数论 exgcd
http://poj.org/problem?id=2891 题意就是孙子算经里那个定理的基础描述不过换了数字和约束条件的个数…… https://blog.csdn.net/HownoneHe/ar ...
- POJ 2891 Strange Way to Express Integers(中国剩余定理)
题目链接 虽然我不懂... #include <cstdio> #include <cstring> #include <map> #include <cma ...
- poj 2981 Strange Way to Express Integers (中国剩余定理不互质)
http://poj.org/problem?id=2891 Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 13 ...
- poj 2891 Strange Way to Express Integers (非互质的中国剩余定理)
Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 9472 ...
- poj——2891 Strange Way to Express Integers
Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 16839 ...
- [POJ 2891] Strange Way to Express Integers
Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 10907 ...
- poj 2891 Strange Way to Express Integers(中国剩余定理)
http://poj.org/problem?id=2891 题意:求解一个数x使得 x%8 = 7,x%11 = 9; 若x存在,输出最小整数解.否则输出-1: ps: 思路:这不是简单的中国剩余定 ...
- [poj 2891] Strange Way to Express Integers 解题报告(excrt扩展中国剩余定理)
题目链接:http://poj.org/problem?id=2891 题目大意: 求解同余方程组,不保证模数互质 题解: 扩展中国剩余定理板子题 #include<algorithm> ...
- poj 2891 Strange Way to Express Integers【扩展中国剩余定理】
扩展中国剩余定理板子 #include<iostream> #include<cstdio> using namespace std; const int N=100005; ...
随机推荐
- js 计算时间差
function GetDateDiff(startDate,endDate) { var startTime = new (<any>Date)(Date.parse(startDate ...
- 在Swift中,如何像Objective-C定义可选接口?
Objective-C中的protocol里存在@optional关键字,被这个关键字修饰的方法并非必须要被实现.我们可以通过接口定义一系列方法,然后由实现接口的类选择性地实现其中几个方法.在Coco ...
- cocos2dx实现单机版三国杀(二)
接上续,东西还没有做完 所以代码免不了改动 之前的头文件现在又改了不少,因为架构也改变了现在主要类就是GameScene.GameUI.PlayInfo.Poker这四个类 前面想的GameLoop ...
- jQuery中国各个省份地图分部代码
jQuery中国各个省份地图分部代码 在线演示本地下载
- 【PostgreSQL-9.6.3】触发器概述(普通触发器)
一个触发器声明了当执行一种特定类型的操作时数据库应该自动执行一个特殊的函数.触发器可以被附加到表.视图和外部表.触发器经常用于做完整性约束,或者某种业务规则的约束. 1. 触发器的创建语法如下: CR ...
- Prime算法生成最小生成树
虽说是生成树,但我只将生成的边输出了.至于怎么用这些边来创建树...我不知道_(:з」∠)_ //Prime方法生成最小生成树 void GraphAdjacencyListWeight::Gener ...
- HDU_1026_Ignatius and the Princess I_BFS(保存路径)
Ignatius and the Princess I Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (J ...
- Cesium学习笔记(四)Camera ----http://blog.csdn.net/hobhunter/article/details/74909641
Cesium 相机控制场景中的视野.操作相机的方法有很多,如旋转,缩放,平移和飞到目的地.Cesium具有默认的鼠标和触摸事件处理程序与相机进行交互,还有一个API以编程方式操纵相机. 我们可以使用该 ...
- IntentService和HandlerThread的使用以及源码阅读
使用MyIntentService.java public class MyIntentService extends IntentService { /** * 是否正在运行 */ private ...
- Notepad++运行JAVA代码
第一种方法: 工具栏->运行 点击后选择运行 1.在运行窗口中输入: cmd /k javac "$(FULL_CURRENT_PATH)" & echo 编译成功 ...