bzoj2190 [SDOI2008]仪仗队

题目链接

如果只考虑y=x直线的上方

每一种斜率的直线只能打到一个人

转化为求一个数x(0~n-1)做分母时，有多少个数y与x互质(欧拉函数)

为什我我要做这种水题？ 因为我太弱了

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<string>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define rre(i,r,l) for(int i=(r);i>=(l);i--)
#define re(i,l,r) for(int i=(l);i<=(r);i++)
#define Clear(a,b) memset(a,b,sizeof(a))
#define inout(x) printf("%d",(x))
#define douin(x) scanf("%lf",&x)
#define strin(x) scanf("%s",(x))
#define LLin(x) scanf("%lld",&x)
#define op operator
#define CSC main
typedef unsigned long long ULL;
typedef const int cint;
typedef long long LL;
using namespace std;
void inin(int &ret)
{
    ret=0;int f=0;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=1;ch=getchar();}
    while(ch>='0'&&ch<='9')ret*=10,ret+=ch-'0',ch=getchar();
    ret=f?-ret:ret;
}
int n;
bool bo[1000010];
int p[100010],x,phi[100010];
void getphi()
{
    phi[1]=1;
    re(i,2,n)
    {
        if(!bo[i])p[++x]=i,phi[i]=i-1;
        for(int j=1;j<=x&&i*p[j]<=n;j++)
        {
            bo[i*p[j]]=1;
            if(i%p[j]==0)
            {
                phi[i*p[j]]=phi[i]*p[j];
                break;
            }
            else phi[i*p[j]]=phi[i]*(p[j]-1);
        }
    }
}
int main()
{
    inin(n);
    getphi();int ans=0;
    re(i,1,n-1)ans+=phi[i];
    cout<<(ans*2+1);
     return 0;
}