1213 解的個數
已知整數x,y滿足如下面的條件:
ax+by+c = 0
p<=x<=q
r<=y<=s
求滿足這些條件的x,y的個數。
第一行有一個整數n（n<=10），表示有n個任務。n<=10
以下有n行，每行有7個整數，分別為：a,b,c,p,q,r,s。均不超過108。
共n行，第i行是第i個任務的解的個數。
2
2 3 -7 0 10 0 10
1 1 1 -10 10 -9 9
1
19
分類標籤 Tags 點此展開
神坑啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊
以下內容摘自某大神的題解
1.首先我們可以很直觀地看出來這是用擴充套件歐幾里得演算法解二元一次方程，但問題是我們所熟悉的擴歐能解的方程都是ax+by=gcd(a,b)形式的，而題目給出的是ax+by=-c形式。舉個例子：
2x+4y=18，首先我們可以換成x+2y=9形式。9不是2和4的最大公約數1，但9是1的倍數，所以如果我們解出一組x,y滿足x+2y=1，那麼x和y都乘上9/1就是原方程的一組解了。如果c/gcd(a,b)==0，那麼就沒有整數解。
2.如今我們得到了一組x,y，根據擴歐定理的後續內容，適合的解系一定是(x+bk,y-ak)，注意現在的a,b是簡化後的方程的係數（拿上面提到的例子講，現在a=1,b=2），列舉找在區間內的解的個數（組數）就好。
3.現在解決各種WA/TLE/RE問題：
（1）方程無解：c/gcd(a,b)==0，直接輸出0；
（2）區間不合法：題目中沒有保證區間左端點小於右端點，所以如果讀入的區間不合法，直接輸出0；
（3）a=0或b=0：
if((a==0)&&(y<r||y>s)) {printf("0\n");continue;}
if((b==0)&&(x<p||x>q)) {printf("0\n");continue;}
因為不管加多少，x/y還是原來的味道……但是如果不加特判可能會導致TLE（這個跟程式碼具體的寫法有關，我後面有用while迴圈，直接卡T了）
（4）a==0&&b==0：
RE的關鍵所在，因為gcd求出來是0……這個需要認真思考一下，如果c!=0，顯然方程不成立，無解；如果c==0，x和y就可以任意取了，由乘法原理可得解的個數就是兩個區間內部整數點的個數的乘積
if (c!=0)printf("0\n");
else
{
ll cnt=(q-p+1)*(s-r+1);
printf("%lld\n",cnt);
}
continue;
（5）記得要開long long
兩個神坑的資料點：
4
0 1 2 0 0 0 2
1 0 2 0 0 0 0
1 0 2 0 2 0 20
2 0 3 -10 10 -10 10
ans：0 0 0 0
4
0 0 0 -1 1 -1 1
0 0 0 1 -1 1 2
0 0 1 1 1 1 1
0 0 0 -3406792423987599 -23487749 23947250
ans：9 0 0 2753863780940000
1 #include<iostream> 2 #include<cstdio> 3 using namespace std; 4 long long int x,y; 5 long long int tot=0; 6 long long int a,b,c,p,q,r,s; 7 long long int gcd(long long int a,long long int b) 8 { 9 if(b==0) 10 return a; 11 else 12 return gcd(b,a%b); 13 } 14 long long int exgcd(long long int a,long long int b,long long int & x,long long int & y) 15 { 16 if(b==0) 17 { 18 x=1; 19 y=0; 20 return a; 21 } 22 long long int r=exgcd(b,a%b,x,y); 23 long long int tmp; 24 tmp=x; 25 x=y; 26 y=tmp-a/b*y; 27 return r; 28 } 29 int main() 30 { 31 32 int n; 33 scanf("%d",&n); 34 for(int i=1;i<=n;i++) 35 { 36 tot=0; 37 //scanf("%lld %lld %lld %lld %lld %lld %lld",&a,&b,&c,&p,&q,&r,&s); 38 cin>>a>>b>>c>>p>>q>>r>>s; 39 c=-c; 40 if((a==0)&&(y<r||y>s)) 41 { 42 printf("0\n"); 43 continue; 44 } 45 if((b==0)&&(x<p||x>q)) 46 { 47 printf("0\n"); 48 continue; 49 } 50 51 if(p>q||r>s) 52 { 53 cout<<0<<endl; 54 continue; 55 } 56 57 int gys=gcd(a,b); 58 if(gys==0) 59 { 60 if (c!=0) 61 { 62 printf("0\n"); 63 continue; 64 } 65 else 66 { 67 tot=(q-p+1)*(s-r+1); 68 printf("%lld\n",tot); 69 continue; 70 } 71 } 72 if(c%gys!=0) 73 { 74 cout<<0<<endl; 75 continue; 76 } 77 exgcd(a,b,x,y); 78 x=x*(c/gys); 79 y=y*(c/gys); 80 a=a/gys; 81 b=b/gys; 82 while(x>=p) 83 { 84 x=x-b; 85 y=y+a; 86 } 87 while(x<p&&b!=0) 88 { 89 x=x+b; 90 y=y-a; 91 } 92 while(x>=p&&x<=q&&y>=r&&y<=s) 93 { 94 tot++; 95 x=x+b; 96 y=y-a; 97 if(x<p||x>q||y<r||y>s) 98 break; 99 } 100 printf("%lld\n",tot); 101 } 102 103 return 0; 104 }