国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

bzoj 2226: [Spoj 5971] LCMSum線性篩歐拉函數

2019-11-06 06:35:08
字體:
來源:轉載
供稿:網友

題目大意:給定n,求1到n中所有數與n的lcm之和 題解:枚舉d=GCD(i,n),令F(n)為n以內與n互質的數之和,則ans=Σ(d|n)d*F(d)*n/d=nΣF(d) F(d)有一個性質,就是與d互質的數一定能兩兩組合成d,可以用輾轉相除法輕松證明,只有1和2特殊,特判即可。

#include<cstdio>#include<cstdlib>#include<iostream>#include<iomanip>#include<ctime>#include<cmath>#include<cstring>#include<string>#include<algorithm>using namespace std;int n;long long f[1010000];bool pd[1010000];int my_stack[1010000];int top=0;void shai(){ f[1]=1; for(int i=2;i<=1000000;i++) { if(!pd[i]) { my_stack[++top]=i; f[i]=i-1; } for(int j=1;i*my_stack[j]<=1000000;j++) { pd[i*my_stack[j]]=true; if(i%my_stack[j]==0) { f[i*my_stack[j]]=f[i]*my_stack[j]; break; } f[i*my_stack[j]]=f[i]*(my_stack[j]-1); } } //for(int i=2;i<=1000000;i++) f[i]=i*f[i]/2;}inline long long F(int x){ if(x==1) return 1; return f[x]*x/2;}int main(){ //freopen("a.in","r",stdin); //freopen("a.out","w",stdout); int T; scanf("%d",&T); shai(); while(T--) { scanf("%d",&n); long long ans=0; int i; for(i=1;i*i<n;i++) { if(n%i==0) { ans+=F(i); ans+=F(n/i); } } if(i*i==n) ans+=f[i]*i/2;
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 札达县| 常德市| 麻栗坡县| 达孜县| 北宁市| 峨山| 陵水| 枝江市| 隆尧县| 陆良县| 西平县| 台北县| 右玉县| 五河县| 彩票| 石河子市| 常德市| 景谷| 镇巴县| 丽江市| 岳阳市| 淮南市| 边坝县| 白玉县| 龙江县| 富源县| 肃北| 新民市| 兴隆县| 临邑县| 邛崃市| 丰都县| 商水县| 凤冈县| 台江县| 汶上县| 密云县| 茶陵县| 华亭县| 清徐县| 天全县|