博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
P1495 【模板】中国剩余定理(CRT)/曹冲养猪
阅读量:2153 次
发布时间:2019-04-30

本文共 1416 字,大约阅读时间需要 4 分钟。

题目描述

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

输入格式

第一行包含一个整数 n——建立猪圈的次数,接下来 n 行,每行两个整数 a i   b i a_i\ b_i ai bi
表示建立了 a i a_i ai 个猪圈,有 b i b_i bi 头猪没有去处。你可以假定 a i a_i ai a j a_j aj互质。

输出格式

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

输入输出样例

输入

33 15 17 2

输出

16

说明/提示

1 ≤ n ≤ 10 , 0 ≤ b i < a i ≤ 1000000 , 1 ≤ ∏ a i ≤ 1 0 18 1 \leq n\le10,0 \leq b_i\lt a_i\le1000000,1 \leq \prod a_i \leq 10^{18} 1n100bi<ai10000001ai1018

题解

不是吧不是吧不会有人看不出来 这是模板CRT,

即做法为套公式:
m = ∏ a i ( 即 最 小 公 倍 数 )          t = m / a [ i ]          x = t − 1 ( m o d    a i ) m=\prod{a_i}(即最小公倍数)\ \ \ \ \ \ \ \ t=m/a[i] \ \ \ \ \ \ \ \ x=t^{-1}(mod\ \ a_i) m=ai()        t=m/a[i]        x=t1(mod  ai)
a n s = ∑ x ∗ t ∗ b [ i ] ans=\sum{x*t*b[i]} ans=xtb[i]
其中逆元用扩展欧几里得求
当然这只是CRT的其中一个解,要求最小解还得模 m m m
具体的CRT证明这里不多讲,本人是看一篇懂的

代码

#include
#include
#define ll long long //不开longlong见祖宗using namespace std;ll n,m=1,ans,t,x,y,a[20],b[20];void exgcd(ll a,ll b)//扩展欧几里得算法求逆元{
if(b==0) {
x=1; y=0; return; } exgcd(b,a%b); int t=x; x=y; y=t-a/b*y;}int main(){
scanf("%lld",&n); for(int i=1; i<=n; i++) {
scanf("%lld%lld",&a[i],&b[i]); m*=a[i];//求最小公倍数 } for(int i=1; i<=n; i++) {
t=m/a[i]; x=y=0; exgcd(t,a[i]); x=(x%a[i]+a[i])%a[i]; ans=(ans+x*t*b[i])%m;//CRT 记得模数 }cout<

转载地址:http://xhpwb.baihongyu.com/

你可能感兴趣的文章
TP5.1项目从windows的Apache服务迁移到linux的Nginx服务需要注意几点。
查看>>
win10安装软件 打开时报错 找不到 msvcp120.dll
查看>>
PHPunit+Xdebug代码覆盖率以及遇到的问题汇总
查看>>
PHPUnit安装及使用
查看>>
PHP项目用xhprof性能分析(安装及应用实例)
查看>>
composer安装YII
查看>>
Sublime text3快捷键演示
查看>>
sublime text3 快捷键修改
查看>>
关于PHP几点建议
查看>>
硬盘的接口、协议
查看>>
VLAN与子网划分区别
查看>>
Cisco Packet Tracer教程
查看>>
02. 交换机的基本配置和管理
查看>>
03. 交换机的Telnet远程登陆配置
查看>>
微信小程序-调用-腾讯视频-解决方案
查看>>
phpStudy安装yaf扩展
查看>>
密码 加密 加盐 常用操作记录
查看>>
TP 分页后,调用指定页。
查看>>
Oracle数据库中的(+)连接
查看>>
java-oracle中几十个实用的PL/SQL
查看>>