QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#504882#9101. Zayin and BusWorldFinalEscaped#TL 0ms0kbC++143.1kb2024-08-04 16:55:352024-08-04 16:55:36

Judging History

你现在查看的是最新测评结果

  • [2024-08-04 16:55:36]
  • 评测
  • 测评结果:TL
  • 用时:0ms
  • 内存:0kb
  • [2024-08-04 16:55:35]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;

#define ll long long

const int N = 105;
const int M = 100005;
const int inf = 1e6;

struct Edge {
    int to, nxt, cap, flow, cost; 
} edge[M];
int head[N], tol;
int pre[N], dis[N];
bool vis[N];

int node, n, m;

void add(int u, int v, int cap, int cost) {
    edge[tol] = {v, head[u], cap, 0, cost};
    head[u] = tol;
    tol++;
}
void addedge(int u, int v, int cap, int cost) {
    add(u, v, cap, cost);
    add(v, u, 0, -cost);    
}
bool spfa(int S, int T) {
    queue<int> q;
    for (int i = 0; i < node; i++) {
        dis[i] = inf;
        vis[i] = false;
        pre[i] = -1;
    }
    dis[S] = 0, vis[S] = true, q.push(S);
    while (!q.empty()) {
        int u = q.front(); q.pop();
        vis[u] = false;
        for (int i = head[u]; i; i = edge[i].nxt) {
            int v = edge[i].to;
            if (edge[i].cap > edge[i].flow && dis[v] > dis[u] + edge[i].cost) {
                dis[v] = dis[u] + edge[i].cost;
                pre[v] = i;
                if (!vis[v]) vis[v] = true, q.push(v);
            }
        }
    }
    if (pre[T] == -1) return false;
    else return true;
}
int mcmf(int S, int T, ll &cost) {
    int flow = 0; cost = 0;
    while (spfa(S, T)) {
        int Min = inf;
        for (int i = pre[i]; i != -1; i = pre[edge[i ^ 1].to]) {
            if (Min > edge[i].cap - edge[i].flow)
                Min = edge[i].cap - edge[i].flow;
        }
        flow += Min;
        for (int i = pre[T]; i != -1; i = pre[edge[i ^ 1].to]) {
            edge[i].flow += Min;
            edge[i ^ 1].flow -= Min;
            cost += edge[i].cost * Min;
        }
    }
    return flow;
}

int a[55], b[55], c[55], id[55];
vector<int> to[55];

void sc() {    
    memset(head, -1, sizeof(head));
    tol = 0;
    node = 0;

    scanf("%d", &n);
    scanf("%d", &m);
    for (int i = 1; i <= m; i++) {
        scanf("%d%d%d", &a[i], &b[i], &c[i]);
        int x; scanf("%d", &x);
        while (x--) {
            int t; scanf("%d", &t);
            to[i].push_back(t);
        }
        id[i] = i;
    }
    sort(id + 1, id + m + 1, [&](int x, int y) {
        int k1 = b[x] % 2 == 0;
        int k2 = b[y] % 2 == 0;
        return k1 != k2 ? k1 > k2 : x < y;
    });
    
    int S = n + m + 1, T = n + m + 2;
    for (int i = 1; i <= n; i++) {
        addedge(S, i, 1, 0);
    }
    int cur = 0;
    for (int _ = 1; _ <= m; _++) {
        int i = id[_];
        for (auto j: to[i])
            addedge(j, n + i, 1, 0);
        
        addedge(i + n, T, a[i], 0);
        for (int j = 1; j <= 2 * b[i]; j++) {
            ++cur;
            addedge(i + n, T, 1, cur * inf + (j & 1 ? 1 : 0));
        }
        addedge(i + n, T, c[i], inf * inf);
    }
    
    cerr << "out!\n";
    
    node = n + m + 2;
    ll costs = 0;
    int flows = mcmf(S, T, costs);
    printf("flows = %d, costs = %lld\n", flows, costs);
}

int main() {
    int T; scanf("%d", &T);
    while (T--) sc();
    return 0;
}

/*
2
5
2
2 3 1 2 3 1 
1 1 1 3 4 5 2
5
2
2 3 1 1 3 
1 0 1 2 1 2
*/

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 0
Time Limit Exceeded

input:

14
1
1
1
2
1
3
2
1
1 1
2
1
1 2
2
1
2 1
2
1
1 3
2
1
3 1
2
1
1 4
2
1
4 1
3
1 1
1 1 1
3
1 2
1 1 1
3
1 1
1 3 2
3
1 2
1 3 2

output:


result: