Disgruntled Judge UVa 10375

最新推荐文章于 2023-07-30 11:33:53 发布

AYZXX

最新推荐文章于 2023-07-30 11:33:53 发布

阅读量249

点赞数

分类专栏：数学数论

本文链接：https://blog.csdn.net/AYZXX/article/details/78745687

版权

数学同时被 2 个专栏收录

32 篇文章 0 订阅

订阅专栏

数论

6 篇文章 0 订阅

订阅专栏

题目链接

思路：这道题是算法竞赛紫书第十章第316页的一个习题，上面有主要的代码；

方法：因为含有阶乘，直接运算，数太大，一定会溢出的，而书上的方法，避免了这种溢出，而是采用分解的方法，利用数论中合数可以分解若干个素数相乘的形式，相乘或者相除时，只需在对应素数的指数相加减，即可。我们可以另开一个数组，用来对应存素数的指数。（代码简单易懂）。

#include <iostream>
#include <string.h>
#include <math.h>
#include <vector>
#include <stdio.h>
#include <algorithm>
using namespace std;
int e[100000];
bool book[100000];
vector<int>prime;
void Prime()
{
    for(int i=2;i<=10000;i++)
    {
        if(book[i]==false)prime.push_back(i);
        for(int j=i*i;j<=10000;j+=i)
            if(book[j]==false)
            book[j]=true;
    }
}
void Add_factorial(int n,int d)
{
    for(int i=0; i<prime.size()&&n>1; i++)
        while(n%prime[i]==0)
        {
            e[i]+=d;
            n=n/prime[i];
        }
}
void add_factorial(int m,int k)
{
    for(int i=1;i<=m;i++)
        Add_factorial(i,k);
    return ;
}
int main()
{
    int p,q,r,s;
    memset(book,false,sizeof(book));
    Prime();
    while(cin>>p>>q>>r>>s)
    {
        memset(e,0,sizeof(e));
        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);
        double cont=1.0;
        for(int i=0;i<prime.size(); i++)
            cont*=pow(prime[i],e[i]);
        printf("%.5lf\n",cont);
    }
    return 0;
}