欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

Codeforces Round #511 (Div. 2), problem: (C) Enlarge GCD 数论

程序员文章站 2022-03-02 23:05:55
...

思路大体是先找到an取值范围内的质数并保存,求出输入的数据中最大公倍数并将a数组都除以最大公倍数,然后遍历an的取值范围(从2开始,1不算),如果是质数就统计这个质数和它的倍数,就是使最小公倍数增加所需要去掉的数的个数,遍历过程中取最小值即可。并注意判定特殊情况。这个方法几乎要超时了似乎有更好的办法,但是我没想到。

ac代码:

#include <bits/stdc++.h>
#define FOR(I,A,B) for(int I = (A); I < (B); I++)
#define FORE(I,A,B) for(int I = (A); I <= (B); I++)
#define PRII pair<int,int> 
#define LL long long 

#define INF 1000000001
#define N 200005
#define M 15000000
int gcd(int a,int b) { return b?gcd(b,a%b):a;}
using namespace std;
int a[300005];
int k[M],ss[M];
int main()
{
	int n;
	int ans = INF;
	cin >> n;
	cin >> a[1] >> a[2];
	bool fy = true;
	int agcd = gcd(a[1],a[2]);
	if(a[1] != a[2]) fy = false;
	FORE(i,3,n){
		cin >> a[i];
		agcd = gcd(a[i],agcd);
		if(a[i] != a[i-1]) fy = false;
	}
	FORE(i,2,int(sqrt(M)+1)){
		if(ss[i] == 0){
			for(int j = i*2;j <= M;j+=i) ss[j] = 1;
		}
	}
	FORE(i,1,n){
		k[a[i]/agcd]++;
	} 
相关标签: codeforces