UVA 10375 唯一分解定理

Choose and divide Time Limit: 3000MS     64bit IO Format: %lld & %llu Submit Status uDebug Description

Problem descriptions: System Crawler 2016-08-19

#include"iostream"
#include"cstring"
#include"math.h"
#include"cstdio"
using namespace std;
int prime[1300];
int e[1300];
int visit[10005];
int ps;
int list()
{
    int cnt=0;
    memset(visit,0,sizeof(visit));
    for(int i=2;i*i<=10003;i++)
    if(!visit[i])
    {
        for(int j=i*i;j<=10003;j+=i)
            if(!visit[j]) visit[j]=1;
    }
    for(int i=2;i<=10003;i++)
        if(!visit[i])
        prime[cnt++]=i;
    return cnt;


}
int fuck(int q,int state) //state=1 表示在分子 -1表示在分母
{
    int i,j;
    for(i=1;i<=q;i++)
    {
        for(j=0;j<ps;j++)
        {
            if(i/prime[j]==0) break;  //当这个素数已经大于 这一项 就可以跳出循环了
            if(i%prime[j]==0) 
            {
                //e[j]+=state;
                int t=i;
                while(t%prime[j]==0) //这个素数的次数 由他 被余的次数决定
                {
                    e[j]+=state;
                    t/=prime[j];
                }
            }

        }
    }
}
int main()
{
    int p,q,r,s;ps=list();
    while(cin>>p>>q>>r>>s)
    {
        memset(e,0,sizeof(e));
        fuck(p,1);fuck(s,1);fuck(r-s,1);  //分子部分 都是乘
        fuck(q,-1);fuck(r,-1);fuck(p-q,-1);//分母部分都是除
        double ans=1;int i;
        for(i=0;i<ps;i++)
            ans*=pow(prime[i],e[i]);
        printf("%.5f\n",ans);
    }
    return 0;
}