結果

問題 No.335 門松宝くじ
コンテスト
ユーザー toma
提出日時 2019-09-24 23:16:06
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 150 ms / 2,000 ms
コード長 2,302 bytes
コンパイル時間 1,950 ms
コンパイル使用メモリ 191,900 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-19 14:24:46
合計ジャッジ時間 3,313 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include"bits/stdc++.h"
using namespace std;
#define REP(k,m,n) for(int (k)=(m);(k)<(n);(k)++)
#define rep(i,n) REP((i),0,(n))
using ll = long long;
constexpr ll INF = 1ll << 60;

template<typename T>
class SegmentTree {
private:
	using F = function<T(T, T)>; // モノイド型
	int n; // 横幅
	F f;   // モノイド
	T e;   // モノイド単位元
	vector<T> data;

public:
	// init忘れに注意
	SegmentTree() {}
	SegmentTree(F f, T e) :f(f), e(e) {}
	void init(int n_) {
		n = 1;
		while (n < n_)n <<= 1;
		data.assign(n << 1, e);
	}
	void build(const vector<T>& v) {
		int n_ = v.size();
		init(n_);
		rep(i, n_)data[n + i] = v[i];
		for (int i = n - 1; i >= 0; i--) {
			data[i] = f(data[(i << 1) | 0], data[(i << 1) | 1]);
		}
	}
	void set_val(int idx, T val) {
		idx += n;
		data[idx] = val;
		while (idx >>= 1) {
			data[idx] = f(data[(idx << 1) | 0], data[(idx << 1) | 1]);
		}
	}
	T query(int a, int b) {
		// [a,b)
		T vl = e, vr = e;
		for (int l = a + n, r = b + n; l < r; l >>= 1, r >>= 1) {
			if (l & 1)vl = f(vl, data[l++]); // unknown
			if (r & 1)vr = f(data[--r], vr); // unknown
		}
		return f(vl, vr);
	}
};

ll estimate(vector<ll>& e) {
	// seg init
	const int N = e.size();
	auto f1 = [](ll a, ll b) {return max(a, b); };
	auto f2 = [](ll a, ll b) {return min(a, b); };
	SegmentTree<ll> seg_max(f1, -INF), seg_min(f2, INF);
	seg_max.build(e), seg_min.build(e);

	// map init
	map<ll, ll> mp;
	rep(i, N)mp[e[i]] = i;

	// query
	ll res = 0;
	REP(b, 1, N + 1)REP(a, 1, b) {
		ll l = mp[a], r = mp[b];
		ll tres = 0;
		if (l < r) {
			ll big = seg_max.query(0, r);
			if (big > b)tres = big;
			else {
				if (seg_max.query(0, l) > a ||
					seg_min.query(r, N) < b ||
					seg_min.query(l, N) < a)
					tres = b;
			}
		}
		else {
			swap(l, r);
			ll big = seg_max.query(l, N);
			if (big > b)tres = big;
			else {
				if (seg_max.query(r, N) > a ||
					seg_min.query(0, l) < b ||
					seg_min.query(0, r) < a)
					tres = b;
			}
		}
		res += tres;
	}
	return res;
}

int main()
{
	// input
	int N, M;
	cin >> N >> M;
	vector<vector<ll>> E(M, vector<ll>(N));
	rep(i, M)rep(j, N)cin >> E[i][j];

	vector<pair<ll, ll>> score;
	rep(i, M) {
		score.push_back({ estimate(E[i]),-i });
	}
	sort(score.begin(), score.end());
	cout << -score.back().second << endl;

	return 0;
}
0