CodeForces - 964C（等比数列求和+逆元）

最新推荐文章于 2020-09-01 16:38:44 发布

原创最新推荐文章于 2020-09-01 16:38:44 发布 · 460 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#ACM #CodeForces

数论专栏收录该内容

3 篇文章

订阅专栏

本文提供了一道CodeForces平台上的题目964C的解答思路及代码实现，该题要求计算特定表达式的模10^9+9的值，通过等比数列求和结合逆元的方法解决。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

链接：CodeForces - 964C

题意：求 $n \sum i = 0 s_{i} a^{n - i} b^{i} by 10^{9} + 9 值。s[i]为+1或-1。$

$题解：可证：每k项的和相差(b/a)^k。等比数列求和加逆元。$

$逆元不一定要最后一步用的，而是每步都可以用。$

#include <bits/stdc++.h>
using namespace std;

const double EPS = 1e-8;
const int mod = 1e9 + 9;
const int INF = 0x3f3f3f3f;
const int maxn = 1e5 + 10;
long long n, a, b, k;
char s[maxn];

long long pow_mod(long long x, long long n)
{
    long long res = 1;
    while(n){
        if(n & 1) res = res * x % mod;
        x = x * x % mod;
        n >>= 1;
    }
    return res;
}

void Exgcd(long long a, long long b, long long& x, long long& y)
{
    if(!b){
        y = 0; x = 1;
        return ;
    }
    Exgcd(b, a % b, y, x);
    y -= a / b * x;
}

long long NY(long long a, long long b)
{
    long long x, y;
    Exgcd(a, b, x, y);
    return (x % mod + mod) % mod;
}

int main()
{
    scanf("%lld%lld%lld%lld", &n, &a, &b, &k);
    scanf("%s", s);

    long long ans = 0;
    for(int i = 0; i < k; i++){
        long long tmp = pow_mod(a, n - i) * pow_mod(b, i) % mod;
        if(s[i] == '+') ans += tmp;
        else ans -= tmp;
        ans = (ans % mod + mod) % mod;
    }

    long long q = pow_mod(b, k) * NY(pow_mod(a, k), mod) % mod;
    long long m = (n + 1) / k;

    if(q == 1) ans = ans * m % mod;
    else ans = ans * (pow_mod(q, m) - 1) % mod * NY(q - 1, mod) % mod;

    ans = (ans % mod + mod) % mod;

    printf("%lld\n", ans);;

    return 0;
}