Submission #1258168

#TimeUsernameProblemLanguageResultExecution timeMemory
1258168allin27xTriple Peaks (IOI25_triples)C++20
70 / 100
139 ms23112 KiB
#include "triples.h"
#include <bits/stdc++.h>
using namespace std;
#define int long long

int ans = 0;
int n;
const int SQRT = 450;

ofstream fout("optfor10");

void count1(vector<signed> &h){
	for (int idx = 0; idx<n; idx++){
		int x = h[idx];
		int idh = idx - x;
		if (idh < 0) continue;
		int y = h[idh] - x;
		if (y <= 0 || y==x) continue;
		int idy = idx + y;
		if (idy >= n) continue;
		if (h[idy] == y) ans++;
	}
}

void count2(vector<signed> &h){
	for (int idx = 0; idx<n; idx++){
		int x = h[idx];
		int idy = idx + x;
		if (idy >= n) continue;
		int y = h[idy];
		int idh = idx - y;
		if (idh < 0) continue;
		if (h[idh] == x + y) ans++;
	}	
}

void count3(vector<signed> &h){
	for (int idx = 0; idx<n; idx++){
		int x = h[idx];
		int idh = idx + x;
		if (idh >= n) continue;
		int y = h[idh] - x;
		if (y <= 0 || y==x) continue;
		int idy = idh + y;
		if (idy >= n) continue;
		if (h[idy] == y) ans++;
	}
}

void count4(vector<signed> &h){
	vector<vector<int>> tl(3*n);
	for (int i=0; i<n; i++){
		int j = i + h[i]; if (j >= 3*n) continue;
		tl[j].push_back(h[i]);
	}
	
	for (int idt = 0; idt<3*n; idt++){
		if (tl[idt].size() > SQRT){
			for (int idx = 0; idx < idt; idx++){
				int x = h[idx];
				int y = idt - idx; y -= x; 
				if (y<=0 || (y&1)) continue;
				y/=2; int idy = idt - y;
				int idh = idt - x - y;
				if (h[idh] == x+y && h[idy] == y) ans++;
			}
		} else {
			for (int y: tl[idt]){
				for (int hv: tl[idt]){
					int x = hv - y; if (x<=0) continue;
					int idx = idt - hv - y;
					if (idx <0) continue;
					if (h[idx] == x) ans++;

				}
			}
		}
	}
}


long long count_triples(std::vector<signed> h) {
	n = h.size();
	auto rh = h; reverse(rh.begin(), rh.end());
	count1(h); count1(rh);
	count2(h); count2(rh);
	count3(h); 
	count4(h); 
	return ans;
}

vector<signed> tr;
vector<vector<signed>> best;
int nr = 0; int last = 0;
void brut(int idx){
	if (idx == 10){
		nr++;
		int nw = count_triples(tr);
		if (nw > ans) best.clear(), ans = nw;
		if (nw==ans && best.size() < 200000) best.push_back(tr);
		int pr = nr / 1000000;
		if (pr >last) last = pr, cout<<last<<endl;
		return;
	}
	for (int i=1; i<=9; i++) tr[idx] = i, brut(idx+1);


}
std::vector<signed> construct_range(signed n, signed K) {
	tr.resize(n);
	return {1,1,1};
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...