CodeVS1990 中国余数定理2
2016-04-09 16:01
197 查看
3990 中国余数定理
2
题目描述 Description
Skytree神犇最近在研究中国博大精深的数学。
这时,Sci蒟蒻前来拜访,于是Skytree给Sci蒟蒻出了一道数学题:
给定n个质数,以及k模这些质数的余数。问:在闭区间[a,b]中,有多少个k?最小的k是多少?
Sci蒟蒻数学能力差了Skytree三条街,所以他只好寻求计算机的帮助。他发邮件给同为oier的你,你能帮他解决这个问题吗?
输入描述 Input Description
输入第一行为三个正整数n、a、b。
第2到n+1行,每行有两个整数,分别代表第n个质数和k模第n个质数的余数。
输出描述 Output Description
输出为两个整数,代表闭区间[a,b]中k的个数和闭区间[a,b]中最小的k。如果k不存在,则输出两个0。
样例输入 Sample Input
样例1:
3 2 28
3 2
5 3
7 2
样例2:
3 24 31
3 2
5 3
7 2
样例输出 Sample Output
样例1:
1
23
样例2:
0
0
数据范围及提示 Data Size & Hint
1<=a<=b<=10^14
n<=10
输入保证所有n个质数的乘积<=10^14
每个质数<=1.5*10^9
请无视通过率(被人黑了。。。)
数据保证不会溢出64bit整数
模线性方程组解的构造过程如下,
令ji=m1*m2*...*mn,
Mi=ji/mi
ti=Mi在模mi意义下的逆元
方程的最小解如下
ans=ai*ti*Mi%ji
方程的通解为 ai * ti * mi % ji + k * ji
其中k是非负整数
2
题目描述 Description
Skytree神犇最近在研究中国博大精深的数学。
这时,Sci蒟蒻前来拜访,于是Skytree给Sci蒟蒻出了一道数学题:
给定n个质数,以及k模这些质数的余数。问:在闭区间[a,b]中,有多少个k?最小的k是多少?
Sci蒟蒻数学能力差了Skytree三条街,所以他只好寻求计算机的帮助。他发邮件给同为oier的你,你能帮他解决这个问题吗?
输入描述 Input Description
输入第一行为三个正整数n、a、b。
第2到n+1行,每行有两个整数,分别代表第n个质数和k模第n个质数的余数。
输出描述 Output Description
输出为两个整数,代表闭区间[a,b]中k的个数和闭区间[a,b]中最小的k。如果k不存在,则输出两个0。
样例输入 Sample Input
样例1:
3 2 28
3 2
5 3
7 2
样例2:
3 24 31
3 2
5 3
7 2
样例输出 Sample Output
样例1:
1
23
样例2:
0
0
数据范围及提示 Data Size & Hint
1<=a<=b<=10^14
n<=10
输入保证所有n个质数的乘积<=10^14
每个质数<=1.5*10^9
请无视通过率(被人黑了。。。)
数据保证不会溢出64bit整数
分析
废话不多说,这是中国剩余定理模板题,直接看过程模线性方程组解的构造过程如下,
令ji=m1*m2*...*mn,
Mi=ji/mi
ti=Mi在模mi意义下的逆元
方程的最小解如下
ans=ai*ti*Mi%ji
方程的通解为 ai * ti * mi % ji + k * ji
其中k是非负整数
代码
//CodeVS3990 中国余数定理2 中国剩余定理 #include <cstdio> #define ll long long using namespace std; ll N, A, B, a[15], m[15], cnt, minn, t[15], M[15], ji, ans, tmp; void exgcd(ll a, ll b, ll &x, ll &y) { if(!b){x=1,y=0;return;} exgcd(b,a%b,x,y); ll t=x;x=y;y=t-a/b*y; } int main() { int i; ll x, y; scanf("%lld%lld%lld",&N,&A,&B); ji=1; for(i=1;i<=N;i++) { scanf("%lld%lld",&m[i],&a[i]); ji*=m[i]; } for(i=1;i<=N;i++) { M[i]=ji/m[i]; exgcd(M[i],m[i],x,y); t[i]=(x+m[i])%m[i]; ans=(ans+a[i]*t[i]%ji*M[i]%ji)%ji; } while(ans<A)ans+=ji; tmp=ans; while(ans<=B)ans+=ji,cnt++; if(cnt==0)tmp=0; printf("%lld\n%lld\n",cnt,tmp); return 0; }
相关文章推荐
- DOM之获取元素,获取和设置属性
- QNX 实时操作系统(Quick Unix)
- 第六周作业
- ajax返回数据类型 text json xml
- 实验七:Linux内核如何装载和启动一个可执行程序
- iOS开发--即时通讯
- Linux 性能检测工具Vmstat命令
- 面积公式小结
- JavaScript 触发事件列表
- ImageView 的 src 属性
- 在CentOS 7上构建RAID5、LVM和SAMBA服务器(4)——逻辑卷(组)扩容
- 状压dp
- 继承与派生定义及相关程序调用举例
- 1012. The Best Rank
- C语言中使用快速排序算法对元素排序的实例详解
- c++ tuple类型
- Qt5.2.0配置sdk
- 自定义view的measureSpec是谁的mode和size
- dapper 的学习
- BZOJ 4311 向量