hdu2841 筛选素因子 容斥定理

hdu2841

#include<stdio.h>
#include<string.h>
#define max 100001

int Prim[max],num[max][20];
void init()
{
	int i,j;
	memset(Prim,0,sizeof(Prim));
    for(i=1; i<=100000; i++)	num[i][0]=0;
	for(i=2; i<=100000; i++)
	if(Prim[i]==0)
	{
		num[i][1]=i;
		num[i][0]++;
		for(j=i*2; j<=100000; j+=i)
		{
			num[j][++num[j][0]]=i;
			Prim[j]=1;
		}
	}
}

int Noprime(int m,int n,int x)
{
	int i,ret=0;
	for(i=x;i<=num[m][0];i++)
		ret+=n/num[m][i]-Noprime(m,n/num[m][i],i+1);
	return ret;
}

int main()
{
	int i,T,m,n;
	__int64 k;
	init();
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&m,&n);
		k=n;
		for(i=2;i<=m;i++)
			k+=n-Noprime(i,n,1);
		printf("%I64d\n",k);
	}
	return 0;
}