其他分享
首页 > 其他分享> > P3868 [TJOI2009]猜数字(CRT板子)

P3868 [TJOI2009]猜数字(CRT板子)

作者:互联网

题目描述

现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示。其中第二组中的数字是两两互素的。求最小的非负整数n,满足对于任意的i,n - ai能被bi整除。

输入输出格式

输入格式:

 

输入数据的第一行是一个整数k,(1 ≤ k ≤ 10)。接下来有两行,第一行是:a1,a2,...,ak,第二行是b1,b2,...,bk

 

输出格式:

 

输出所求的整数n。

 

输入输出样例

输入样例#1: 复制
3
1 2 3
2 3 5
输出样例#1: 复制
23

说明

所有数据中,第一组数字的绝对值不超过109(可能为负数),第二组数字均为不超过6000的正整数,且第二组里所有数的乘积不超过1018

每个测试点时限1秒

注意:对于C/C++语言,对64位整型数应声明为long long,如使用scanf, printf函数(以及fscanf, fprintf等),应采用%lld标识符。


 


 


 


 

 

 

 

 

 1 #include"bits/stdc++.h"
 2 using namespace std;
 3 typedef long long ll;
 4 
 5 ll k,a[11],b[11],s=1,ans=0;
 6 
 7 
 8 void exgcd(ll a,ll b,ll &x, ll &y)
 9 {
10     if (!b){
11         x=1,y=0 ;return ;
12     }
13     
14     exgcd(b,a%b,x,y); ll t;
15     t=x;
16     x=y;
17     y=t-(a/b)*y;
18 }
19 
20 inline ll mul(ll a,ll b)
21 {
22     ll r=0;
23     while (b)
24     {
25         if (b&1)r+=a%=s;
26         b>>=1;
27         a+=a%=s;
28     }return r;
29 }
30 int main()
31 {
32     cin>>k;
33     for (int i=1;i<=k;i++)cin>>a[i];
34     for (int i=1;i<=k;i++)cin>>b[i],s*=b[i];
35     for(int i=1;i<=k;i++)
36     {
37         ll x,y;
38         exgcd(s/b[i],b[i],x,y);
39         x=(x%b[i]+b[i])%b[i];(1)
40     //     cout<<i<<" "<<x<<endl;
41         ans+=mul(s/b[i]*x,((a[i]%b[i]+b[i])%b[i]));//快速乘,否则爆long long(2)
// 一和二式都取模意义下的最小正整数
42 ans%=s; 43 } 44 cout<<ans; 45 }

 

 

 

 


 


 

标签:...,数字,CRT,int,ll,long,P3868,TJOI2009,第二组
来源: https://www.cnblogs.com/zhangbuang/p/10338369.html