牛客 计数器(裴蜀定理)

博客围绕一道题展开,需用裴蜀定理求解加若干次a1,a2...an后模m有多少不同值。将问题转换为a1*x1+a2*x2+...+an*xn的形式,依据裴蜀定理推广得出其值为gcd(a1,a2,...,an)的倍数,进一步推导得出最终结果与gcd(d,m)有关。

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

思路:这个题我们需要用到裴蜀定理,题目要求加若干次a1,a2...an然后模m有多少个不同的值,可以转换为a1*x1+a2*x2+...+an*xn有多少个不同的值,我将a1*x1+a2*x2+...+an*xn的值设为c,很显然,这个式子是裴蜀定理的推广,裴蜀定理(或贝祖定理,Bézout's identity)得名于法国数学家艾蒂安·裴蜀,说明了对任何整数a、b和它们的最大公约数d,关于未知数x和y的线性不定方程(称为裴蜀等式):若a,b是整数,且gcd(a,b)=d,那么对于任意的整数x,y,ax+by都一定是d的倍数,特别地,一定存在整数x,y,使ax+by=d成立,其推广则是对于任意整数x1,x2,....,xn,a1*x1+a2*x2+...+an*xn都一定是gcd(a1,a2,...,an)的倍数。对于这道题,我们设d=gcd(a1,a2,...,an),通过定理不难知道a1*x1+a2*x2+...+an*xn值为k*d(k是个不定值),而我们要求的是(k*d)%m有多少个不同的值,可以设(d*k)%m=p(p由d,k,m决定),然后可以将这个式子转化为d*k-m*t=p(p<m),从而变成了a*x+b*y的形式,从而利用上述定理推出p=gcd(d,m)*k'.。

代码:

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
const int maxn=1e2+9;
ll a[maxn];
ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    ll i,j,k,n,m;
    cin>>n>>m;
    for(i=0;i<n;i++){
        cin>>a[i];
        a[i]=((a[i]%m)==0?m:a[i]%m);
    }
    ll d=m;
    for(i=0;i<n;i++){
        d=gcd(d,a[i]);
    }
    ll ans=m/d;
    cout<<ans<<endl;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值