BZOJ2190: [SDOI2008]仪仗队

有两种做法一个是欧拉函数的前缀和另一个是莫比乌斯反演
貌似数据水能让莫比乌斯反演过去
我用的是莫比乌斯反演

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int maxn=40001;
int mu[maxn],pre[maxn],prime[maxn];
bool check[maxn];
int tot;
int main()
{
     int i,j,k;
     mu[1]=1;
     for(i=2;i<maxn;i++)    
      {
        if(!check[i])
             prime[++tot]=i,mu[i]=-1;
        for(j=1;j<=tot;j++)
          {
              k=prime[j]*i;
              if(k>=maxn)break;
              if(i%prime[j]==0){check[k]=true,mu[k]=0;break;}
              mu[k]=-mu[i];check[k]=true;
          }
      }
    for(i=1;i<maxn;i++)pre[i]=pre[i-1]+mu[i];
    int n;
    scanf("%d",&n);
    n--;
    int T=1,Div1,Div2,l;
    int ans=0;  
    while(T<=n)
    {
        Div1=n/T;
        l=n/Div1+1;
        ans+=Div1*Div1*(pre[l-1]-pre[T-1]);
        T=l;
    }
    ans+=2;
    printf("%d\n",ans);
    return 0;
}