SPOJ-NUMTRYE 质因子分解 + 欧拉函数应用

最新推荐文章于 2017-10-18 12:42:39 发布

数论只会GCD

最新推荐文章于 2017-10-18 12:42:39 发布

阅读量635

点赞数

分类专栏：数论-欧拉函数 ---------ACM-数学------------- SPOJ 文章标签： ACM 数论欧拉函数质因子

本文链接：https://blog.csdn.net/ADjky/article/details/69544877

版权

---------ACM-数学------------- 同时被 3 个专栏收录

16 篇文章 0 订阅

订阅专栏

数论-欧拉函数

10 篇文章 1 订阅

订阅专栏

SPOJ

3 篇文章 0 订阅

订阅专栏

传送门：SPOJ -NUMTRYE

题解:

首先 $\sum_{i=1}^ngcd(n, i) = \sum_{d|n}d\phi(\frac{n}{d})$
所以 $g(n) = \sum_{i=1}^nn/gcd(n, i) = \sum_{d|n}\frac{n}{d}\phi(\frac{n}{d}) = \sum_{d|n}d\phi(d)$
根据欧拉函数定义和欧拉函数性质的证明过程：http://blog.csdn.net/adjky/article/details/61198817
枚举n的因子利用欧拉函数积性可得:
设 $n = p_1^{e_1}p_2^{e_2}···p_k^{e_k}$
那么
$设d|n, d = p_1^{a_1}p_2^{a_2}···p_k^{a_k} 所以\phi（d）= \prod_{i = 1}^kp_i^{a_i}- p_i^{a_i - 1}$
$那么d*\phi(d) = \prod_{i = 1}^kp_i^{2a_i}- p_i^{2a_i - 1}$
因此:
$g (n) = \sum d | n d ϕ (d) = \prod i = 1 k \sum j = 0 e i ϕ (p j i) = \prod i = 1 k (1 + \sum j = 1 e i p 2 j i - p 2 j - 1 i) = \prod i = 1 k p 2 e i + 1 i + 1 p i + 1$ $g(n) = \sum_{d|n}d\phi(d) = \prod_{i = 1}^k\sum_{j = 0}^{e_i}\phi(p_i^j) = \prod_{i = 1}^k(1 + \sum_{j = 1}^{e_i}p_i^{2j}- p_i^{2j - 1}) \\ = \prod_{i = 1}^k\frac{p_i^{2e_i+1} + 1} {p_i + 1}$
所以 $f(n)/g(n) = \prod_{i=1}^k(p_i+1)$

分解质因子即可:
注意最多有10000个cas所以要先处理1e6以内的质因子优化分解:

code:

#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
const ll p = 1e9 + 7;
ll n;
bool vis[1000010];
ll prime[1000010];
int cnt = 0;
void getPrime(){
    memset(vis, 0, sizeof vis);
    for(int i = 2; i <= 1000000; ++i){
        if(!vis[i]) prime[++cnt] = i;
        for(int j = 1; j <= cnt && prime[j] * i <= 1000000; ++j){
            vis[i * prime[j]] = 1;
            if(i % prime[j] == 0) break;
        }
    }
}
ll solve(){
    ll res = 1;
    for(int i = 1; i <= cnt && prime[i] * prime[i] <= n; ++i){
        if(n % prime[i] == 0){
            while(n % prime[i] == 0){
                n /= prime[i];
            }
            res = res * (prime[i] + 1) % p;
        }
    }
    if(n > 1) res = res * (n + 1) % p;
    return res;
}

int main(){
    getPrime();
    int t;
    cin >> t;
    while(t--){
        cin >> n;
        cout << solve() << endl;
    }

    return 0;
}