Atcoder Beginer Contest 162 E - Sum of gcd of Tuples (Hard)

Atcoder Beginer Contest 162 E - Sum of gcd of Tuples (Hard)

一.题目大意
有$N$个数，每一个数都在$[1,K]$之间，共有$K^N$种不同的排列，对于每一种排列求$gcd$，把每个$gcd$加起来，模$10^9+7$。
二.大体思路
不能暴力，否则~~~
考虑对于$x$的情况：
$$gcd(A_1,A_2,...,A_N)=x$$
满足($k_i表示x是A_i的k_i倍$)：
$$A_1=x{k}_1,A_2=x{k}_2,...,A_N=x{k}_N$$
每一个$k_i$共有$[K/x]$种取值，一共有$N$个，所以有$[K/x{]}^N$中排列方式。
不过，这样会挂~~
因为不能全是$x$的倍数的倍数，所以要减掉$x=x的倍数的倍数的情况$（懵 ）
即若$y$是$x$的倍数（不是$x$），则要减掉$y$的情况数（可以用$an{s}_y$存储）
三.上代码！！！

#include "bits/stdc++.h"
using namespace std;
const int maxN=1e5+10;
const long long MOD=1e9+7;
int n;
long long k;
long long ans[maxN];
long long res=0LL;
long long fast_mi(long long a,int p)
{
	if(p==0)
	{
		return 1LL;
	}
	long long tmp=fast_mi(a,p>>1);
	tmp=tmp*tmp%MOD;
	if(p&1)
	{
		tmp=tmp*a%MOD;
	}
	return tmp;
}
int main()
{
	cin>>n>>k;
	for(long long x=k;x>=1LL;x--)
	{
		long long multi_num=k/x;
		long long num=fast_mi(multi_num,n);
		for(int i=x+x;i<=k;i+=x)
		{
			num=(num-ans[i]+MOD)%MOD;
		}
		ans[x]=num;
		res=(res+num*x)%MOD;
	}
	cout<<res<<endl;
	return 0;
}

EOF