codeforces-359-C：Prime Number

A_Thinking_Reed_

于 2019-09-12 15:59:45 发布

阅读量170

点赞数

分类专栏：模拟/贪心数学 | 思维

本文链接：https://blog.csdn.net/A_Thinking_Reed_/article/details/100774326

版权

数学 | 思维同时被 2 个专栏收录

12 篇文章 0 订阅

订阅专栏

模拟/贪心

10 篇文章 0 订阅

订阅专栏

题目：http://codeforces.com/problemset/problem/359/C

题意：给非增序列a[]和x，求 $\sum 1/x^{a_{i}}$ 的分子和分母的gcd。

思路： $\sum_{i=1}^{n} \frac{1}{x^{a_{i}}} = \frac{\sum_{i=1}^{n}x^{s-a_{i}}}{x^{s}}\cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot(s=\sum_{i=1}^{n}a_{i})$ ，显然 $gcd = x^\min_{i=1}^{n}(s-a_{i})$ ，不过当最小的个数超过x时，会对x+1造成一个贡献。当最小的个数能被x整除时，需要再向更大的考虑。注意分子有可能大于分母，例如样例8 2， 2 2 2 2 2 2 2 2，ans = 65536。

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn = 1e5+5;
const LL mod = 1e9+7;

LL quick_pow(LL a, LL b){
    LL sum = 1;
    while(b){
        if(b&1) sum = sum*a%mod;
        a = a*a%mod; b>>=1;
    }return sum;
}
LL n, x, sum, a[maxn];
int main()
{
    cin >> n >> x;
    for(int i=1; i<=n; i++) scanf("%lld", &a[i]), sum += a[i];
    for(int i=1; i<=n; i++) a[i] = sum-a[i];
    LL p, num = 1;
    for(int i=n-1; i>=0; i--){
        if(a[i] == a[i+1]) num ++;
        else {
            if(num % x == 0) num = num/x, a[i+1]++, i++;
            else {p = a[i+1]; break;}
        }
    }
    cout << quick_pow(x, min(p,sum)) << endl;
}