NYIST 489 哭泣天使

哭泣天使
时间限制：1000 ms | 内存限制：65535 KB
难度：5

描述
Doctor Who乘着Tardis带着Amy来到了一个星球，一开Tadis大门，发现这个星球上有个壮观的石像群，全是一些天使石像，有的石像在哭泣，有的石像像在微笑，共有m行n列，Doctor用“音速起子”扫描了一下整个石像群，得到了每行天使中在哭泣的天使的个数。当他与Amy在这里行走了一段时间之后，Doctor忽然想起了什么，怀疑这些石像是不是传说中的一种黑暗生物——“哭泣天使”——一种看似石像，却会在人不看它的时候移动，会强制把人送回某个过去的时间点，并借此汲取时间能量的生物。Doctor可不想自己和Amy迷失在一个未知的时间点里，于是Doctor立刻用“音速起子”又扫描了整个石像群，想再看看每行的在哭泣的天使个数与刚才是否相符，但是，越急就越容易出错，他一不小心扫描错了，扫描出了每列中哭泣的天使的个数。现在，由于音速起子的能量不足了，他不能够再次扫描，他想根据已有的数据判断出是否有天使改变了自己的表情，从哭泣变成不哭泣或者从不哭泣变成哭泣了。

输入
第一行是一个整数T，表示共有T组测试数据(T<=50)
每组测试数据第一行是两个整数m,n(0<m,n<=300)分别表示行数和列数
随后的两行分别有m个数和n个数分别表示对应m行中哭泣的天使石像的个数与对应n个列中哭泣的天使石像的个数。

输出
如果能根据已有信息判断出必然有石像改变了表情，则输出Terrible
如果根据已有信息无法确定石像发生了改变，则输出Not Sure (有时，你确定两次扫描时状态相同，但由于不确定之间是否发生过改变，故也输出Not Sure)

样例输入
2
2 3
1 1
1 1 0
3 3
0 1 2
3 0 0

样例输出
Not Sure
Terrible

来源
[张云聪]原创

上传者
张云聪

 

解题：最大流，居然卡dinic，把源点汇点反转，建立逆图，果断AC了。

 1 #include <iostream>

 2 #include <cstdio>

 3 #include <cstring>

 4 #include <cmath>

 5 #include <algorithm>

 6 #include <climits>

 7 #include <vector>

 8 #include <queue>

 9 #include <cstdlib>

10 #include <string>

11 #include <set>

12 #include <stack>

13 #define LL long long

14 #define pii pair<int,int>

15 #define INF 0x3f3f3f3f

16 using namespace std;

17 const int maxn = 810;

18 struct arc {

19     int to,flow,next;

20     arc(int x = 0,int y = 0,int z = -1) {

21         to = x;

22         flow = y;

23         next = z;

24     }

25 };

26 arc e[maxn*maxn];

27 int head[maxn],d[maxn],cur[maxn];

28 int tot,n,m,S,T;

29 void add(int u,int v,int flow) {

30     e[tot] = arc(v,flow,head[u]);

31     head[u] = tot++;

32     e[tot] = arc(u,0,head[v]);

33     head[v] = tot++;

34 }

35 bool bfs() {

36     queue<int>q;

37     memset(d,-1,sizeof(d));

38     d[T] = 1;

39     q.push(T);

40     while(!q.empty()) {

41         int u = q.front();

42         q.pop();

43         for(int i = head[u]; ~i; i = e[i].next) {

44             if(e[i].flow && d[e[i].to] == -1) {

45                 d[e[i].to] = d[u] + 1;

46                 q.push(e[i].to);

47             }

48         }

49     }

50     return d[S] > -1;

51 }

52 int dfs(int u,int low) {

53     if(u == S) return low;

54     int tmp = 0,a;

55     for(int &i = cur[u]; ~i; i = e[i].next) {

56         if(e[i].flow > 0 && d[e[i].to] == d[u] + 1&&(a=dfs(e[i].to,min(low,e[i].flow)))) {

57             e[i].flow -= a;

58             e[i^1].flow += a;

59             tmp += a;

60             low -= a;

61             if(!low) break;

62         }

63     }

64     if(!tmp) d[u] = -1;

65     return tmp;

66 }

67 int dinic(int sum) {

68     int ans = 0;

69     while(bfs()) {

70         memcpy(cur,head,sizeof(head));

71         ans += dfs(T,INF);

72         if(ans > sum) break;

73     }

74     return ans;

75 }

76 int main() {

77     int tmp,sum,t,sum2;

78     scanf("%d",&t);

79     while(t--) {

80         scanf("%d %d",&n,&m);

81         memset(head,-1,sizeof(head));

82         sum2 = sum = S = tot = 0;

83         T = n+m+1;

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

85             scanf("%d",&tmp);

86             add(i,S,tmp);

87             sum += tmp;

88             for(int j = n+1; j <= n+m; ++j) add(j,i,1);

89         }

90         for(int i = 1; i <= m; ++i){

91             scanf("%d",&tmp);

92             add(T,n+i,tmp);

93             sum2 += tmp;

94         }

95         printf("%s\n",sum == sum2&&dinic(sum) == sum?"Not Sure":"Terrible");

96     }

97     return 0;

98 }

View Code