BZOJ大视野 2190: [SDOI2008]仪仗队

最新推荐文章于 2018-08-09 17:42:00 发布

L_0_Forever_LF

最新推荐文章于 2018-08-09 17:42:00 发布

阅读量895

点赞数

分类专栏： BZOJ 数学文章标签： BZOJ大视野数学

本文链接：https://blog.csdn.net/L_0_Forever_LF/article/details/43601125

版权

BZOJ 同时被 2 个专栏收录

374 篇文章 1 订阅

订阅专栏

数学

32 篇文章 0 订阅

订阅专栏

经过观察，发现一个规律，n>2时

无视原图第一列第n行，剩下一个(n-1)*(n-1)的矩阵，m=n-1

剩下的图中，1——m列，m到1行，每列、行可见的数量为 <=i 且与 i 互质的数的数量

所以用到欧拉函数

code：

/**************************************************************
    Problem: 2190
    Language: C++
    Result: Accepted
    Time:40 ms
    Memory:11028 kb
****************************************************************/
 
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int prime[1110000],pr,phi[1110000];
bool v[1110000];
int n;
 
void getprime()
{
    memset( v,true,sizeof v );
    pr = 0;
    int i,j;
    for( i=2;i<=n;i++ )
    {
        if( v[i] )
        {
            prime[++pr]=i;
            phi[i] = i-1;
        }
        for( j=1; j<=pr && prime[j]*i<=n; j++ )
        {
            v[ prime[j]*i ] = false;
            if( i%prime[j] == 0 )
            {
                phi[ i*prime[j] ] = phi[i] * prime[j];
                break;
            }
            else
                phi[ i*prime[j] ] = phi[i] * phi[prime[j]];
        }
    }
}
int main()
{
    int i,k;
    n=50000;
    getprime();
     
    scanf("%d",&n);
      
    if( n==1 )printf("0\n");
    else if( n==2 ) printf("1\n");
    else
    {
        n--;
        k=3;
        for(i=2;i<=n;i++)
            k+= 2*phi[i];
        printf("%d\n",k);
    }
     
    return 0;
}