bzoj 2186: [Sdoi2008]沙拉公主的困惑

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

Sample Output

1

数据范围：
对于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;
}