题解 P1264 【K-联赛】
看一个队是否有希望夺冠,即这个队伍将所有能赢的比赛都赢下来,算出这个队最多拿多少分,再判断其他队的积分是否有可能都低于这个最高分。
将所有球队作为图中结点,每两队比赛也作为图中结点。从S引一条边到比赛结点,容量为两队还剩的比赛数,从每队结点连一条边到T,容量为这个队和当前最高分的差,如果一个比赛结点代表的是A和B的比赛,那么从这个结点连两条边分别到A的结点和B的结点 容量为INF
如果这个网络的最大流等于所有未进行的比赛场次之和,则当前球队可以成为冠军
把每个队枚举一遍即可
#include <bits/stdc++.h>
using namespace std;
const int maxn = 700;
const int INF = 1000000000;
struct Edge {
int from, to, cap, flow;
};
bool operator < (const Edge& a, const Edge& b) {
return a.from < b.from || (a.from == b.from && a.to < b.to);
}
struct Dinic {
int n, m, s, t;
vector<Edge> edges; // 边数的两倍
vector<int> G[maxn]; // 邻接表,G[i][j]表示结点i的第j条边在e数组中的序号
bool vis[maxn]; // BFS使用
int d[maxn]; // 从起点到i的距离
int cur[maxn]; // 当前弧指针
void init(int n) {
for(int i = 0; i < n; i++) G[i].clear();
edges.clear();
}
void AddEdge(int from, int to, int cap) {
edges.push_back((Edge){from, to, cap, 0});
edges.push_back((Edge){to, from, 0, 0});
m = edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool BFS() {
memset(vis, 0, sizeof(vis));
queue<int> Q;
Q.push(s);
vis[s] = 1;
d[s] = 0;
while(!Q.empty()) {
int x = Q.front(); Q.pop();
for(int i = 0; i < G[x].size(); i++) {
Edge& e = edges[G[x][i]];
if(!vis[e.to] && e.cap > e.flow) {
vis[e.to] = 1;
d[e.to] = d[x] + 1;
Q.push(e.to);
}
}
}
return vis[t];
}
int DFS(int x, int a) {
if(x == t || a == 0) return a;
int flow = 0, f;
for(int& i = cur[x]; i < G[x].size(); i++) {
Edge& e = edges[G[x][i]];
if(d[x] + 1 == d[e.to] && (f = DFS(e.to, min(a, e.cap-e.flow))) > 0) {
e.flow += f;
edges[G[x][i]^1].flow -= f;
flow += f;
a -= f;
if(a == 0) break;
}
}
return flow;
}
int Maxflow(int s, int t) {
this->s = s; this->t = t;
int flow = 0;
while(BFS()) {
memset(cur, 0, sizeof(cur));
flow += DFS(s, INF);
}
return flow;
}
};
Dinic g;
const int maxt = 25 + 5;
int n, w[maxt], d[maxt], a[maxt][maxt];
inline int ID(int u, int v) { return u*n+v+1; }
inline int ID(int u) { return n*n+u+1; }
bool canWin(int team) {
// 计算team全胜后的总胜利场数
int total = w[team];
for(int i = 0; i < n; i++)
total += a[team][i];
for(int i = 0; i < n; i++)
if(w[i] > total) return false;
// 构图。s=0, 结点(u,v)的编号为u*n+v+1, 结点u的编号为n^2+u+1, t=n^2+n+1
g.init(n*n+n+2);
int full = 0;
int s = 0, t = n*n+n+1;
for(int u = 0; u < n; u++) {
for(int v = u+1; v < n; v++) {
if(a[u][v] > 0) g.AddEdge(s, ID(u,v), a[u][v]); // S到(u,v)的弧
full += a[u][v];
g.AddEdge(ID(u,v), ID(u), INF); // (u,v)到u的弧
g.AddEdge(ID(u,v), ID(v), INF); // (u,v)到v的弧
}
if(w[u] < total) g.AddEdge(ID(u), t, total-w[u]); // u到T的弧
}
return g.Maxflow(s, t) == full;
}
int main() {
scanf("%d", &n);
for(int i = 0; i < n; i++) scanf("%d%d", &w[i], &d[i]);
for(int i = 0; i < n; i++)
for(int j = 0; j < n; j++)
scanf("%d", &a[i][j]);
bool first = true;
for(int i = 0; i < n; i++)
if(canWin(i)) {
if(first) first = false; else printf(" ");
printf("%d", i+1);
}
printf("\n");
return 0;
}