Poj2480欧拉函数

枚举n的约数d，∑d*phi(d) 就是所求答案，剩下的就是参考别人的证明。

化简 p^i*phi(p^(k-i)) 可得 p^k - p^(k-1) ，注意特判 k==i的情况，注意LL。

#define _CRT_SECURE_NO_WARNINGS
#pragma comment(linker, "/STACK:102400000,102400000")
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <string>
#include <vector>
#include <cmath>
#include <queue>
#include <map>
#include <set>
using namespace std;typedef long long LL;LL gao( LL sum, LL k, LL p)
{LL ans = 0;LL t = sum;ans += k*(t - t/p);ans += t;return ans;
}int main()
{LL n;
//    n = 1<<31;
//    cout<<n<<endl;while(scanf("%I64d",&n)!=EOF){LL t = n; LL ans = 1;for( LL i = 2;i*i<=t;i++){if(t%i) continue;LL cnt =0;  LL sum = 1;while(t%i==0){t/=i;cnt++;sum*=i;}ans *= gao(sum,cnt,i);}if(t>1) ans*=gao(t,1,t);printf("%I64d\n",ans);}return 0;
}

转载于:https://www.cnblogs.com/yigexigua/p/4756909.html

Poj2480欧拉函数相关推荐

poj2480（利用欧拉函数的积性求解）
题目链接: http://poj.org/problem?id=2480 题意:∑gcd(i, N) 1<=i <=N,就这个公式,给你一个n,让你求sum=gcd(1,n)+gcd(2 ...
POJ-2480 Longge's problem 欧拉函数
详见代码: #include <cstdlib> #include <cstring> #include <cstdio> #include <cmath&g ...
poj2154-color-polyan次二面体+欧拉函数优化
N<=1e9,O(nlogn)的做法会超时.从枚举置换转变为枚举轮换长度,然后可以利用欧拉函数,把复杂度变为O(√n * logn) 1 /*-------------------------- ...
hdu 1286( 欧拉函数）
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1286 数学题真的是有点吃不消了... View Code 1 #include<iostream ...
费马定理中值定理_数论-欧拉函数、欧拉定理
欧拉函数积性函数满足 ( 互质) 定义对于正整数 ,欧拉函数是小于等于的所有数中与互质的数的个数. 欧拉函数是积性函数(这个证明不是很显然,这个链接里面有很多种证明方法) 记作: 公式 , ...
poj2154(Polya+欧拉函数优化模版）
#include <cstdio> #include <cstring> #include<iostream> using namespace std; const ...
【数学专题】约数个数与欧拉函数
整理的算法模板合集: ACM模板目录一.约数个数 1. AcWing 1291. 轻拍牛头 2. AcWing 1294. 樱花 2.1 AcWing 197. 阶乘分解 3. AcWing 19 ...
【数学知识】三种方法求 [1,n] 中所有数欧拉函数（线性筛欧拉函数优化至 O(n) ）
整理的算法模板合集: ACM模板 ①直接求小于或等于n,且与n互质的数个数(求[1,n]中所有数的欧拉函数时间复杂度:O(nn)O(n\sqrt{n})O(nn)) ②求[1,n]之间每个数的质因数 ...
bzoj 1409 Password 矩阵快速幂+欧拉函数
可以发现,该数组的mi就是斐波那契数列所以要矩阵快速幂搞出第n位但是斐波那契数列上涨的很快,这就需要欧拉定理了 p^phi(q)%q=1(gcd(p,q)==1) p是素数,所以可以用然后需要5 ...

Poj2480欧拉函数

Poj2480欧拉函数相关推荐

最新文章

热门文章