代码
ull ans[MAXN];
void db()
{
memset(ans,0,sizeof(ans));
ans[1]=1;
int i,j;
for(i=2;i<=MAXN;i++)
{
if(!ans[i])
{
for(j=i;j<=MAXN;j+=i)
{
if(!ans[j])
ans[j]=j;
ans[j]=ans[j]/i*(i-1);
}
}
}
}