bzoj 2705: [SDOI2012]Longge的问题—

Description

Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。

Input

一个整数，为N。

Output

一个整数，为所求的答案。

Sample Input

Sample Output

HINT

【数据范围】

对于60%的数据，0<N<=2^16。

对于100%的数据，0<N<=2^32。

———————————————————————

这道题如果一个数x gcd(n,x)==y 那么gcd(b/y,x/y)==1

所以我们枚举因数d 求一下1-n/d有多少个数和n/d的gcd为1 这个可以用欧拉函数

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
const int M=1e3+7;
LL read(){LL ans=0,f=1,c=getchar();while(c<'0'||c>'9'){if(c=='-') f=-1; c=getchar();}while(c>='0'&&c<='9'){ans=ans*10+(c-'0'); c=getchar();}return ans*f;
}
LL n,v,ans;
int p[M],cnt;
LL f(LL x){for(int i=1;i<=cnt;i++)if(x%p[i]==0) x=x/p[i]*(p[i]-1); return x;}
int main(){n=read(); v=n;for(LL x=2;x*x<=v;x++)if(v%x==0){p[++cnt]=x;while(v%x==0) v/=x;}if(v!=1) p[++cnt]=v;for(LL x=1;x*x<=n;x++)if(n%x==0){LL y=n/x;ans=ans+y*f(x);if(x!=y) ans=ans+x*f(y);}printf("%lld\n",ans);return 0;
}

View Code

转载于:https://www.cnblogs.com/lyzuikeai/p/7664991.html

bzoj 2705: [SDOI2012]Longge的问题——欧拉定理相关推荐

bzoj 2705: [SDOI2012]Longge的问题（欧拉函数）
2705: [SDOI2012]Longge的问题 Time Limit: 3 Sec Memory Limit: 128 MB Submit: 3077 Solved: 1914 [Submit ...
P2303 [SDOI2012] Longge 的问题
P2303 [SDOI2012] Longge 的问题思路我们显然可以枚举每一对数的gcdgcdgcd进行求解,进而我们有如下推导: =>∑i=1ngcd(i,n)=>\sum _{i ...
P2303 [SDOI2012] Longge(数论/欧拉函数）
P2303 [SDOI2012] Longge 一道看似非常基础的数论题,但是蕴含了非常多的知识,求解 ∑i=1ngcd(i,n)\sum_{i=1}^ngcd(i,n) i=1∑ngcd(i,n) ...
BZOJ 2707: [SDOI2012]走迷宫 [高斯消元 scc缩点]
2707: [SDOI2012]走迷宫题意:求s走到t期望步数,\(n \le 10^4\),保证\(|SCC| \le 100\) 求scc缩点,每个scc高斯消元,scc之间直接DP 注意每次清 ...
BZOJ.3004.[SDOI2012]吊灯(结论)
题目链接 BZOJ 洛谷题意: 将树划分为k个连通块,要求每个连通块大小相同.输出可能的大小. 结论: 满足条件时颜色的连通块数为k,当且仅当有 \(n/k\) 个节点满足它的子树是k的倍数(显然还 ...
BZOJ2705 [SDOI2012]Longge的问题欧拉函数
Time Limit: 3 Sec Memory Limit: 128 MB Submit: 3874 Solved: 2469 [Submit][Status][Discuss] Descrip ...
【bzoj2705】[SDOI2012]Longge的问题欧拉函数
题目描述 Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题.现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N). 输入一个整数,为N. 输出 ...
BZOJ 2707: [SDOI2012]走迷宫( tarjan + 高斯消元 )
数据范围太大不能直接高斯消元, tarjan缩点然后按拓扑逆序对每个强连通分量高斯消元就可以了. E(u) = 1 + Σ E(v) / degree(u) 对拍时发现网上2个程序的INF判断和我不一 ...
bzoj 2707: [SDOI2012]走迷宫（Trajan+高斯消元+Dp）
2707: [SDOI2012]走迷宫 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 862 Solved: 328 [Submit][Statu ...

bzoj 2705: [SDOI2012]Longge的问题——欧拉定理