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 0x3F3F3F3F typedef 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; }