求幂（京东2018秋招）

东东对幂运算很感兴趣,在学习的过程中东东发现了一些有趣的性质: 9^3 = 27^2, 2^10 = 32^2
东东对这个性质充满了好奇,东东现在给出一个整数n,希望你能帮助他求出满足 a^b = c^d(1 ≤ a,b,c,d ≤ n)的式子有多少个。
例如当n = 2: 1^1=1^1
1^1=1^2
1^2=1^1
1^2=1^2
2^1=2^1
2^2=2^2
一共有6个满足要求的式子

输入描述:

输入包括一个整数n(1 ≤ n ≤ 10^6)

输出描述:

输出一个整数,表示满足要求的式子个数。因为答案可能很大,输出对1000000007求模的结果

示例1

输入

2

输出

6

思路：暴力肯定会超时，有什么技巧呢？a^b=c^d可以写成这样的形式(i^x)^b=(i^y)^d，暴力i，同时遍历x和y（只要满足i^x和i^y不大于n即可），遍历x和y 的时候把最大公约数拿出来放到底数上，那么剩余的部分分别和b和d相乘，剩余部分最大的*（b或d）不大于n，所有可以求出一个最大取值范围。

/**
a^c=b^d----->(i^x)^c = (i^y)^d
遍历i，且i^x<=n、j^y<=n,求出x和y的范围
遍历x和y，取出x和y的最大公约数放到底数做i的指数，
x和y的最大值除以最大公约数m和一个参数(c或d)相乘的到不大于n的数，即n/m
*/
#include<iostream>
#include<set>
using namespace std;
int MOD = 1000000000 + 7;
int gcd(int a,int b){
    return (a % b == 0) ? b : gcd(b,a%b);
}
int main() {
    int n;
    cin>>n;
    //int会越界，且还需要强制类型转换
    long long ans =(long long) 1*n*(n*2-1) % MOD;
    set<int> s;
    //遍历i
    for (int i = 2; i*i <= n; i++){
        if ( s.find(i)!=s.end()) continue;
        int tmp = i;
        int cnt = 0;
        //满足条件的底数，同时求x和y的范围
        while(tmp <= n) {
            s.insert(tmp);
            tmp = tmp * i;
            cnt++;
        }
        //相同的公约数，放到下面，反正不会大于n
        for(int k = 1; k <= cnt; k++) {
            for(int j = k + 1; j <= cnt; j++) {
                ans = (ans + n / (j / gcd(k, j) ) * (long)2 ) % MOD;
            }
        }
    }
    cout<<ans<<endl;
    return 0;
}

/**
#特殊数据：不满足幂关系的（底数为1，底数不为1全相等）
①底数为1的单独处理，指数为1的无所谓
1^ =1^ (n个数中选择2个)(前后可以互换位置)
C(2,n)=n*(n-1)/2
A(2,n)=n*(n-1)
②底数不为1，只能是相等（(n-1)*n）
*/
/**
①1^1=1^1、1^1=1^2、1^2=1^1、1^2=1^2
②2^1=2^1、2^2=2^2
*/