UVa 10375 唯一分解定理

来源：互联网发布：手机淘宝叫快递编辑：程序博客网时间：2024/06/05 16:24

#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>#include<cstdlib>#include<vector>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fod(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=1e5+10;int prime[N],e[N],p,q,r,s,tot=0;void Euler_prime() {    bool vis[N]; memset(vis,0,sizeof(vis));    for(int i=2;i<=N;i++) {        if(!vis[N]) prime[++tot]=i;        for(int j=1;j<=tot;j++) {            if(i*prime[j]>N) break;            vis[i*prime[j]]=1;            if(i%prime[j]==0) break;        }    }}void add_integer(int n,int d) {    for(int i=1;i<=tot;i++) {        while(n%prime[i]==0) {            e[i]+=d;            n/=prime[i];        }        if(n==1) break;    }}void add_factorial(int n,int d) {    for(int i=1;i<=n;i++)        add_integer(i,d);}void init() {    memset(e,0,sizeof(e));}void solve() {    add_factorial(p,1);    add_factorial(q,-1);    add_factorial(p-q,-1);    add_factorial(r,-1);    add_factorial(s,1);    add_factorial(r-s,1);}int main() {    Euler_prime();    while(~scanf("%d%d%d%d",&p,&q,&r,&s)) {        init();        solve();        double ans=1.00000;        for(int i=1;i<=tot;i++)            ans*=pow(prime[i],e[i]);        printf("%.5lf\n",ans);    }    return 0;}

欧拉筛在n<=1e5时跑的没有埃式筛快~~orz

阅读全文

0 0