国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

51Nod - 1352 拓展gcd

2019-11-10 20:19:35
字體:
來源:轉載
供稿:網友

題意:

給出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示例
25 2 410 2 3Output示例
12

思路:

根據題意可以列出一個不定方程a*x+b*y=n+1,利用拓展gcd求一組可行解。要求出大于0的最小x,直接用所求出的一組可行解(x%(b/g) + (b/g))%(b/g)即可(其中g是a和b的最大公約數)。這樣就求得了最小的a*x,那么剩下的就是n-a*x個數,x的通解形式是x+k*(b/g),那么a*x的通解就是a*x+k*(a*b/g)==a*x+k*lcm,所以直接拿n-a*x除以lcm(a,b)再算上一開始的那個最小的a*x就是最終答案。

代碼:

#include <bits/stdc++.h>using namespace std;typedef long long ll;ll extgcd(ll a, ll b, ll &x, ll &y) {    ll d = a;    if (b != 0) {        d = extgcd(b, a % b, y, x);        y -= (a / b) * x;    }    else {        x = 1; y = 0;    }    return d;}int main() {    int T;    scanf("%d", &T);    while (T--) {        ll n, a, b, x, y;        scanf("%I64d%I64d%I64d", &n, &a, &b);        ++n;        ll g = extgcd(a, b, x, y);        if (n % g) {            puts("0");            continue;        }        x *= n / g;        ll bb = b / g;        x = (x % bb + bb) % bb;        if (x == 0) x += bb;        if (a * x >= n) {            puts("0");            continue;        }        else {            ll lcm = a / g * b;            PRintf("%I64d/n", (n - 1 - x * a) / lcm + 1);        }    }    return 0;}
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 峨眉山市| 商洛市| 二连浩特市| 鄂伦春自治旗| 宜兴市| 瑞安市| 武乡县| 夏河县| 遂昌县| 施秉县| 叶城县| 五台县| 淮安市| 洪泽县| 靖江市| 五大连池市| 乃东县| 黄梅县| 五常市| 兴隆县| 绿春县| 海门市| 碌曲县| 军事| 崇左市| 灯塔市| 临夏市| 莱阳市| 平武县| 土默特左旗| 长岛县| 满城县| 西乡县| 竹北市| 乐都县| 筠连县| 镶黄旗| 茌平县| 都匀市| 河间市| 峨山|