很有趣的一个数论题。
题目要求的是分子和分母的最大公约数,但因为x是素数,所以其实就是要将分子和分母都表示成(x^k)*m,m不能被x整除,然后取分子和分母k的最小值。
对于分母,题目已经限定好就是x^s,s=a1+a2+...+an
所以关键点就是处理分子那一部分。
首先对于每个1/x^ai,通分后分子的形式就是x^(s-ai),我们的任务就是要求出分子中所有项的最大公约数(就是我们要的x^k),那肯定是取次数最小的那一项了。
但由于其中某些项相加有可能是系数能被x整除,
比如2^2+2^2=2*2^2=2^3,
对于这种情况分子求出来的次数是3而不是2
所以我们还要对分子的项进行处理。
处理的方法是将相同次数的项加到一块,这里我用map处理
然后每次取次数最小的t,判断系数r是否能被x整除,
如果不可以表示这个t就是我们要找的分子的t,
否则就将r中一个x乘到x^t上去,也就是将r/x加到x^(t+1)的系数上去,直到找到一个最小的t并且r%x!=0
这个过程就像那种三个空瓶可以换一瓶新饮料那样,喝完后新的空瓶加上原来的空瓶再去换。
这样我们就找到分子和分母的最大次数,取两者的最小值,再用快速幂就能求出答案了。
#include<cstdio>
#include<map>
using namespace std;
#define MOD 1000000007
#define LL __int64
int n, i;
LL a[100000], x, s, t, r;
map<LL,LL> MP;
map<LL,LL>::iterator it;
void solve(){
LL ans=1;
while(s>0){
if(s%2==1) ans = (ans*x)%MOD;
x = (x*x)%MOD;
s/=2;
}
printf("%I64d\n", ans);
}
int main(){
while(~scanf("%d %I64d", &n, &x)){
MP.clear();
s=0;
for(i=0; i<n; i++){
scanf("%I64d", &a[i]);
s+=a[i];
}
for(i=0; i<n; i++){
MP[s-a[i]]++;
}
while(1){
it = MP.begin();
t = it->first;
r = it->second;
MP.erase(t);
if(r%x>0){
if(t<s) s=t;
solve();
break;
}
r/=x;
t++;
MP[t]+=r;
}
}
return 0;
}