[LGOJ1495]曹冲养猪——[CRT]

【题目描述】

自从曹冲搞定了大象以后,曹操就开始琢磨让儿子干些事业,于是派他到中原养猪场养猪,可是曹冲很不高兴,于是在工作中马马虎虎,有一次曹操想知道母猪的数量,于是曹冲想狠狠耍曹操一把。
举个例子,假如有16头母猪,如果建了 3个猪圈,剩下1 头猪就没有地方安家了;如果建造了5 个猪圈,但是仍然有 1头猪没有地方去;如果建造了 7个猪圈,还有2 头没有地方去。你作为曹总的私人秘书理所当然要将准确的猪数报给曹总,你该怎么办?

【输入格式】

第一行包含一个整数n,表示建立猪圈的次数;

接下来n 行,每行两个整数 ,表示建立了 ai 个猪圈,有bi 头猪没有去处。你可以假定 ai,aj 互质。

【输出格式】

输出仅包含一个正整数,即为曹冲至少养猪的数目。

S a m p l e    I n p u t Sample~~Input Sample  Input

3
3 1
5 1
7 2

S a m p l e    O u t p u t Sample~~Output Sample  Output

16

【题意分析】

因为 ∀ i , j ∈ [ 1 , n ] , a i \forall i,j \in[1,n],a_i i,j[1,n],ai a j a_j aj互质,那么我们可以用CRT(中国剩余定理)解决此题

转化题意,CRT是求得同余方程组

{ x ≡ b 1 ( m o d   a 1 ) x ≡ b 2 ( m o d   a 2 ) . . . x ≡ b n ( m o d   a n ) \begin{cases} x \equiv b_1 (mod~a_1) \\ x \equiv b_2 (mod~a_2) \\ ...\\ x \equiv b_n (mod~a_n)\end{cases} xb1(mod a1)xb2(mod a2)...xbn(mod an)

的最小正整数解

先预处理出这n个模数的LCM,即

L C M = ∏ i = 1 n a i LCM= \prod_{i=1}^n a_i LCM=i=1nai

考虑把问题分割,对于第i个同余方程,令
u = L C M a i u=\frac{LCM}{a_i} u=aiLCM

先求得一个x,它是u的倍数,而且模ai等于bi
举个栗子,同余方程组

{ x ≡ 1   ( m o d   3 ) x ≡ 1   ( m o d   5 ) x ≡ 2   ( m o d   7 ) \begin{cases} x \equiv 1~(mod~3) \\ x \equiv 1~ (mod~5) \\ x \equiv 2~ (mod~7)\end{cases} x1 (mod 3)x1 (mod 5)x2 (mod 7)

那么LCM是3x5x7=105
先看第三个方程,u是105/7=15,那我们只要找到一个最小的15的倍数,它模7余2,是30,即方程

15 x ≡ 2   ( m o d   7 ) 15x\equiv2~(mod~7) 15x2 (mod 7)

的正最小整数解,等价于

15 x + 7 y = 2 15x+7y=2 15x+7y=2

显然珂以利用扩欧求解方程

15 x + 7 y = 1 15x+7y=1 15x+7y=1

的最小正整数解,那么将答案再乘上右边的系数2即可求得30。

如果我们把由第三个方程搞出的这个30记做 λ 3 \lambda_3 λ3,那么这个问题的最后答案就是

a n s = ( ∑ i = 1 n λ i ) % L C M ans=(\sum_{i=1}^n \lambda_i)\%LCM ans=(i=1nλi)%LCM

如果为负,加上若干倍LCM再取模即可。

Code:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#define MAXN 100
#define int long long
using namespace std;

int a[MAXN], b[MAXN], alpha[MAXN], x, y, n, LCM = 1, ans;

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

void CRT () {
	for (register int i = 1; i <= n; i++, x = 0, y = 0) {
		EXGCD (x, y, LCM / a[i], a[i]);
		ans += (x * (LCM / a[i]) * b[i]), ans = (ans % LCM + LCM) % LCM;
	}
	printf ("%lld\n", ans);
}

signed main () {
	scanf ("%lld", &n);
	for (register int i = 1; i <= n; i++) {
		scanf ("%lld%lld", &a[i], &b[i]); LCM *= a[i];
	}
	CRT ();
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值