/ SeriousOJ /

Record Detail

Accepted


  
# Status Time Cost Memory Cost
#1 Accepted 66ms 16.625 MiB
#2 Accepted 75ms 16.469 MiB
#3 Accepted 71ms 16.105 MiB
#4 Accepted 84ms 16.191 MiB
#5 Accepted 68ms 16.863 MiB
#6 Accepted 71ms 16.844 MiB
#7 Accepted 81ms 16.203 MiB
#8 Accepted 68ms 16.488 MiB
#9 Accepted 83ms 16.809 MiB
#10 Accepted 86ms 16.602 MiB
#11 Accepted 69ms 16.871 MiB
#12 Accepted 65ms 16.406 MiB
#13 Accepted 77ms 16.488 MiB
#14 Accepted 76ms 16.828 MiB
#15 Accepted 135ms 16.773 MiB
#16 Accepted 143ms 17.695 MiB
#17 Accepted 351ms 17.23 MiB
#18 Accepted 172ms 16.734 MiB
#19 Accepted 243ms 16.672 MiB
#20 Accepted 271ms 17.289 MiB
#21 Accepted 241ms 17.227 MiB
#22 Accepted 242ms 17.164 MiB
#23 Accepted 258ms 16.672 MiB
#24 Accepted 314ms 17.07 MiB
#25 Accepted 200ms 17.031 MiB
#26 Accepted 262ms 16.516 MiB
#27 Accepted 543ms 16.512 MiB
#28 Accepted 217ms 17.445 MiB
#29 Accepted 303ms 16.547 MiB
#30 Accepted 685ms 17.098 MiB
#31 Accepted 590ms 17.098 MiB
#32 Accepted 179ms 16.484 MiB
#33 Accepted 139ms 17.613 MiB
#34 Accepted 159ms 17.363 MiB
#35 Accepted 993ms 17.098 MiB

Code

#pragma GCC optimize("O3,unroll-loops,Ofast")
// #pragma GCC target("avx2")
#include<bits/stdc++.h>
// #include <ext/pb_ds/assoc_container.hpp>
// #include <ext/pb_ds/tree_policy.hpp>
// using namespace __gnu_pbds;
#define int long long
#define endl '\n'

using namespace std;
using pii = pair<int, int>;
using tup = tuple<int, int, int>;

// mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
// template <class T> using ordered_set = tree<T, null_type,
//                          less<T>, rb_tree_tag, tree_order_statistics_node_update>;

const int inf = (1e15)+5;
const int mod = 1000000007;
const double eps = 1e-6;
const int N = 100001;

vector<int> divs[N];
void preprocess() {
	for (int i = 1; i < N; i++) {
		for (int j = i; j < N; j += i) {
			divs[j].push_back(i);
		}
	}
}

int n, a[N];
int cnt[N];


int calc(int x) {
	int tot = 0, out = 0;
	for(int i=1; i<=n; i++) {
		if(a[i] % x == 0) tot++;
		else out = gcd(out, a[i]);
	}

	int onno = n - tot;
	int lage;

	if(n % 2 == 0) {
		if(tot < n/2) return 0;
		lage = n / 2;
	}

	else {
		if(tot >= n/2+1) lage = n/2;
		else if(tot == n/2) lage = n/2+1;
		else return 0;
	}

	lage -= onno;

	for(int i=1; i<=n; i++) {
		if(a[i] % x) continue;
		for(int j : divs[a[i]])
			cnt[j]++;
	}

	// cout << tot << ' ' << onno << ' ' << lage << ' ' << out << endl;
	int ret = 0;
	for(int i=1; i<N; i++) {
		if(out % i or cnt[i] < lage) continue;
		ret = max(ret, i);
	}

	for(int i=1; i<=n; i++) {
		if(a[i] % x) continue;
		for(int j : divs[a[i]])
			cnt[j]--;
	}
	// if(out % x == 0) ret *= x;
	
	// cout << x << ": " << ret + x << endl;
	return ret + x;
}

void solve(int tc) {
	cin >> n;
	for(int i=1; i<=n; i++) cin >> a[i];

	bool f = 1;
	for(int i=2; i<=n; i++) f &= (a[i] == a[i-1]);
	if(f) {
		cout << 2 * a[1] << endl;
		return;
	}

	int who = a[1];
	for(int i=2; i<=n; i++) {
		if(divs[a[i]].size() < divs[who].size())
			who = a[i];
	}
	int ans = 0;
	for(int i : divs[who]) ans = max(ans, calc(i));

	cout << ans << endl;
}
    
int32_t main() {
    cin.tie(NULL)->sync_with_stdio(false);
    cout.precision(10);

    preprocess();

    int T = 1;
    cin >> T;

    for (int i = 1; i <= T; i++) {
        solve(i);
    }

    return 0;
}

Information

Submit By
Type
Submission
Problem
P1077 Even Odd GCD (Hard Version)
Contest
Bangladesh 2.0
Language
C++17 (G++ 13.2.0)
Submit At
2024-08-16 17:18:54
Judged At
2024-08-16 17:18:54
Judged By
Score
100
Total Time
993ms
Peak Memory
17.695 MiB