PACMAN吃毛豆

前几天做这个题，由于写的太差了，导致我写了9个版本，直到最后一个才过掉……

 

数一下我的几个版本：

ver1：考试时写出的第一版，用的是一个动态规划思想，但是调不对，放弃……

ver2：考试时写的第二版，用网络流建模，建模方法：若从i能到j，就从i向j连一条费用为-1，容量为1的边。每个点向汇连一条容量为1、费用为0的边，源向每个点连一条费用为-1、容量为1的边，超级源向源连一条费用为0、容量为2的边，求超级源到汇的费用流。WA掉。

ver3：看了解题报告后写的一个版本，改了建模方法，拆点，若从i能到j，就从i+n向j连一条费用为0，容量为1的边。i向i+n连一条费用为-1、容量为1的边，其余相同。

ver4：改了费用流算法，由原始对偶转为isap增广。

ver5：又改了费用流算法，转为push-relabel增广。

ver6：还改费用流算法，直接spfa增广。

ver7：骗分策略，先求一条最长上升子序列，再在这条线左右分别求一条，WA掉。

ver8：回归最大流，直接手动增广，先求一条最长上升子序列，模拟费用流增广，直接构图，不含0边，然后用spfa求最短路，用的ver2的构图方式，WA掉。

ver9：改变构图，用ver3的构图，终于AC……

 

附上源代码：

 

 

 

 

#include <stdio.h>#include <memory.h>#include <stdlib.h>#define maxn 4002#define maxe 2009000#define oo 0x3F3F3F3Ftypedef struct node{ int x, y;} node;node p[maxn + 3];int edge[maxn + 30];int next[maxe * 2 + 30];int point[maxe * 2 + 30];int cost[maxe * 2 + 30];int src, sink, v, phi, mincost, etop = 1;char used[maxn + 32];int queue[maxn + 32];int dist[maxn + 32];int spfa(){ memset(used, 0, sizeof(used)), memset(dist, 0x3F, sizeof(dist)), dist[src] = 0; memset(queue, 0, sizeof(queue)); int head = src, tail = src, i, e, dtmp; for (; head; head = queue[i], queue[i] = 0) for (used[i = head] = 0, e = edge[i]; e; e = next[e]) if ((dtmp = dist[i] + cost[e]) < dist[point[e]]) { dist[point[e]] = dtmp; if (!used[point[e]]) used[point[e]] = 1, queue[tail] = point[e], tail = point[e]; } return dist[sink];}void add(int s, int t, int w){ point[++etop] = t; next[etop] = edge[s]; edge[s] = etop; cost[etop] = w;}int cmper(const void *i, const void *j){ node *p = (node *)i, *q = (node *)j; if (p->x == q->x) return p->y - q->y; return p->x - q->x;}int lis(int n){ memset(used, 0, sizeof(used)); int tail = 0, j = 1, i; for (; j <= n; ++j) { int l = 0, r = tail, mid = (r + 1) >> 1; for(; l != r; mid = (l + r + 1) >> 1) (p[queue[mid]].y <= p[j].y) ? l = mid : (r = mid - 1); dist[j] = queue[r], queue[++r] = j, r > tail ? ++tail : 0; } for (j = queue[tail], tail = 0; j; j = dist[j]) used[j] = 1, ++tail; for (i = 1; i <= n; ++i) for (j = i + 1; j <= n; ++j) if (used[i] && used[j] && dist[j] == i) add(j, i + n, 0); else if ((!used[i] || !used[j]) && p[i].y <= p[j].y) add(i + n, j, 0); for (i = 1; i <= n; ++i) if (!used[i]) add(src, i, 0), add(i, sink, 0), add(i, i + n, -1); else add(i + n, i, 1); return tail;}int main(int argc, char *argv[]){ FILE *fin = fopen("pacman.in", "r"); FILE *fout = fopen("pacman.out", "w"); int n, i, j, k, ans = 0; fscanf(fin, "%d", &n), src = 2 * n + 1, sink = src + 1, v = sink; for (i = 1; i <= n; ++i) fscanf(fin, "%d%d", &p[i].x, &p[i].y); qsort(p + 1, n, sizeof(p[0]), cmper), ans = lis(n); fprintf(fout, "%d/n", ans - spfa() + 1); fclose(fin); fclose(fout); return 0;}