[JLOI2013]地形生成

JLOI2013过了好长时间，才写第四题。。

第一问比较好想。

第二问我想到了n^3次方的做法，但是数据。。。。于是没敢写，然后上网查了一下题解，居然是O(n^3)过的，数据这么弱。。。

/* 

 * Problem: JLOI2013-Terrain

 * Author: Shun Yao

 */



#include <string.h>

#include <stdlib.h>

#include <limits.h>

#include <assert.h>

#include <stdio.h>

#include <ctype.h>

#include <math.h>

#include <time.h>



#include <map>

#include <set>

#include <list>

#include <stack>

#include <queue>

#include <deque>

#include <string>

#include <vector>

#include <bitset>

#include <utility>

#include <iomanip>

#include <numeric>

#include <sstream>

#include <iostream>

#include <algorithm>

#include <functional>



//using namespace std;



const int MAXN = 1010;

const int MOD = 2011;



int n;

std::map<int, int> ocr;



struct Data {

	int h, k;

} a[MAXN];



bool cmp(Data a, Data b) {

	return a.h != b.h ? a.h > b.h : a.k < b.k;

}



int main(/*int argc, char **argv*/) {

	int i, j, t, f[2][MAXN], ans1, ans2, cur, k;

	

	freopen("terrain.in", "r", stdin);

	freopen("terrain.out", "w", stdout);

	

	scanf("%d", &n);

	for (i = 1; i <= n; ++i) {

		scanf("%d%d", &a[i].h, &a[i].k);

		++ocr[a[i].h];

	}

	std::sort(a + 1, a + n + 1, cmp);

	ans1 = ans2 = 1;

	for (i = 1; i <= n; ) {

		t = 0;

		for (j = 0; j < ocr[a[i].h]; ++j)

			ans1 = ans1 * (std::min(i, a[i + j].k) + t++) % MOD;

		memset(f[0], 0, sizeof f[0]);

		f[0][1] = 1;

		cur = 0;

		for (j = 0; j < ocr[a[i].h]; ++j) {

			memset(f[1 - cur], 0, sizeof f[1 - cur]);

			for (k = 1; k <= std::min(i, a[i + j].k); ++k)

				f[1 - cur][k] = (f[1 - cur][k - 1] + f[cur][k]) % MOD;

			cur = 1 - cur;

		}

		t = 0;

		for (j = 1; j <= i; ++j)

			t = (t + f[cur][j]) % MOD;

		ans2 = ans2 * t % MOD;

		i += ocr[a[i].h];

	}

	printf("%d %d", ans1, ans2);

	

	fclose(stdin);

	fclose(stdout);

	return 0;

}