| # | Time | Username | Problem | Language | Result | Execution time | Memory |
|---|---|---|---|---|---|---|---|
| 37287 | DoanPhuDuc | Sailing Race (CEOI12_race) | C++98 | 1719 ms | 11460 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <bits/stdc++.h>
#define FOR(x, a, b) for (int x = a; x <= b; ++x)
#define FOD(x, a, b) for (int x = a; x >= b; --x)
#define REP(x, a, b) for (int x = a; x < b; ++x)
#define DEBUG(X) { cout << #X << " = " << X << endl; }
#define PR(A, n) { cout << #A << " = "; FOR(_, 1, n) cout << A[_] << " "; cout << endl; }
#define PR0(A, n) { cout << #A << " = "; REP(_, 0, n) cout << A[_] << " "; cout << endl; }
#define BitCount(x) __builtin_popcount(x)
using namespace std;
typedef long long LL;
typedef pair <int, int> II;
const int N = 5e2 + 10;
const int INF = 0x3f3f3f3f;
int n, request;
int CL1[N][N], CL2[N][N];
int dp[N][N][2], f[N][N][2], g[N][N][2];
bool e[N][N];
bool fre[N][N][2];
vector <int> adj[N], rev[N];
int O;
int CCW(int i, int j, int k) {
if (i == k || j == k) return -1;
if (i < j) {
return (j < k && k < n) || (k < i);
} else return (j < k && k < i);
}
bool CMP1(int x, int y) {
return CCW(O, x, y);
}
bool CMP2(int x, int y) {
return !CCW(O, x, y);
}
int DP(int l, int r, int k) {
if (l == r) return dp[l][r][k] = 0;
int &ans = dp[l][r][k];
if (ans != -1) return ans;
ans = 0;
REP(i, 0, adj[r].size()) {
int v = adj[r][i];
int d = CCW(l, r, v); if (d == -1) continue;
if (d == k) {
ans = max(ans, 1 + DP(l, v, k));
ans = max(ans, 1 + DP(r, v, k ^ 1));
}
}
return ans;
}
int Compute(int l, int r, int k) {
if (l == r) return 0;
int &ans = f[l][r][k];
if (ans != -1) return ans;
ans = e[l][r] ? 1 : -INF;
REP(i, 0, adj[l].size()) {
int u = adj[l][i];
if (CCW(l, r, u) == (k ^ 1))
ans = max(ans, 1 + Compute(u, r, k));
}
return ans;
}
void Init() {
memset(CL1, -1, sizeof CL1);
memset(CL2, -1, sizeof CL2);
REP(i, 0, n) {
O = i;
sort(rev[i].begin(), rev[i].end(), CMP1);
int p = 0;
for (int j = (i + 1) % n; j != i; j = (j + 1) % n) {
while (p < rev[i].size() && CCW(i, j, rev[i][p]) != 1) ++p;
if (p < rev[i].size()) CL1[i][j] = rev[i][p];
}
sort(rev[i].begin(), rev[i].end(), CMP2);
p = 0;
for (int j = (i - 1 + n) % n; j != i; j = (j - 1 + n) % n) {
while (p < rev[i].size() && CCW(i, j, rev[i][p]) != 0) ++p;
if (p < rev[i].size()) CL2[i][j] = rev[i][p];
}
}
}
int main() {
#ifdef LOCAL
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
#endif // LOCAL
scanf("%d%d", &n, &request);
REP(u, 0, n) {
int v;
while (scanf("%d", &v) == 1) {
if (v == 0) break;
e[u][v - 1] = true;
adj[u].push_back(v - 1);
rev[v - 1].push_back(u);
}
}
Init();
memset(dp, -1, sizeof dp);
int ans = 0, num = -1;
REP(u, 0, n) {
REP(k, 0, adj[u].size()) {
int v = adj[u][k];
int can = 1 + DP(u, v, 0);
g[u][v][0] = can;
if (can > ans) {
ans = can;
num = u;
}
can = 1 + DP(u, v, 1);
g[u][v][1] = can;
if (can > ans) {
ans = can;
num = u;
}
}
}
REP(s, 0, n) {
int p = s;
for (int t = (s - 1 + n) % n; t != s; t = (t - 1 + n) % n) {
g[s][t][1] = max(g[s][t][1], g[s][p][1]);
p = t;
}
p = s;
for (int t = (s + 1) % n; t != s; t = (t + 1) % n) {
g[s][t][0] = max(g[s][t][0], g[s][p][0]);
p = t;
}
}
if (request == 0) {
printf("%d\n%d", ans, num + 1);
return 0;
}
memset(f, -1, sizeof f);
REP(i, 0, n) {
for (int j = (i + 1) % n; j != i; j = (j + 1) % n) {
int x = Compute(i, j, 1); if (x == 0) continue;
int k = CL1[i][j]; if (k == -1) continue;
REP(p, 0, adj[j].size()) {
int v = adj[j][p];
if (CCW(j, k, v) == 1 && CCW(v, i, j) == 1) {
/* int can = 2 + x + max(g[v][(k + 1) % n][1],
g[v][(i - 1 + n) % n][0]);
if (can > ans) {
ans = can;
num = k;
}*/
int can = 2 + x + max(g[v][k][1], g[v][i][0]);
if (can > ans) {
ans = can;
num = k;
}
}
}
}
for (int j = (i - 1 + n) % n; j != i; j = (j - 1 + n) % n) {
int x = Compute(i, j, 0); if (x == 0) continue;
int k = CL2[i][j]; if (k == -1) continue;
REP(p, 0, adj[j].size()) {
int v = adj[j][p];
if (CCW(j, k, v) == 0 && CCW(v, i, j) == 0) {
/* int can = 2 + x + max(g[v][(i + 1) % n][1],
g[v][(k - 1 + n) % n][0]);
if (can > ans) {
ans = can;
num = k;
}*/
int can = 2 + x + max(g[v][i][1], g[v][k][0]);
if (can > ans) {
ans = can;
num = k;
}
}
}
}
}
printf("%d\n%d", ans, num + 1);
return 0;
}
Compilation message (stderr)
| # | Verdict | Execution time | Memory | Grader output |
|---|---|---|---|---|
| Fetching results... | ||||