Description
大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。
Input
第一行为两个整数T,R。R<=10^9+10,T<=10000,表示该组中测试数据数目,R为模后面T行,每行一对整数N,M,见题目描述 m<=n
Output
共T行,对于每一对N,M,输出1至N!中与M!素质的数的数量对R取模后的值
Sample Input
1 11
4 2
4 2
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
数据范围:
对于100%的数据,1 < = N , M < = 10000000
又!是!数!论!题!
看这个吧。。http://www.cnblogs.com/zyfzyf/p/3997986.html
#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
int mu[50001],prime[50001];
bool check[50001];
int tot;
int sum[50001];
inline int min(int x,int y)
{
if(x<y)
return x;
return y;
}
inline void findmu()
{
memset(check,false,sizeof(check));
mu[1]=1;
int i,j;
for(i=2;i<=50000;i++)
{
if(!check[i])
{
tot++;
prime[tot]=i;
mu[i]=-1;
}
for(j=1;j<=tot;j++)
{
if(i*prime[j]>50000)
break;
check[i*prime[j]]=true;
if(i%prime[j]==0)
{
mu[i*prime[j]]=0;
break;
}
else
mu[i*prime[j]]=-mu[i];
}
}
for(i=1;i<=50000;i++)
sum[i]=sum[i-1]+mu[i];
}
int main()
{
// freopen("zap.in","r",stdin);
// freopen("zap.out","w",stdout);
findmu();
int b,d,k;
int T;
scanf("%d",&T);
while(T!=0)
{
T--;
scanf("%d%d%d",&b,&d,&k);
if(b>d)
{
int tt=b;
b=d;
d=tt;
}
int ans=0;
int i;
int la=0;
b=b/k;
d=d/k;
for(i=1;i<=b;i=la+1)
{
la=min(b/(b/i),d/(d/i));
ans+=(b/i)*(d/i)*(sum[la]-sum[i-1]);
}
printf("%d\n",ans);
}
return 0;
}