qdu 2017级纳新题(扩展欧几里得)

在你面前撒个娇 哎呦喵喵喵喵喵

Description

我们一起学猫叫

一起喵喵喵喵喵

在你面前撒个娇

哎呦喵喵喵喵喵

我的心脏砰砰跳

迷恋上你的坏笑

你不说爱我我就喵喵喵

每当xjy和hqy一起唱起这首歌时,就会吸引好多猫群来听歌,这天他们又吸引来n群猫。

每个猫群里都分为雄猫和雌猫,第一个猫群里有1只雄猫和n只雌猫,第二个猫群里有2只雄猫和n-1只雌猫....第n个猫群里有n只雄猫和1只雌猫。

xjy和hqy都很喜欢猫,他们对猫分别有好感度a和b,当一个猫群的雄猫数是a的倍数时,这个猫群会对xjy产生好感,同样的,当一个猫群的雌猫数是b的倍数时,这个猫群会对hqy产生好感。

演唱结束后,他们想带一些猫群回他们自己的小窝,但是只有同时对xjy和hqy都产生好感的猫群才能被带走。

他们想知道能带走多少猫群,但xjy在hqy面前通常比较笨,所有只能求救于你了。

Input

第一行一个整数T,(1<=T<=50000),表示有T组测试样例。

之后T行,每行三个整数分别表示n,a,b(1<=n,a,b<=2147483647)

Output

共输出T行,每行对应一个测试数据,表示能带走的猫群数量。

Sample Input 1 

2
5 2 4
10 2 3

Sample Output 1

1
2

Hint

对于第二个测试样例,可以带走的是2只雄猫9只雌猫的猫群和8只雄猫3只雌猫的猫群。

题意:表示选拔的时候死活一直错很迷,后来看了学长的代码,换了一种写法就对了,好迷呀,中文题,不解释啦~~

题解:这个题就是求方程 ax+by=n+1 的正整数解的个数,先跑一遍扩展欧几里得求方程 ax+by=gcd(a,b)的解,然后根据通解公式 x=x0*(n+1)/gcd+b*t/gcd  ,   y=y0*(n+1)/gcd-a*t/gcd(t 为整数) 可以计算出 x,y 的最小正整数解  x=(x0%(b/gcd)+(b/gcd) )%(b/gcd)  ( x0 可能为负数,所以要加一个(b/gcd) )   y=(y0%(a/gcd)+(a/gcd) )%(a/gcd)  ( y0 可能为负数,所以要加一个(a/gcd)) 不要问我为什么,记住就好了!然后就是,x的最小值对应y的最大值,根据 by=(n+1)-(a*x)求出y的最大值,然后根据通解公式可以知道y最大-y最小=y0*(n+1)/gcd-(y0*(n+1)/gcd-a*t最大/gcd=a*t最大/gcd 求出t最大,然后结果加一就好了。详细请看代码:

#include <iostream>
#include <cstring>
#include <queue>
using namespace std;
typedef long long ll;
ll n,a,b,x,y;
ll ex_gcd(ll a, ll b, ll &x, ll &y){//扩展欧几里得模板
	if(b==0){
		x=1,y=0;
		return a;
	}
	ll p=ex_gcd(b,a%b,y,x);
	y-=a/b*x;
	return p;
}
int main(){
	int t;
	cin >> t;
	while(t--){
		cin >> n >> a >> b;
		ll gcd=ex_gcd(a,b,x,y);
		ll ans=0;
		if((n+1)%gcd==0){
			ll k=(n+1)/gcd;
			x*=k;//别忘记乘以(n+1)/gcd,因为所求x为ax+by=gcd(a,b)的解
			y*=k;//别忘记乘以(n+1)/gcd,因为所求y为ax+by=gcd(a,b)的解
			ll r1=b/gcd;		
			x=(x%r1+r1)%r1;
			if(x==0) x=r1;//求出非0的x的最小解的	
			ll re=n+1-x*a;
			if(re<=0) ans=0;
			else{
				ll yy=re/b;//y的最大解
				ll r2=a/gcd;
				y=(y%r2+r2)%r2;//求出非0的y的最小解
			    if(y==0) y=r2;
				ans=(yy-y)/r2+1;//最大减最小除以a/gcd+1别忘记加一
			}		
		}
		cout << ans << endl;
	}
	return 0;
}

猜你喜欢

转载自blog.csdn.net/lgz0921/article/details/82108921