2020.02.19普及C组模拟赛8（第二题）

2.家庭做业（zadaca）

题目描述

Mirko最近收到了一个家庭做业，做业的任务是计算两个数A和B的最大公约数。因为这两个数太大了，咱们给出了N个数，他们的乘积是A，给出M个数，他们的乘积是B。

Mirko想要验算本身的答案，因此他想找你写一个程序来解决这个问题。若是这个最大公约数超过了9位数，那么只须要输出最后9位就能够了。

输入

第一行包含一个正整数N，范围是1到1000。第二行是N个用空格隔开的正整数(小于10亿），他们的乘积是A。第三行包含一个正整数M，范围是1到1000。第四行是M个用空格隔开的正整数（小于10亿），他们的乘积是B。

输出

输出有且只有一行，表示A和B的最大公约数，若是结果超过了9位数，输出最后9位数就能够了。

样例输入和输出

正解
用b去和每一个a[i]找最大公因数，将这些最大公因数相乘，b和a[i]都要除以这个最大公因数。
输出就用
printf("%09lld",s);
AC代码

#include<iostream>
#include<cstdio>
using namespace std;
int n,m,f;
long long b,k,s,a[1005];
long long gys(long long x,long long y)//找最大公因数（展转相除法）
{
	long long z;
	while(x%y!=0)
	{
		z=x%y;
		x=y;
		y=z;
	}
	return y;
}
int main()
{
	freopen("zadaca.in","r",stdin);
	freopen("zadaca.out","w",stdout);
	s=1;
	cin>>n;
	for(int i=1;i<=n;i++)
	 cin>>a[i];
	cin>>m;
	for(int i=1;i<=m;i++) 
	{
		cin>>b;
		for(int j=1;j<=n;j++)
		{
			k=gys(a[j],b);
			a[j]/=k;//除以最大公因数
			b/=k;
			if(s*k>=1000000000)f=1;//特判
			s=s*k%1000000000;
		}
	}
	if(f==1)printf("%09lld",s);
	else cout<<s;
	return 0;
}

下面附本次比赛的其余题目

2020.02.19普及C组模拟赛8（第一题）
2020.02.19普及C组模拟赛8（第二题）
2020.02.19普及C组模拟赛8（第三题）
2020.02.19普及C组模拟赛8（第四题）
2020.02.19普及C组模拟赛8（总结）

谢谢