以前的做法

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
LL ai[4],ri[4],M;
void Exgcd(LL a,LL b,LL& d,LL& x,LL& y)
{if(b == 0) { d = a, x = 1, y = 0;}else {Exgcd(b,a%b,d,y,x); y -= x*(a / b);}
}
LL gcd(LL a,LL b)
{return b == 0 ? a : gcd(b,a%b);
}
LL solve(int n)
{LL a1,r1;LL a, b, c;LL x, y, d;a1 = ai[1], r1 = ri[1];for(int i = 2; i <= n; ++i){a = a1, b = ai[i], c = ri[i] - r1;Exgcd(a,b,d,x,y);if(c % d) return -1;LL t = b / d;x = (x * (c / d)% t + t) % t;r1 = a1 * x + r1;a1 = a1 *(ai[i] / d); }return r1;
}
int main()
{LL p,e,k,d;int Kase = 1;while(cin >> p >> e >>k >>d){ai[1] = 23, ai[2] = 28, ai[3] = 33;ri[1] = p , ri[2] = e, ri[3] = k;if(p == -1 && e == -1 && k == -1 && d == -1) break;LL ans = solve(3);M = 1;for(int i = 1; i <= 3; ++i) M = M * ai[i] / gcd(M,ai[i]);if(ans >= M) ans %= M; // 求最小的值. while(ans <= d) ans += M;printf("Case %d: the next triple peak occurs in %lld days.\n",Kase++,ans-d);}
}

中国剩余定理撸一发

不过中国剩余定理 ai[i] 需要互质

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
LL ai[4],ri[4],M;
void Exgcd(LL a,LL b,LL& d,LL& x,LL& y)
{if(b == 0) { d = a, x = 1, y = 0;}else {Exgcd(b,a%b,d,y,x); y -= x*(a / b);}
}
LL gcd(LL a,LL b)
{return b == 0 ? a : gcd(b,a%b);
}
LL China(int n)
{M = 1;LL Mi,x,y,d,ans = 0;for(int i = 1; i <= n; ++i) M *= ai[i];for(int i = 1; i <= n; ++i){Mi = M / ai[i];Exgcd(Mi,ai[i],d,x,y);ans = (ans + Mi * x * ri[i]) % M;}if(ans < 0) ans += M;return ans ;
}
int main()
{LL p,e,k,d;int Kase = 1;while(cin >> p >> e >>k >>d){ai[1] = 23, ai[2] = 28, ai[3] = 33;ri[1] = p , ri[2] = e, ri[3] = k;if(p == -1 && e == -1 && k == -1 && d == -1) break;LL ans = China(3);/*M = 1;for(int i = 1; i <= 3; ++i) M = M * ai[i] / gcd(M,ai[i]);*///if(ans >= M) ans %= M; // 求最小的值. while(ans <= d) ans += M;printf("Case %d: the next triple peak occurs in %lld days.\n",Kase++,ans-d);}
}

转载于:https://www.cnblogs.com/aoxuets/p/5506906.html

POJ 1006 同余方程组相关推荐

  1. 初等数论--同余方程--同余方程组:中国剩余定理

    初等数论--同余方程--同余方程组:中国剩余定理 博主是初学初等数论(整除+同余+原根),本意是想整理一些较难理解的定理.算法,加深记忆也方便日后查找:如果有错,欢迎指正. 我整理成一个系列:初等数论 ...

  2. [XSY] 智慧树(线性同余方程组,线段树/树状数组)

    智慧树 解决此题有两个要点: 如何判断一个线性同余方程组有没有解 如何统计合法子序列数目 先看第2点: 若一个序列是合法的,则这个序列的所有子序列都是合法的 考虑对∀1≤i≤n\forall 1\le ...

  3. 【中国剩余定理】POJ 1006 HDU 1370 Biorhythms

    题目链接: http://poj.org/problem?id=1006 http://acm.hdu.edu.cn/showproblem.php?pid=1370 题目大意: (X+d)%23=a ...

  4. poj 1006 生理周期

    题目链接:http://poj.org/problem?id=1006 题意:中文题. 中国剩余定理: 1 #include <cstdio> 2 #include <cmath&g ...

  5. 生理周期,POJ(1006)

    题目链接:http://poj.org/problem?id=1006 解题报告: 1.枚举天数的时候可以根据前面的结果直接跳过一些错误的答案. ///三个周期是23,28,33, #include ...

  6. Poj 1006 / OpenJudge 2977 1006 Biorhythms/生理周期

    1.链接地址: http://poj.org/problem?id=1006 http://bailian.openjudge.cn/practice/2977 2.题目: Biorhythms Ti ...

  7. POJ 1006 Biorhythms

    题目大意: 人生来就有三个生理周期,分别为体力.感情和智力周期,它们的周期长度为23天.28天和33天每一个周期.中有一天是高峰.在高峰这天,人会在相应的方面表现出色.例如,智力周期的高峰,人会思维敏 ...

  8. poj 1006(中国剩余定理)

    中国剩余定理: <孙子算经>中有"物不知数"问题:"今有物不知其数,三三数之余二 ,五五数之余三 ,七七数之余二,问物几何?"答为"23& ...

  9. poj 1006 java_POJ 1006 Java:中国剩余定理

    题目描述: 人生来就有三个生理周期,分别为体力.感情和智力周期,它们的周期长度为23天.28天和33天.每一个周期中有一天是高峰.在高峰这天,人会在相应的方面表现出色.例如,智力周期的高峰,人会思维敏 ...

最新文章

  1. 使用Teams接收Azure服务运行状况通知
  2. 哈希扩展长度攻击_哈希长度扩展攻击
  3. LeetCode 130. 被围绕的区域(图的BFS/DFS)
  4. php sqrt函数,sqrt函数怎么使用
  5. HAProxy负载均衡代理
  6. php中crypt怎么还原,PHP crypt() 函数 - PHP 教程 - 菜鸟学堂-脚本之家
  7. turtlebot3_teleop_key 键盘控制程序,使之0速度时不会一直发送 topic,通过topic代替键盘控制小车运动。
  8. 学习笔记:Oracle的trace文件可见性
  9. Android项目:通过ant重新打包proguard混淆器jar文件
  10. java 反射调用方法_Java 反射详解,重要方法解析
  11. php订单号生成规则,项目笔记之订单号生成规则以及方法,第一篇!
  12. SQL SERVER2008查询分析器的基本使用
  13. 机器学习入门之莺尾花训练
  14. 著者四角号码查询_著者姓名汉语拼音与四角号码数字混编书次号的应用
  15. 在Win10 WSL中安装百度Apollo
  16. python bar函数循环_python bar函数怎么使用
  17. 这是一个谷歌抄腾讯的时代
  18. 通达OA v11.7后台SQL注入
  19. arcgis表转excel一直失败_arcGIS中属性表怎么导出到EXCEL里|excel表格转shp数据丢失...
  20. 贪心法——活动安排问题

热门文章

  1. cPanel附加域名出现Error from park wrapper: 使用带以下 IP 的命名服务器:
  2. vc实现html显示,在VC对话框中用ChtmlView控件显示HTML(2)
  3. 使用python制作ArcGIS插件(5)其他技巧
  4. Go语言中的复合类型及面向对象思想
  5. DXUT框架剖析(4)
  6. 【java学习之路】(java框架)008.JdbcTemplate
  7. oracle非延迟约束,Oracle可延迟约束Deferable的使用
  8. java 实例域_Java实例域初始化
  9. 2018 蓝桥杯省赛 B 组模拟赛(五)题 B题
  10. LeetCode81. 搜索旋转排序数组 II(二分查找)