51nod 1352 集合计数 扩展欧几里得

Posted xjhz

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了51nod 1352 集合计数 扩展欧几里得相关的知识,希望对你有一定的参考价值。

基准时间限制:1 秒 空间限制:131072 KB 
 

给出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

AX+BY=N+1;

解方程;

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<queue>
#include<algorithm>
#include<stack>
#include<cstring>
#include<vector>
#include<list>
#include<bitset>
#include<set>
#include<map>
#include<time.h>
using namespace std;
#define LL long long
#define bug(x)  cout<<"bug"<<x<<endl;
const int N=1e5+10,M=1e6+10,inf=1e9+10,MOD=1e9+7;
const LL INF=1e18+10,mod=1e9+7;
const double eps=(1e-8),pi=(4*atan(1.0));

void exgcd(LL a, LL b, LL &x, LL &y)
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return;
    }
    exgcd(b, a % b, x, y);
    LL tmp = x;
    x = y;
    y = tmp - (a / b) * y;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        LL n,a,b;
        scanf("%lld%lld%lld",&n,&a,&b);
        LL x,y;
        LL gc=__gcd(a,b);
        a/=gc;
        b/=gc;
        exgcd(a,b,x,y);
        n++;
        if(n%gc==0)
        {
            n/=gc;x*=n;y*=n;
            LL xx=(x-(x/b)*b);
            LL yy=(y+(x/b)*a);
            if(xx<=0)xx+=b,yy-=a;
            //cout<<xx<<" "<<yy<<endl;
            if(xx>n||yy<0)printf("0\n");
            else printf("%lld\n",min((n-xx)/b+1,yy/a+((yy%a)?1:0)));
        }
        else printf("0\n");
    }
    return 0;
}
/*
p = p0 + b/Gcd(a, b) * t
q = q0 - a/Gcd(a, b) * t(其中t为任意整数)
*/

 

以上是关于51nod 1352 集合计数 扩展欧几里得的主要内容,如果未能解决你的问题,请参考以下文章

51Nod 1352 集合计数 扩展欧几里得

51 Nod 1352 集合计数(中国剩余定理+扩展欧几里得)

扩展欧几里德解的数量(51nod 1352)

51nod 1352 集合计数

51nod_1352_集合计数 (exgcd)

51nod 1256 扩展欧几里得