WIKIOI 1213 解的个数 题解与分析

1213 解的个数

题目描述 Description

已知整数x,y满足如下面的条件:

 

ax+by+c = 0

p<=x<=q

r<=y<=s

 

求满足这些条件的x,y的个数。

输入描述 Input Description

第一行有一个整数n（n<=10），表示有n个任务。n<=10

以下有n行，每行有7个整数，分别为：a,b,c,p,q,r,s。均不超过10⁸。

输出描述 Output Description

共n行，第i行是第i个任务的解的个数。

样例输入 Sample Input

2

2 3 -7 0 10 0 10

1 1 1 -10 10 -9 9

样例输出 Sample Output

1

19

【算法】：拓展欧几里得

【分析】：

          求ax+by+c=的方程处在p<=x<=q和r<=y<=s中的解个数，加上各种特判。

          首先，转化方程ax+by=-c,求出ax+by=gcd(a,b)的基本解x，y<用拓展欧几里得>，当且仅当gcd(a,b)|-c时，

          因此特判一：当gcd(a,b)|-c不成立时，输出0.然后转化解x，y.对于方程ax+by=-c，其解x1，y1为x1=x/gcd(a,b)*-C,y1=y/gcd(a,b)*-C

         对于其他解，则由方程ax+by=-c→a(x+bt)+b(y-at)=-c,则x2=x1+bt,y2=y1-at <t∈(-∞,∞)且t∈z>，先将基础解的x降至p以下<p>x>，然后以x不越界为标准用while调整t求解即可。

         特判二：若a，b均为0且c为0，则原方程有无数解，然后在范围内的解有(q-p+1)*(s-r+1)。

         特判三：若a，b均为0且c不为0，则无解.

         特判四：除上面情况外若a为0,直接按by=-c来做.

         特判五：除上面情况外若b为0,直接按ax=-c来做。

         注意：①需考虑q>=p且s>=r，否则无解

                     ②：在进入过程前应将a，b除以gcd(a,b)，减小步长

【提交地址】：http://wikioi.com/problem/1213/

【代码】：

[cpp]  view plain copy

1. #include<stdio.h>  
2. #include<string.h>  
3. #include<stdlib.h>  
4. #include<algorithm>  
5. using namespace std;  
6. long long N,A,B,C,P,Q,R,S,X,Y;  
7. long long gcd(long long x,long long y){return x%y==0 ? y : gcd(y,x%y);};  
8. void exgcd(long long a,long long b)  
9. {  
10.       if(b==0)  
11.       {  
12.             X=1;  
13.             Y=0;  
14.             return;  
15.       }  
16.       exgcd(b,a%b);  
17.       long long tmp;  
18.       tmp=X;  
19.       X=Y;  
20.       Y=tmp-(a/b)*Y;  
21. }  
22. int main()  
23. {  
24.       //freopen("input.in","r",stdin);  
25.       //freopen("output.out","w",stdout);   
26.       scanf("%lld",&N);  
27.       for(int i=1;i<=N;i++)  
28.       {  
29.             scanf("%lld%lld%lld%lld%lld%lld%lld",&A,&B,&C,&P,&Q,&R,&S);  
30.             if(A==0 && B==0 && C!=0){printf("0\n");continue;}  
31.             if(A==0 && B==0 && C==0){printf("%lld\n",max((long long)0,Q-P+1)*max((long long)0,S-R+1));continue;}//当范围的第二个数-第一个数为负,则无解   
32.             if(A==0)  
33.             {  
34.                   long long sign=0;  
35.                   if(R<=-C/B && -C/B<=S && -C%B==0)   sign=1;  
36.                   printf("%lld\n",max((long long)0,Q-P+1)*sign);  
37.                   continue;  
38.             }  
39.             if(B==0)  
40.             {  
41.                   long long sign=0;  
42.                   if(P<=-C/A && -C/A<=Q && -C%A==0)   sign=1;  
43.                   printf("%lld\n",max((long long)0,S-R+1)*sign);  
44.                   continue;  
45.             }              
46.             exgcd(A,B);  
47.             long long ans=0,GCD=gcd(A,B);  
48.               
49.             if(-C%GCD!=0){printf("0\n");continue;}  
50.                
51.             X=X*-C/GCD;  
52.             Y=Y*-C/GCD;  
53.             A/=GCD;//减小步长   
54.             B/=GCD;  
55.             if(B<0)   B=-B,A=-A;  
56.             while(P<=X)  
57.             {  
58.                   X-=B;  
59.                   Y+=A;  
60.             }  
61.             while(X<=Q)  
62.             {       
63.                   if(P<=X && X<=Q && R<=Y && Y<=S)   ans++;  
64.                   X+=B;  
65.                   Y-=A;  
66.             }  
67.             printf("%lld\n",ans);  
68.       }  
69.       return 0;  
70. }