Submission #489904

#	Time	Username	Problem	Language	Result	Execution time	Memory
489904		fhvirus	Sirni (COCI17_sirni)	C++17	140 / 140	1704 ms	534524 KiB

sirni

// Knapsack DP is harder than FFT.
#include<bits/stdc++.h>
using namespace std;
typedef long long ll; typedef pair<int,int> pii; typedef pair<ll,ll> pll;
#define ff first
#define ss second
#define pb emplace_back
#define AI(x) begin(x),end(x)
template<class I>bool chmax(I&a,I b){return a<b?(a=b,true):false;}
template<class I>bool chmin(I&a,I b){return b<a?(a=b,true):false;} 
#ifdef OWO
#define debug(args...) SDF(#args, args)
#define OIU(args...) ostream& operator<<(ostream&O,args)
#define LKJ(S,B,E,F) template<class...T>OIU(S<T...>s){O<<B;int c=0;for(auto i:s)O<<(c++?", ":"")<<F;return O<<E;}
LKJ(vector,'[',']',i)LKJ(deque,'[',']',i)LKJ(set,'{','}',i)LKJ(multiset,'{','}',i)LKJ(unordered_set,'{','}',i)LKJ(map,'{','}',i.ff<<':'<<i.ss)LKJ(unordered_map,'{','}',i.ff<<':'<<i.ss)
template<class...T>void SDF(const char* s,T...a){int c=sizeof...(T);if(!c){cerr<<"\033[1;32mvoid\033[0m\n";return;}(cerr<<"\033[1;32m("<<s<<") = (",...,(cerr<<a<<(--c?", ":")\033[0m\n")));}
template<class T,size_t N>OIU(array<T,N>a){return O<<vector<T>(AI(a));}template<class...T>OIU(pair<T...>p){return O<<'('<<p.ff<<','<<p.ss<<')';}template<class...T>OIU(tuple<T...>t){return O<<'(',apply([&O](T...s){int c=0;(...,(O<<(c++?", ":"")<<s));},t),O<<')';}
#else
#pragma GCC optimize("Ofast")
#define debug(...) ((void)0)
#endif
 
struct DSU {
	int n; vector<int> f;
	DSU (int nn) : n(nn), f(nn) { iota(AI(f), 0); }
	int F(int u) { return u == f[u] ? u : f[u] = F(f[u]); }
};

const int kC = 1e7;
vector<pii> bkt[kC];
 
signed main(){
	ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
 
	int N; cin >> N;
	vector<int> P(N); for(int &i: P) cin >> i;
 
	sort(AI(P));
	P.erase(unique(AI(P)), end(P));
	N = P.size();
 
	if(N <= 1000){
		for(int i = 0; i < N; ++i)
			for(int j = i+1; j < N; ++j)
				bkt[min(P[j] % P[i], P[i])].pb(i, j);
	} else {
		for(int i = 0; i < N-1; ++i){
			for(int j = P[i], lst = i+1; j <= P.back(); j += P[i]){
				int p = lower_bound(begin(P)+lst, end(P), j) - begin(P);
				int r = P[p] / P[i];
				bkt[P[p] - P[i] * r].pb(i, p);
				j = P[i] * r;
				lst = p+1;
			}
		}
	}
 
	DSU dsu(N);
	int ans = 0;
	int cnt = N;
	for(int w = 0; w < P.back(); ++w){
		for(auto &[u, v]: bkt[w]){
			u = dsu.F(u);
			v = dsu.F(v);
			if(u != v){
				dsu.f[v] = u;
				ans += w;
				--cnt;
			}
			if(cnt == 1) goto jizz;
		}
	}

jizz:
	cout << ans << '\n';
 
	return 0;
}

• Subtask 1: 14 / 14
• Subtask 2: 14 / 14
• Subtask 3: 14 / 14
• Subtask 4: 14 / 14
• Subtask 5: 14 / 14
• Subtask 6: 14 / 14
• Subtask 7: 14 / 14
• Subtask 8: 14 / 14
• Subtask 9: 14 / 14
• Subtask 10: 14 / 14