日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

【SDOI2008】仪仗队

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

題面

題解

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

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

或者用歐拉函數(shù)

代碼

#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

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎勵來咯,堅持創(chuàng)作打卡瓜分現(xiàn)金大獎

總結

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

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