bzoj 2190 [SDOI2008]仪仗队（欧拉函数）

 

【题目链接】

 

    http://www.lydsy.com/JudgeOnline/problem.php?id=2190

 

【题意】

   

    n*n的正方形，在(0,0)格点可以看到的格子数目。

 

【思路】

   

    预处理出欧拉函数。

    (x,y)=1,1<=y<=n,x<y的数对为t=sigma{ phi(i) } 1<=i<=n，则答案为2*t+1。

    当然也可以用莫比乌斯反演求，懒得写了 (〜￣△￣)〜

 

【代码】

 

#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;

const int N = 5e4+10;

int n;
int phi[N];

void get_phi(int n)
{
    phi[1]=1;
    for(int i=2;i<=n;i++) if(!phi[i]) 
    {
        for(int j=i;j<=n;j+=i) 
        {
            if(!phi[j]) phi[j]=j;
            phi[j]=phi[j]/i*(i-1);
        }
    }
}

int main()
{
    scanf("%d",&n);
    get_phi(n);
    int ans=0;
    for(int i=1;i<n;i++) ans+=phi[i];
    printf("%d",ans*2+1);
    return 0;
}