Codeforces —— 359C Prime Number

很有趣的一个数论题。

题目要求的是分子和分母的最大公约数，但因为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;
}