結果

問題 No.335 門松宝くじ
コンテスト
ユーザー airis
提出日時 2016-01-16 12:55:44
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 323 ms / 2,000 ms
コード長 4,051 bytes
コンパイル時間 989 ms
コンパイル使用メモリ 94,552 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-19 19:56:05
合計ジャッジ時間 3,656 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <numeric>
#include <bitset>
#include <complex>
#define rep(x, to) for (int x = 0; x < (to); x++)
#define REP(x, a, to) for (int x = (a); x < (to); x++)
#define foreach(itr, x) for (typeof((x).begin()) itr = (x).begin(); itr != (x).end(); itr++)
#define EPS (1e-14)

using namespace std;

typedef long long ll;
typedef pair<int, int> PII;
typedef pair<ll, ll> PLL;
typedef complex<double> Complex;
typedef vector< vector<int> > Mat;

int N, M;
int E[10][805];
int ans_mx, ans_id;

bool is_kadomatsu(int a, int b, int c) {
	if (a < b && b > c) return true;
	if (a > b && b < c) return true;
	else return false;
}

struct SegmentTreeMin {
		int n;
		int dat[100005];
		void init(int m) {
			n = 1;
			while (n < m) {
				n *= 2;
			}
			for (int i = 0; i < 2 * n - 1; i++) {
				dat[i] = (int)(1e+9 + 7);
			}
		}
		void update(int i, int value) {
			i += n - 1;
			dat[i] = value;
			while (i > 0) {
				i = (i - 1) / 2;
				int chl = 2 * i + 1;
				int chr = 2 * i + 2;
				dat[i] = min(dat[chl], dat[chr]);
			}
		}
		//[a, b) node-k [l, r)
		int query(int a, int b, int k, int l, int r) {
			if (a <= l && r <= b) return dat[k];
			if (b <= l || r <= a) return (int)(1e+9 + 7);
			int chl = 2 * k + 1;
			int chr = 2 * k + 2;
			return min(query(a, b, chl, l, (l + r)/2), query(a, b, chr, (l + r)/2, r));
		}
};

struct SegmentTreeMax {
		int n;
		int dat[100005];
		string type;
		void init(int m) {
			n = 1;
			while (n < m) {
				n *= 2;
			}
			for (int i = 0; i < 2 * n - 1; i++) {
				dat[i] = 0;
			}
		}
		void update(int i, int value) {
			i += n - 1;
			dat[i] = value;
			while (i > 0) {
				i = (i - 1) / 2;
				int chl = 2 * i + 1;
				int chr = 2 * i + 2;
				dat[i] = max(dat[chl], dat[chr]);
			}
		}
		//[a, b) node-k [l, r)
		int query(int a, int b, int k, int l, int r) {
			if (a <= l && r <= b) return dat[k];
			if (b <= l || r <= a) return 0;
			int chl = 2 * k + 1;
			int chr = 2 * k + 2;
			return max(query(a, b, chl, l, (l + r)/2), query(a, b, chr, (l + r)/2, r));
		}
};

SegmentTreeMin seg_min;
SegmentTreeMax seg_max;

ll calc(int p) {
	int res = 0;
	seg_min.init(N);
	seg_max.init(N);
	for (int i = 0; i < N; i++) {
		seg_min.update(i, E[p][i]);
		seg_max.update(i, E[p][i]);
	}
#if 0
	for (int i = 0; i < N; i++) {
		for (int j = i + 1; j < N; j++) {
			printf("[%d,%d) min %d\n", i, j+1, seg_min.query(i, j+1, 0, 0, seg_min.n));
			printf("[%d,%d) max %d\n", i, j+1, seg_max.query(i, j+1, 0, 0, seg_max.n));
		}
	}
#endif
	for (int i = 0; i < N; i++) {
		for (int j = i + 1; j < N; j++) {
			int mx = 0;
			int x = 0;
			//[i, j)
			x = seg_min.query(i, j, 0, 0, seg_min.n);
			if (is_kadomatsu(E[p][i], x, E[p][j])) {
				mx = max(mx, max(x, max(E[p][i], E[p][j])));
			}
			x = seg_max.query(i, j, 0, 0, seg_max.n);
			if (is_kadomatsu(E[p][i], x, E[p][j])) {
				mx = max(mx, max(x, max(E[p][i], E[p][j])));
			}
			if (E[p][i] < E[p][j]) {
				//[0, i)
				x = seg_max.query(0, i, 0, 0, seg_max.n);
				if (is_kadomatsu(x, E[p][i], E[p][j])) {
					mx = max(mx, max(x, E[p][j]));
				}
				//[j, N)
				x = seg_min.query(j, N, 0, 0, seg_min.n);
				if (is_kadomatsu(E[p][i], E[p][j], x)) {
					mx = max(mx, E[p][j]);
				}
			} else {
				// E[p][i] > E[p][j]
				// [0, i)
				x = seg_min.query(0, i, 0, 0, seg_min.n);
				if (is_kadomatsu(x, E[p][i], E[p][j])) {
					mx = max(mx, E[p][i]);
				}
				// [j, N)
				x = seg_max.query(j, N, 0, 0, seg_max.n);
				if (is_kadomatsu(E[p][i], E[p][j], x)) {
					mx = max(mx, max(E[p][i], x));
				}
			}
#if 0
			printf("%d,%d => %d\n", i, j, mx);
#endif
			res += mx;
		}
	}
	return res;
}

void solve() {
	for (int i = 0; i < M; i++) {
		int mx = calc(i);
		if (ans_mx < mx) {
			ans_mx = mx;
			ans_id = i;
		}
	}
	cout<< ans_id << endl;
}

int main() {
	cin >> N >> M;
	rep(i, M) rep(j, N) {
		cin >> E[i][j];
	}
	solve();
	return 0;
}
0