2190: [SDOI2008]儀仗隊

Time Limit:?10 Sec??Memory Limit:?259 MB
[Submit][Status][Discuss]

Description

　　作為體育委員，C君負責(zé)這次運動會儀仗隊的訓(xùn)練。儀仗隊是由學(xué)生組成的N * N的方陣，為了保證隊伍在行進中整齊劃一，C君會跟在儀仗隊的左后方，根據(jù)其視線所及的學(xué)生人數(shù)來判斷隊伍是否整齊(如下圖)。 　　??　　現(xiàn)在，C君希望你告訴他隊伍整齊時能看到的學(xué)生人數(shù)。

Input

　　共一個數(shù)N。

Output

　　共一個數(shù)，即C君應(yīng)看到的學(xué)生人數(shù)。

Sample Input

　　4

Sample Output

　　9

HINT

?

【數(shù)據(jù)規(guī)模和約定】 　　對于 100% 的數(shù)據(jù)，1 ≤ N ≤ 40000

思路：求gcd（x,y)==1 ?1<=x,y<=N;

　　 ? ans=2*sigma(phi[i])+2-1;

　　　ps:我以為N==1的時候為0；然而并沒有1這樣的數(shù)據(jù)；

#include<bits/stdc++.h> using namespace std; #define ll long long #define esp 1e-10 const int N=1e5+10,M=1e6+10,mod=1e9+7,inf=1e9+10; ll p[M],ji; bool vis[M]; ll phi[M]; void get_eular(int n) {ji = 0;memset(vis, true, sizeof(vis));for(int i = 2; i <= n; i++){if(vis[i]){p[ji ++] = i;phi[i] = i - 1;}for(int j = 0; j < ji && i * p[j] <= n; j++){vis[i * p[j]] = false;if(i % p[j] == 0){phi[i * p[j]] = phi[i] * p[j];break;}elsephi[i * p[j]] = phi[i] * phi[p[j]];}} } ll sumphi[M]; int main() {ll x,y,z,i,t;get_eular(40100);phi[1]=1;sumphi[0]=0;for(i=1;i<=40000;i++)sumphi[i]=sumphi[i-1]+phi[i];while(~scanf("%lld",&x)){if(x==1)printf("0\n");elseprintf("%lld\n",2*sumphi[x-1]+1);}return 0; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/jhz033/p/5702839.html

總結(jié)

以上是生活随笔為你收集整理的bzoj 2190: [SDOI2008]仪仗队 线性欧拉函数的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。