题目连接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=1577
这个题,我先开始用暴力,当然不是纯暴力,当场就TLE了。。
PS:后来我A了这个题之后又在自己学校的OJ加这个题,自己造的数据不管有多强大,有一种暴力始终可以过。。。后来我加了好多万组数据,最后连我在ZOJ上提交AC的代码都要250MS过的时候,那个暴力最后还是可以在800MS左右过掉题目。。。但是在ZOJ上,我的程序0MS,而暴力的那个是TLE。。可见ZOJ的数据真的做得很好。。。
先说说解法吧。
假设g=gcd(p,q),l=lcm(p,q),key=l/g
我们可以想象一下:
那个key是什么?想清楚这个问题就很简单了。
l是p和q公共的素因子的并集,而g则是交集(仔细理解这句话!!)
于是key的意思就是l刚好比g多的一部分素因子
于是我们可以选择这些素因子和g组合得到了p
另一部分和g组合得到了q
假设key的素因子个数(重复的不算)有n个
那么选择方案就有:C(n,0)+C(n,1)+C(n,2)+...+C(n,n)
熟悉组合学公式的童鞋肯定很快就可以反应过来
上面那个式子的和为2^n,所以剩下的工作就非常的简单了
我的代码:
#include<stdio.h> #include<algorithm> using namespace std; int prime[1100]; bool flag[1100]; void init() { int i,j,num=0; for(i=2;i<1100;i++) { if(!flag[i]) { prime[num++]=i; for(j=i*i;j<1100;j=j+i) flag[j]=true; } } } void swap(int &a,int &b) { int temp; temp=a; a=b; b=temp; } int power(int n) { int i,s=1; for(i=1;i<=n;i++) s=s*2; return s; } void slove(int n) { int i,ans=0; for(i=0;prime[i]*prime[i]<=n;i++) { if(n%prime[i]==0) { ans++; n=n/prime[i]; while(n%prime[i]==0) n=n/prime[i]; } if(n==1) break; } if(n>1) ans++; printf("%d/n",power(ans)); } int main() { int x,y,n; init(); while(scanf("%d%d",&x,&y)!=EOF) { if(y%x!=0) { printf("0/n"); continue; } n=y/x; slove(n); } return 0; }