日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

【SDOI2008】仪仗队

發布時間:2023/12/9 53 豆豆
生活随笔 收集整理的這篇文章主要介紹了 【SDOI2008】仪仗队 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

題面

題解

當$(x,y)$能被看到時,$gcd(x,y)=1$,

所以可以求$\sum_{i=0}^n\sum_{j=0}^n[gcd(x,y)=1]$

或者用歐拉函數

代碼

#include<bits/stdc++.h> #define RG register #define clear(x, y) memset(x, y, sizeof(x)); using namespace std;template<typename T> inline T read() {T data=0, w=1;char ch=getchar();while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();if(ch=='-') w=-1, ch=getchar();while(ch>='0'&&ch<='9') data=(data<<3)+(data<<1)+(ch^48), ch=getchar();return data*w; }const int maxn(40010); int phi[maxn], prime[maxn], cnt; bool is_prime[maxn];int getphi(int n) {for(RG int i=2;i<=n;i++){if(!is_prime[i]){prime[++cnt]=i;phi[i]=i-1;}for(RG int j=1;j<=cnt;j++){if(prime[j]*i>n) break;is_prime[prime[j]*i]=true;if(!(i%prime[j])) {phi[i*prime[j]]=phi[i]*prime[j];break;}else phi[i*prime[j]]=phi[i]*(prime[j]-1);}} }int n, ans; int main() {n=read<int>();getphi(n);if(n==1) return printf("0\n")&0;for(RG int i=3;i<=n;i++) ans+=phi[i-1];printf("%d\n", (ans<<1)+3);return 0; }

轉載于:https://www.cnblogs.com/cj-xxz/p/10185834.html

創作挑戰賽新人創作獎勵來咯,堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的【SDOI2008】仪仗队的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。