• Prime Number CodeForces


    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 ≤ 1052 ≤ 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).

    Example
    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, as351 = 13·27729 = 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.

    #include<queue>
    #include<stack>
    #include<vector>
    #include<math.h>
    #include<stdio.h>
    #include<numeric>//STL数值算法头文件
    #include<stdlib.h>
    #include<string.h>
    #include<iostream>
    #include<algorithm>
    #include<functional>//模板类头文件
    using namespace std;
    
    const long long INF=1e9+7;
    const long long maxn=101000;
    
    long long n,x;
    long long a[maxn];
    
    long long quick_mod(long long a,long long b)
    {
        long long ans=1;
        a=a%INF;
        while(b)
        {
            if(b&1)
                ans=ans*a%INF;
            a=a*a%INF;
            b>>=1;
        }
        return ans;
    }
    
    int main()
    {
        while(~scanf("%d %d",&n,&x))
        {
            long long sum1=0;
            for(long long i=0; i<n; i++)
            {
                scanf("%d",&a[i]);
                sum1+=a[i];
            }
            for(long long i=0; i<n; i++)
                a[i]=sum1-a[i];
            sort(a,a+n);
            long long ans,j=1,cot=1,t;
            for(j=1; j<=n; j++)
            {
                if(a[j]!=a[j-1])
                {
                    if(cot%x)
                    {
                        ans=a[j-1];
                        break;
                    }
                    else
                    {
                        cot/=x;
                        a[j-1]+=1;
                        j--;
                    }
                }
                else cot++;
            }
            printf("%d
    ",quick_mod(x,min(ans,sum1)));
        }
        return 0;
    }
    
    
    
    
    #include<queue>
    #include<stack>
    #include<vector>
    #include<math.h>
    #include<stdio.h>
    #include<numeric>//STL数值算法头文件
    #include<stdlib.h>
    #include<string.h>
    #include<iostream>
    #include<algorithm>
    #include<functional>//模板类头文件
    using namespace std;
    
    const long long INF=1e9+7;
    const long long maxn=101000;
    
    long long n,x;
    long long a[maxn];
    
    long long gcd(long long a,long long b)
    {
        long long ans=1;
        ans=ans%INF;
        while(b)
        {
            if(b&1)
                ans=ans*a%INF;
            a=a*a%INF;
            b>>=1;
        }
        return ans;
    }
    
    int main()
    {
        while(~scanf("%I64d %I64d",&n,&x))
        {
            long long sum1=0;
            for(long long i=0; i<n; i++)
            {
                scanf("%I64d",&a[i]);
                sum1+=a[i];
            }
            for(long long i=0; i<n; i++)
                a[i]=sum1-a[i];
            sort(a,a+n);
            long long ans,j=1,cot=1,t;
            while(j<=n)
            {
                if(a[j]!=a[j-1])
                {
                    if(cot%x)
                    {
                        ans=a[j-1];
                        break;
                    }
                    long long f=a[j-1]+1;
                    t=cot/x;
                    for(long long k=j-1,s=t; s>0; s--,k--)
                        a[k]=f;
                    j-=t;
                    j++;
                    cot=1;
                }
                else cot++,j++;
            }
            printf("%I64d
    ",gcd(x,min(ans,sum1)));
        }
        return 0;
    }
    













  • 相关阅读:
    (转)Unity3D 开发优秀技术资源汇总
    (转)Unity3d通过Action注册事件,回调方法
    (转)Unity3D研究院之游戏架构脚本该如何来写(三十九)
    (转)Unity3D研究院之异步加载游戏场景与异步加载游戏资源进度条(三十一)
    Unity3D的主要类图
    C# 事件和Unity3D
    unity3d 场景配置文件生成代码
    497. Random Point in Non-overlapping Rectangles
    478. Generate Random Point in a Circle
    470. Implement Rand10() Using Rand7() (拒绝采样Reject Sampling)
  • 原文地址:https://www.cnblogs.com/nyist-xsk/p/7264860.html
Copyright © 2020-2023  润新知