Chinese remainder theorem again(最小公倍数)

Chinese remainder theorem again

 我知道部分同学最近在看中国剩余定理,就这个定理本身,还是比较简单的:
假设m1,m2,…,mk两两互素,则下面同余方程组:
x≡a1(mod m1)
x≡a2(mod m2)
…
x≡ak(mod mk)
在0<=<m1m2…mk内有唯一解。
记Mi=M/mi(1<=i<=k),因为(Mi,mi)=1,故有二个整数pi,qi满足Mipi+miqi=1,如果记ei=Mi/pi,那么会有:
ei≡0(mod mj),j!=i
ei≡1(mod mj),j=i
很显然,e1a1+e2a2+…+ekak就是方程组的一个解,这个解加减M的整数倍后就可以得到最小非负整数解。
这就是中国剩余定理及其求解过程。
现在有一个问题是这样的:
一个正整数N除以M1余(M1 - a),除以M2余(M2-a), 除以M3余(M3-a),总之, 除以MI余(MI-a),其中(a<Mi<100 i=1,2,…I),求满足条件的最小的数。

Input
输入数据包含多组测试实例,每个实例的第一行是两个整数I(1 < I < 10)和a,其中,I表示M的个数,a的含义如上所述,紧接着的一行是I个整数M1,M1…MI,I=0 并且a=0结束输入,不处理。
Output
对于每个测试实例,请在一行内输出满足条件的最小的数。每个实例的输出占一行。
Sample Input

2 1
2 3
0 0

Sample Output

5

题意:

让我们求一个数N,使得 N%mi=mia N % m i = m i − a

分析:

这是一道披着中国剩余定理外衣的最小公倍数的题

我们分析式子

N%mi=mia N % m i = m i − a

可以写成同余式

Nmia (mod mi) N ≡ m i − a   ( m o d   m i )

N+ami (mod mi) N + a ≡ m i   ( m o d   m i )

即:

N+a0 (mod mi) N + a ≡ 0   ( m o d   m i )

所以N+a是 mi m i 的倍数,所以N+a应该是 m1,m2mi m 1 , m 2 ⋯ m i 的公倍数,因为题目要求最小整数
所以就是求 m1,m2mi m 1 , m 2 ⋯ m i 的最小公倍数

注意这里在求的过程中要边除去gcd边相乘,否则会爆long long

code:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b){
    return b ? gcd(b,a%b) : a;
}
int main(){
    ll n,a;
    ll m[20];
    while(~scanf("%lld%lld",&n,&a)){
        if(n == 0 && a == 0) break;
        ll ans = 1;
        for(int i = 0; i < n; i++){
            ll num;
            scanf("%lld",&num);
            ans = ans / gcd(ans,num) * num;
        }
        printf("%lld\n",ans-a);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值