题解 P1072 【Hankson 的趣味题】

zzlzk

2017-08-26 21:09:20

Solution

- LaTex写公式有点麻烦,所以我用以前写好的代替了 - 想看原版可以戳[这里](http://blog.csdn.net/nuclearsubmarines/article/details/77603154) - 首先来分析一下这个题目 ![](https://i.loli.net/2017/08/26/59a16f6ad2018.png) 证明: ![](https://i.loli.net/2017/08/26/59a16fb56c248.png) - 把上面的结论推广一下,得到结论$P$ >对于两个正整数$a,b$,设$gcd(a,b)=k$,则存在$gcd(a/k,b/k)=1$ - 应用结论$P$ ![](https://i.loli.net/2017/08/26/59a170dc98ec0.png) - 整理一下式子 ![](https://i.loli.net/2017/08/26/59a1711685e4f.png) 用心体会这两个式子,你会发现$x$是$a_1$的整数倍而且是$b_1$的因子 ~~好像这个由gcd和lcm也可以得到?~~嗯,就这样 于是得到一种解题思路 >$\sqrt b_1$枚举$b_1$的因子(也就是$x$),如果这个数是$a_1$的整数倍并且满足那两个式子,则$ans++$ - code ```cpp #include<cstdio> using namespace std; int gcd(int a,int b) { return b==0?a:gcd(b,a%b); } int main() { int T; scanf("%d",&T); while(T--) { int a0,a1,b0,b1; scanf("%d%d%d%d",&a0,&a1,&b0,&b1); int p=a0/a1,q=b1/b0,ans=0; for(int x=1;x*x<=b1;x++) if(b1%x==0){ if(x%a1==0&&gcd(x/a1,p)==1&&gcd(q,b1/x)==1) ans++; int y=b1/x;//得到另一个因子 if(x==y) continue; if(y%a1==0&&gcd(y/a1,p)==1&&gcd(q,b1/y)==1) ans++; } printf("%d\n",ans); } return 0; } ```