Blog 5:筛法构造素数表，约分，数论基础

Uva 10375
汝佳还是很聪明的嘛，这个代码写的很好啊
啥时候能把书中的错误改好就行了
唉坑死读者
知识点1：
x以内的素数有
x/lnx的最高位加1个
例如
10000以内
用10000/log(10000)算出
1086
内部素数个数就设置为
2086
知识点2：
筛法构造素数表

int vis[10000];
	int n=10000;
    int m=sqrt(n+0.5);
    for(int i=2;i<=m;i++)
    if(!vis[i])
    for(int j=i*i;j<=n;j+=i)
    vis[j]=1;
	for(int i=2;i<=10000;i++)
    {
		if(!vis[i])
		{
		prime[primesize]=i;
        primesize++; 
		}
	}

知识点三
利用唯一分解定理
将每一个数分解分解结果保存在指数数组里

void add_integer(int n,int d)
{
	for(int i=0;i<primesize;i++)
	{
	while(n%prime[i]==0)
	{
	e[i]+=d;
	n/=prime[i];	
	}	
	if(n==1)
	break;
	}
}

总代码：

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
using namespace std;
int primesize=0;
int prime[2000];
int e[2000];
void add_integer(int n,int d)
{
	for(int i=0;i<primesize;i++)
	{
	while(n%prime[i]==0)
	{
	e[i]+=d;
	n/=prime[i];	
	}	
	if(n==1)
	break;
	}
}
void add_fac(int n,int d)
{
	for(int i=1;i<=n;i++)
	add_integer(i,d);
}
int main()
{
	int p,q,r,s;
    int vis[10000];
	int n=10000;
    int m=sqrt(n+0.5);
    for(int i=2;i<=m;i++)
    if(!vis[i])
    for(int j=i*i;j<=n;j+=i)
    vis[j]=1;
	for(int i=2;i<=10000;i++)
    {
		if(!vis[i])
		{
		prime[primesize]=i;
        primesize++; 
		}
	}
	//构造初始化素数表及初始化各个数据
	while(cin>>p>>q>>r>>s)
	{
	memset(e,0,sizeof(e));
	double ans=1;
	add_fac(p,1);
	add_fac(q,-1);
	add_fac(p-q,-1);
	add_fac(r,-1);
	add_fac(s,1);
	add_fac(r-s,1);
	for(int i=0;i<primesize;i++)
	ans*=pow((double)prime[i],(double)e[i]);
	printf("%.5lf\n",ans);		
    } 
	return 0;
}