uva10375 Choose and divide
Posted invoid
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了uva10375 Choose and divide相关的知识,希望对你有一定的参考价值。
唯一分解定理。
挨个记录下每个质数的指数。
#include<cstdio> #include<algorithm> #include<cstring> #include<cmath> using namespace std; const int maxn = 100000 + 10; typedef long long LL; double ans; int p,q,s,t; bool mark[maxn]; int prime[maxn],cnt; int f[maxn]; void get_prime() { for(int i=2;i<=10000;i++) { if(!mark[i]) prime[++cnt]=i; for(int j=1;j<=cnt;j++) { int t=prime[j]*i; if(t>10000) break; mark[t]=1; if(i%prime[j]==0) break; } } } inline void work(int n,int d) { for(int i=1;i<=cnt;i++) if(!(n%prime[i])) { while(!(n%prime[i])) { n/=prime[i]; f[i]+=d; } if(n==1) break; } } inline void add(int m,int n,int d) { for(int i=m;i<=n;i++) work(i,d); } int main() { get_prime(); while(scanf("%d%d%d%d",&p,&q,&s,&t)==4) { memset(f,0,sizeof(f)); add(q+1,p,1); add(1,p-q,-1); add(t+1,s,-1); add(1,s-t,1); ans=1; //for(int i=1;i<=cnt;i++) if(f[i]) { // if(f[i]>0) for(int j=1;j<=f[i];j++) ans*=prime[i]; // else for(int j=1;j<=-f[i];j++) ans/=prime[i]; // } for(int i=1;i<=cnt;i++) if(f[i]) ans*=pow(prime[i],f[i]); printf("%.5lf\n",ans); } return 0; }
以上是关于uva10375 Choose and divide的主要内容,如果未能解决你的问题,请参考以下文章