POJ 3680 Intervals
题意:给定一些区间,每个区间有一个权值,要求用这些区间去覆盖,每个点最多覆盖k次,问最多得到权值多少
思路:典型的区间k覆盖问题,区间连边容量1,代价-w,然后其他点相邻两两连边,容量k,代价0,跑一下费用流即可
代码:
#include <cstdio> #include <cstring> #include <vector> #include <queue> #include <algorithm> using namespace std; const int MAXNODE = 100005; const int MAXEDGE = 2005; typedef int Type; const Type INF = 0x3f3f3f3f; struct Edge { int u, v; Type cap, flow, cost; Edge() {} Edge(int u, int v, Type cap, Type flow, Type cost) { this->u = u; this->v = v; this->cap = cap; this->flow = flow; this->cost = cost; } }; struct MCFC { int n, m, s, t; Edge edges[MAXEDGE]; int first[MAXNODE]; int next[MAXEDGE]; int inq[MAXNODE]; Type d[MAXNODE]; int p[MAXNODE]; Type a[MAXNODE]; void init(int n) { this->n = n; memset(first, -1, sizeof(first)); m = 0; } void add_Edge(int u, int v, Type cap, Type cost) { edges[m] = Edge(u, v, cap, 0, cost); next[m] = first[u]; first[u] = m++; edges[m] = Edge(v, u, 0, 0, -cost); next[m] = first[v]; first[v] = m++; } bool bellmanford(int s, int t, Type &flow, Type &cost) { for (int i = 0; i < n; i++) d[i] = INF; memset(inq, false, sizeof(inq)); d[s] = 0; inq[s] = true; p[s] = s; a[s] = INF; queue<int> Q; Q.push(s); while (!Q.empty()) { int u = Q.front(); Q.pop(); inq[u] = false; for (int i = first[u]; i != -1; i = next[i]) { Edge& e = edges[i]; if (e.cap > e.flow && d[e.v] > d[u] + e.cost) { d[e.v] = d[u] + e.cost; p[e.v] = i; a[e.v] = min(a[u], e.cap - e.flow); if (!inq[e.v]) {Q.push(e.v); inq[e.v] = true;} } } } if (d[t] == INF) return false; flow += a[t]; cost += d[t] * a[t]; int u = t; while (u != s) { edges[p[u]].flow += a[t]; edges[p[u]^1].flow -= a[t]; u = edges[p[u]].u; } return true; } Type Mincost(int s, int t) { Type flow = 0, cost = 0; while (bellmanford(s, t, flow, cost)); return cost; } } gao; const int N = 405; int t, n, k, save[N], sn; int main() { scanf("%d", &t); while (t--) { scanf("%d%d", &n, &k); int l, r, w; sn = 0; gao.init(100005); while (n--) { scanf("%d%d%d", &l, &r, &w); save[sn++] = l; save[sn++] = r; gao.add_Edge(l, r, 1, -w); } save[sn++] = 0; save[sn++] = 100001; sort(save, save + sn); for (int i = 1; i < sn; i++) gao.add_Edge(save[i - 1], save[i], k, 0); printf("%d\n", -gao.Mincost(0, 100001)); } return 0; }
时间: 2024-12-29 17:40:01