제출 #470978

#	제출 시각	아이디	문제	언어	결과	실행 시간	메모리
470978		Namnamseo	Naan (JOI19_naan)	C++17	29 / 100	183 ms	9064 KiB

naan

#include <algorithm>
#include <iostream>
#include <numeric>
using namespace std;
using ll=long long;
const int maxn = int(2e3) + 10;

int n, l;
ll v[maxn][maxn];
ll vs[maxn][maxn];

ll gcd(ll a, ll b) { while (b) a %= b, swap(a, b); return a; }

struct frac {
	ll m, j;
	frac(ll j_=0, ll m_=1) : m(m_), j(j_) {
		if (m < 0) m=-m, j=-j;
		ll g = gcd(m, abs(j));
		m /= g; j /= g;
	}
	explicit operator int() { return j/m; }
	explicit operator ll() { return j/m; }
	bool operator<(frac r) const { return j*r.m < m*r.j; }
	bool operator>(frac r) const { return j*r.m > m*r.j; }
	bool operator<=(frac r) const { return j*r.m <= m*r.j; }
	bool operator>=(frac r) const { return j*r.m >= m*r.j; }
	frac operator+(frac r) { return frac(j*r.m + m*r.j, m*r.m); }
	frac operator-(frac r) { return frac(j*r.m - m*r.j, m*r.m); }
	frac operator*(frac r) { return frac(j*r.j, m*r.m); }
	frac operator/(frac r) { return frac(j*r.m, m*r.j); }
	frac ceil_to(ll nm) { return frac((j*nm+m-1)/m, nm); }
};

int ansp[maxn];
bool done[maxn];

frac byn(int i, frac s) {
	frac tgt(vs[i][l], n);
	int sj = int(s)+1;
	frac lwidth = frac(sj) - s;
	frac lweight = frac(v[i][sj]) * lwidth;
	if (lweight >= tgt) {
		return (s + tgt / frac(v[i][sj])).ceil_to(n * v[i][sj]);
	}
	tgt = tgt - lweight; s = frac(sj);
	ll isang = vs[i][sj] + ll(tgt);
	int nj = upper_bound(vs[i]+1, vs[i]+l+1, isang) - vs[i] - 1;
	tgt = tgt - (vs[i][nj] - vs[i][sj]); s = nj;
	return s + (tgt / frac(v[i][nj+1]));
}

int main() { cin.tie(0)->sync_with_stdio(0);
	cin >> n >> l;
	for (int i=1; i<=n; ++i) for (int j=1; j<=l; ++j) cin >> v[i][j];
	for (int i=1; i<=n; ++i) partial_sum(v[i]+1, v[i]+l+1, vs[i]+1);

	frac cx{};
	for (int pi=1; pi<=n; ++pi) {
		frac mx; int mxi = -1;
		for (int i=1; i<=n; ++i) if (!done[i]) {
			frac tx = byn(i, cx);
			if (mxi == -1 || tx < mx) mx = tx, mxi = i;
		}
		cx = mx;
		if (pi < n) cout << mx.j << ' ' << mx.m << '\n';
		ansp[pi] = mxi;
		done[mxi] = true;
	}

	for (int i=1; i<=n; ++i) cout << ansp[i] << " \n"[i == n];

	return 0;
}

• Subtask 1: 5 / 5
• Subtask 2: 24 / 24
• Subtask 3: 0 / 71