POJ1573(X问题中国剩余定律)

最新推荐文章于 2019-08-03 23:37:18 发布

原创最新推荐文章于 2019-08-03 23:37:18 发布 · 437 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#多元线性同余方程 #中国剩余定律

数论专栏收录该内容

19 篇文章

订阅专栏

本文概述了AI音视频处理领域的关键技术，包括视频分割、语义识别、自动驾驶、AR、SLAM等，并探讨了其在实际应用中的作用。

X问题

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4234 Accepted Submission(s): 1351

Problem Description

求在小于等于N的正整数中有多少个X满足：X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 < a[i] <= 10)。

Input

输入数据的第一行为一个正整数T，表示有T组测试数据。每组测试数据的第一行为两个正整数N，M (0 < N <= 1000,000,000 , 0 < M <= 10)，表示X小于等于N，数组a和b中各有M个元素。接下来两行，每行各有M个正整数，分别为a和b中的元素。

Output

对应每一组输入，在独立一行中输出一个正整数，表示满足条件的X的个数。

Sample Input

3
10 3
1 2 3
0 1 2
100 7
3 4 5 6 7 8 9
1 2 3 4 5 6 7
10000 10
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 7 8 9

Sample Output

1
0
3

AC代码:

#include<iostream>
#include<cstdio>
using namespace std;

typedef long long ll;

const int maxn = 11;

int a[maxn],b[maxn];

void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if(!b)
	{
		x = 1;
		y = 0;
		d = a;
	}
	else
	{
		ex_gcd(b,a%b,d,y,x);
		y -= x * (a / b);
	}
}

long long gcd(long long a,long long b)
{
	if(b == 0)
		return a;
	return gcd(b,a % b);
}

int main()
{
	ll i,a1,r1,a2,r2,ans,first,second,third,temp,d,x0,y0;
	int T;
	int m,n;
	//freopen("1.txt","r",stdin);
	scanf("%d",&T);
	while(T--)
	{
		bool ifhave = 1;
		scanf("%d%d",&m,&n);
		for(i=0;i<n;i++)
		{
			scanf("%d",&a[i]);
		}
		for(i=0;i<n;i++)
		{
			scanf("%d",&b[i]);
		}
		first = a[0];
		r1 = b[0];
		for(i=1;i<n;i++)
		{
			second = a[i];
			third = b[i] - r1;
			ex_gcd(first,second,d,x0,y0);
			if(third % d != 0)
			{
				ifhave = 0;
				break;
			}
			int t = second / d;
			x0 = (x0 * (third / d) % t + t) % t;
			r1 = first * x0 + r1;
			first = first * (second / d);
		} //r1为最小的解
		long long temp1,temp2,ans;
		int lcm = 1;
		for(i=0;i<n;i++)
		{
			lcm = lcm / gcd(lcm,a[i]) * a[i];
			//cout<<lcm<<endl;
		}
		//cout<<"lcm:"<<lcm<<endl;
		if(!ifhave)
		{
			printf("0\n");
			continue;
		}
		long long res = 0;
		if(r1 <= m)
		{
			res = 1 + (m - r1) / lcm;
		}
		if(res && r1 == 0)
		{
			res--;
		}
		printf("%lld\n",res);
	}
	return 0;	
}