Codechef November Challenge 2015 Simple Sum

题意：

求 $\sum_{K=1}^N \frac{N}{gcd(N,K)}$. 有 $10^6$ 个询问，$N \le 10^7$ 。

思路：

首先证明出
答案应该是$\sum d\cdot\varphi(d)$ 对所有的 $d$ 为 $N$ 的约数。

这个具体的证明可以参考 SPOJ LCMSUM (数论)，然后如果我们先线性预处理出了$1$ - $10^7$ 内的欧拉函数 $phi[i]$，后面再枚举约数，考虑每个约数的贡献算出所有的答案，这样的话最短时间也是$3s$ 左右，但是题目要求$1s$.

注意到关于欧拉函数的下列性质：

所以我们可以利用性质1和2，做到$log N$内回答每一个询问：
$N = k * p ^e$， p为最大的质数，这个可以预处理的。
由性质1:

$$ans = \sum_{d|N} d\varphi(d)=\sum_{\alpha=0}^e \sum_{d|k} dp^\alpha \varphi(dp^\alpha)=\sum_{\alpha=0}^e \sum_{d|k} d\varphi(d) p^\alpha\varphi(p^\alpha) = \left(\sum_{\alpha=0}^e p^\alpha\varphi(p^\alpha)\right) \left(\sum_{d|k} d\varphi(d)\right)\,,$$

第一部分可以由性质2得到，第二部分就直接计算咯（也可以实现预处理$ans[i]$，$1\le i \le 10^6$）！

代码：

#include <bits/stdc++.h>
#define dibs reserve
#define OVER9000 1234567890
#define ALL_THE(CAKE,LIE) for(auto LIE =CAKE.begin(); LIE != CAKE.end(); LIE++)
#define tisic 47
#define soclose 1e-8
#define chocolate win
// so much chocolate
#define patkan 9
#define ff first
#define ss second
#define abs(x) ((x < 0)?-(x):x)
#define uint unsigned int
#define dbl long double
#define pi 3.14159265358979323846
using namespace std;
// mylittledoge

int main()
{
    cin.sync_with_stdio(0);
    cin.tie(0);
    int mx =10000000;
    vector<int> minD(mx+tisic,0);
    for(int i =2; i <= mx; i++) if(minD[i] == 0)
        for(int j =1; j <= mx/i; j++) minD[i*j] =i;
    int T;
    scanf(" %d",&T);
    while(T--) {
        int N;
        long long ans =1;
        scanf(" %d",&N);
        while(N > 1) {
            int p =minD[N], ppw =1;
            while(N%p == 0) {
                ppw *=p;
                N /=p;
            }
            ans *=1LL*ppw*ppw-(1LL*ppw*ppw-1)/(p+1);
        }
            printf("%lld\n",ans);
    }
    return 0;
}