CodeForces - 359C Prime Number（数论）

Prime Number

Simon has a prime number x and an array of non-negative integers a1, a2, ..., an.

Simon loves fractions very much. Today he wrote out number  on a piece of paper. After Simon led all fractions to a common denominator and summed them up, he got a fraction: , where number t equals xa1 + a2 + ... + an. Now Simon wants to reduce the resulting fraction.

Help him, find the greatest common divisor of numbers s and t. As GCD can be rather large, print it as a remainder after dividing it by number 1000000007 (109 + 7).

Input

The first line contains two positive integers n and x (1 ≤ n ≤ 105, 2 ≤ x ≤ 109) — the size of the array and the prime number.

The second line contains n space-separated integers a1, a2, ..., an (0 ≤ a1 ≤ a2 ≤ ... ≤ an ≤ 109).

Output

Print a single number — the answer to the problem modulo 1000000007 (109 + 7).

Examples

input

2 2
2 2

output

8

input

3 3
1 2 3

output

27

input

2 2
29 29

output

73741817

input

4 5
0 0 0 0

output

1

Note

In the first sample . Thus, the answer to the problem is 8.

In the second sample, . The answer to the problem is 27, as 351 = 13·27, 729 = 27·27.

In the third sample the answer to the problem is 1073741824 mod 1000000007 = 73741817.

In the fourth sample . Thus, the answer to the problem is 1.

参考博客： http://blog.csdn.net/hongrock/article/details/14163805

题意：有一个素数x，另外有一个序列a1,a2,.....an，让x^a[i]相加求和得到一个分数，分母为x^a[1]*x^a[2]*.....x^a[n]，最后求分子和分母的最大公约数

很有想法的一道题，看大神的题解费了老半天劲才懂 ，好像离散数学前几章节就讲了类似的，（不过好像我并没有听0.0）

思路：

首先我们想啊，x是素数，所以最后分子和分母一定可以写成m*(x^k)的形式（m不能被x整除），所以最大公约数就是x^k，k为他们两个中比较小的那一个k

分母的k很容易找到，就是a[1]+a[2]+....+a[n]，主要是分子的k

我们令sum=a[1]+a[2]+....+a[n]，那么分子和分母可以写成这样的形式：（x^(sum-a[1])+x^(sum-a[2])+....+x^(sum-a[n])）/(x^sum)

很明显分子的k就是x^(sum-a[i])中最小的那一个sum-a[i]，但是由于分子中可能会出现两个相同的次数k

比如   2^2+2^2=2*2^2=2^3，

那么这种情况下的k就为3而不是2了

所以我们需要判断一下，假如将分子的次数k从小到大排序后，出现两个相同的k，那么就将这两项加到一块，让其系数加一

如果前后两项并不相同，那么我们需要先判断一下它（次数k小的那一项）的系数m能否被素数x整除，如果能整除说明次数k可以增大（让其+1），并让m/=x

如果不能整除，那么这一项的k就是我们需要的分子的k了（k已经从小到大排过序，所以第一次找到的肯定就是答案）

然后和分母的次数k比较一下取较小的，最后快速幂就ok了

代码：

#include<stdio.h>
#include<algorithm>
using namespace std;

#define min(a,b) (a<b?a:b)
#define LL long long int
#define maxn 100000+10
#define g 1000000007
LL a[maxn];

LL quick_mod(LL a,LL b)
{
    LL ans=1;
    a=a%g;
    while(b)
    {
        if(b&1)
            ans=ans*a%g;
        a=a*a%g;
        b>>=1;
    }
    return ans;
}

int main()
{
    LL x,n,i,sum=0;
    scanf("%lld%lld",&n,&x);
    for(i=0;i<n;i++)
    {
        scanf("%lld",&a[i]);
        sum+=a[i];//分母的次数k
    }
    for(i=0;i<n;i++)
        a[i]=sum-a[i];//分子的次数k
    sort(a,a+n);
    LL ans,cnt=1;//cnt为系数
    a[n]=-1;
    for(i=1;i<=n;i++)
    {
        if(a[i]!=a[i-1])
        {
            if(cnt%x)
            {
                ans=a[i-1];
                break;
            }
            else
            {
                cnt/=x;
                a[i-1]+=1;
                i--;
            }
        }
        else cnt++;
    }
    ans=min(ans,sum);
    printf("%lld\n",quick_mod(x,ans));
    return 0;
}