bzoj2190 [SDOI2008]仪仗队(欧拉函数)

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj2190 [SDOI2008]仪仗队(欧拉函数)相关的知识,希望对你有一定的参考价值。

2190: [SDOI2008]仪仗队

Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 3203  Solved: 2062
[Submit][Status][Discuss]

Description

  作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。    技术分享   现在,C君希望你告诉他队伍整齐时能看到的学生人数。

Input

  共一个数N。

Output

  共一个数,即C君应看到的学生人数。

Sample Input

  4

Sample Output

  9


HINT

【数据规模和约定】   对于 100% 的数据,1 ≤ N ≤ 40000

 

/*
可以将最左下的点标为(0,0), 那么显然如果存在一个点(x,y), 且有gcd(x; y) = k (k .= 1),
那么点(x/k,y/k) .定会将点(x,y) 挡住.
如果k = 1, 那么点(x, y) .定会被看到
所以点(x, y)被看到的充分必要条件是Gcd(x, y) == 1;
我们考察矩阵的下三角形,考察他的每一行,可以发现,这一行能够被看到的点的数目就是phi(x)
答案不难发现是∑(phi[x])*2+1(容斥原理)
*/ 

#include<iostream>
#include<cstdio>
#include<cstring>

#define N 50007

using namespace std;
int prime[N];
bool not_prime[N];
int n,ans,tot=1;

void pr_()
{
    prime[1]=2;not_prime[1]=true;
    for(int i=2;i<=N;i++)
    {
        if(!not_prime[i]) prime[++tot]=i;
        for(int j=1;j<=tot && i*prime[j]<=N;j++)
        {
            not_prime[i*prime[j]]=true;
            if(i%prime[j]) break;
        }
    }
}

int get_phi(int x)
{
    int ret=1;
    for(int i=1;prime[i]*prime[i]<=x;i++)
    {
        if(x%prime[i]==0)
        {
            ret*=prime[i]-1;x/=prime[i];
            while(x%prime[i]==0) x/=prime[i],ret*=prime[i];
        }
    }
    if(x>1) ret*=x-1;
    return ret;
}

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

 

以上是关于bzoj2190 [SDOI2008]仪仗队(欧拉函数)的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ-2190: [SDOI2008]仪仗队 (欧拉函数)

BZOJ 2190SDOI 2008仪仗队 欧拉筛

bzoj 2190 [SDOI2008]仪仗队(欧拉函数)

欧拉函数 BZOJ2190 [SDOI2008] 仪仗队

BZOJ2190 [SDOI2008]仪仗队(欧拉函数)

BZOJ2190 SDOI2008 仪仗队 gcd,欧拉函数