【hdu1573-X问题】拓展欧几里得-同余方程组
http://acm.hdu.edu.cn/showproblem.php?pid=1573
求小于等于N的正整数中有多少个X满足:
X mod a0 = b0
X mod a1 = b1
……
X mod ai = bi
(0<ai<=10)
输入:第一行为一个正整数T,表示有T组测试数据。每组测试数据的第一行为两个正整数N,M (0 < N <= 1000,000,000 , 0 < M <= 10),表示X小于等于N,数组a和b中各有M个元素。接下来两行,每行各有M个正整数,分别为a和b中的元素。
输出:对应每一组输入,在独立一行中输出一个正整数,表示满足条件的X的个数。
Sample Input
3
10 3
1 2 3
0 1 2
100 7
3 4 5 6 7 8 9
1 2 3 4 5 6 7
10000 10
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 7 8 9
Sample Output
1
0
3
题解:同余方程组的裸题。尤其要注意的是所求的是正整数,也就是0的时候要特殊处理。
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
#include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> using namespace std; typedef long long LL; const LL M=20; LL a[M],b[M]; LL exgcd(LL u,LL v,LL &x,LL &y) { if(v==0){ x=1;y=0;return u;} else { LL tx,ty; LL d=exgcd(v,u%v,tx,ty); x=ty; y=tx-(u/v)*ty; return d; } } int main() { // freopen("a.in","r",stdin); // freopen("a.out","w",stdout); LL T; scanf("%I64d",&T); while(T--) { LL n,m; bool bk=1; scanf("%I64d%I64d",&n,&m); for(LL i=1;i<=m;i++) scanf("%I64d",&a[i]); for(LL i=1;i<=m;i++) scanf("%I64d",&b[i]); LL A,B,C,g,x,tx,ty,a1,b1; a1=a[1],b1=b[1]; for(LL i=2;i<=m;i++) { A=a1;B=a[i];C=b[i]-b1; g=exgcd(A,B,tx,ty); if(C%g) {bk=0;break;} x=((tx*C/g)%(B/g)+(B/g))%(B/g); b1=a1*x+b1; a1=a1/g*a[i]; } if(bk) { /* 得出不定方程a1x+b1=P后,x=0时b1=P。 因为x可以等于0:ans=(n-b1)/a1+1; 若P=0,则b1=x=0,不满足正整数,所以减去。 */ LL ans=0; if(n>=b1) ans=(n-b1)/a1+1; if(ans && b1==0) ans--; printf("%I64d\n",ans); } else printf("0\n"); } return 0; }