51Nod-集合计数（拓展欧几里得+中国剩余定理）



1352 集合计数


基准时间限制：1 秒 空间限制：131072 KB 分值: 20 难度：3级算法题


 收藏


 关注

给出N个固定集合{1，N},{2,N-1},{3,N-2},...,{N-1,2},{N,1}.求出有多少个集合满足：第一个元素是A的倍数且第二个元素是B的倍数。
提示：
对于第二组测试数据，集合分别是：{1,10},{2,9},{3,8},{4,7},{5,6},{6,5},{7,4},{8,3},{9,2},{10,1}.满足条件的是第2个和第8个。

Input

第1行：1个整数T（1<=T<=50000)，表示有多少组测试数据。
第2 - T+1行：每行三个整数N，A,B（1<=N，A，B<=2147483647)

Output

对于每组测试数据输出一个数表示满足条件的集合的数量，占一行。

Input示例

2
5 2 4
10 2 3

Output示例

1
2

System Message (题目提供者)

题目考点很容易看出，就是求解ax+by=n+1中解的个数

这里不能完全明白的地方是为什么除以LCM便是剩下可能的情况（巨巨们说是中国剩余定理）

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<limits.h>
#include<string>
#include<queue>
#include<vector>
#include<stack>
#include<math.h>
#include<map>
#include<algorithm>
using namespace std;
typedef long long ll;
#define maxn 1000005
#define Mod 1000000007
const int inf=1e10;
ll gcd(ll x,ll y)
{
if(y==0)return x;
return gcd(y,x%y);
}
void ex_gcd(ll a,ll b,ll &x,ll &y)
{
if(b==0)
{
x=1;
y=0;
return;
}
ll x1,y1;
ex_gcd(b,a%b,x1,y1);
x=y1;
y=x1-a/b*y1;
}
int  main()
{
ll n,a,b,i,j,T;
scanf("%lld",&T);
while(T--)
{
ll x,y;
scanf("%lld%lld%lld",&n,&a,&b);
n++;
ll temp=gcd(a,b),LCM=a/temp*b;
if(n%temp)
{
printf("0\n");
continue;
}
ll dn=n/temp;
ll da=a/temp;
ll db=b/temp;
ex_gcd(da,db,x,y);
x*=dn;
x=(x%db+db)%db;//求最小的非负x的解
if(x==0)
x+=db;//保证x是整数
ll ans=n-1-x*a;
//printf("%lld",a);
if(ans<0)printf("0\n");
else  printf("%lld\n",1+ans/LCM);
}
}