前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >洛谷P1072 Hankson 的趣味题(数学)

洛谷P1072 Hankson 的趣味题(数学)

作者头像
attack
发布2018-11-29 17:16:01
3770
发布2018-11-29 17:16:01
举报

题意

题目链接

Sol

充满套路的数学题。。

如果你学过莫比乌斯反演的话不难得到两个等式

\[gcd(\frac{x}{a_1}, \frac{a_0}{a_1}) = 1\]

\[gcd(\frac{b_1}{b_0}, \frac{b_1}{x}) = 1\]

然后枚举\(b_1\)的约数就做完了。。

代码语言:javascript
复制
// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int MAXN = 1e6;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int T, a0, a1, b0, b1, ans;
int gcd(int a, int b) {
    return b == 0 ? a : gcd(b, a % b);
}
void check(int x) {
    if(x % a1) return ;
    ans += (gcd(x / a1, a0 / a1) == 1 && gcd(b1 / b0, b1 / x) == 1);
}
int main() {
    T = read();
    while(T--) {
        a0 = read(), a1 = read(), b0 = read(), b1 = read(); ans = 0;
        for(int x = 1; x * x <= b1; x++) {
            if(b1 % x == 0) {
                check(x);
                if(b1 != x) check(b1 / x);
            }
        }
        cout << ans << endl;
    }
    return 0; 
}
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2018-10-27 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 题意
  • Sol
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档