题意:给定个一个DAG,每条边至少被覆盖一次,且有一个代价,问从起点开始的路径覆盖所有边至少一次的最小代价。
思路:上下界最小费用可行流。
#include <bits/stdc++.h>
using namespace std;
typedef int lint;
typedef long long LL;
const int inf = 0x3f3f3f3f;
const int maxn = 505;
struct EDGE {
int from, to, next, cost, cap; // 如果需要修改 cost为LL
};
namespace MFMC {
const static int maxn = 100011;
const static int maxm = 500005;
EDGE edge[maxm];
int tot, he[maxn], n;
void init(int _n) {
tot = 0;
n = _n + 1;
memset(he, -1, n * sizeof(int));
}
void Add(int u, int v, int cap,int cost) { // 如果需要修改 cost为LL
edge[tot] = EDGE{u, v, he[u], cost, cap};
he[u] = tot++;
}
void add(int u, int v, int cap,int cost) { // 如果需要修改 cost为LL
Add(u, v, cap,cost);
Add(v, u, 0,-cost);
}
//O(VE)
//record_e[i]是fa[i]->i的边的编号
template<typename DT>
void spfa(int s, DT dist[], int rec[]) {
queue<int> q;
static bool inq[maxn];
memset(dist, 0x3f, n * sizeof(DT));
memset(inq, 0, n * sizeof(bool));
memset(rec, -1, n * sizeof(int));
q.push(s);
dist[s] = 0;
while (!q.empty()) {
int u = q.front();
q.pop();
inq[u] = 0;
for (int e = he[u]; ~e; e = edge[e].next) {
if (0 == edge[e].cap)
continue;
int v = edge[e].to;
if (dist[v] > dist[u] + edge[e].cost) {
dist[v] = dist[u] + edge[e].cost;
rec[v] = e;
if (!inq[v]) {
q.push(v);
inq[v] = 1;
}
}
}
}
}
template<typename DT>
void dijkstra(int s, DT dist[], int rec[]) {
priority_queue<pair<DT, int> > q;//-dist, vertex
memset(dist, 0x3f, n * sizeof(DT));
memset(rec, -1, n * sizeof(int));
dist[s] = 0;
q.push(make_pair(0, s));
while (!q.empty()) {
s = q.top().second;
DT c = -q.top().first;
q.pop();
if (c != dist[s]) continue;
for (int e = he[s]; ~e; e = edge[e].next) {
if (0 == edge[e].cap) continue;
int v = edge[e].to;
if (dist[v] > c + edge[e].cost) {
dist[v] = c + edge[e].cost;
rec[v] = e;
q.push(make_pair(-dist[v], v));
}
}
}
}
//Need dijkstra_GRAPH_EDGES_PQ
//O(FE log(E)),F is the maximum flow
template<typename FT, typename CT>
void mfmc(int s, int t, FT &maxflow, CT &mincost) {
CT inf;
memset(&inf, 0x3f, sizeof(CT));
static CT dist[maxn];
static int rec_e[maxn];
maxflow = mincost = 0;
CT realdist = 0; //real distance from s to t
bool first = true;
while (1) {
if (first) {
spfa( s, dist, rec_e);
first = false;
} else {
dijkstra( s, dist, rec_e);
}
if (inf == dist[t])
break;
FT minF = numeric_limits<FT>::max();
for (int e = rec_e[t]; ~e; e = rec_e[edge[e].from])
minF = min(minF, (FT) edge[e].cap);
maxflow += minF;
realdist += dist[t];
mincost += minF * realdist;
for (int e = rec_e[t]; ~e; e = rec_e[edge[e].from]) {
edge[e].cap -= minF;
edge[e ^ 1].cap += minF;
}
for (int e = 0; e < tot; ++e) {
EDGE &ed = edge[e];
ed.cost += dist[ed.from] - dist[ed.to];
}
}
}
};
int in[maxn],ou[maxn];
int main(){
int n;
scanf("%d",&n);
int S = 1,T = n+1,SS = n+2,TT = n+3;
MFMC::init(TT);
int cnt,b,t;
int ans = 0;
for( int i = 1;i <= n;i++ ){
scanf("%d",&cnt);
for( int j = 1;j <= cnt;j++ ){
scanf("%d%d",&b,&t);
ans += t;
MFMC::add( i,b,inf,t );
in[b]++;ou[i]++;
}
}
MFMC::add( T,S,inf,0 );
for( int i = 1;i <= n;i++ ){
MFMC::add( i,T,inf,0 );
if( in[i] > ou[i] ){
MFMC::add( SS,i,in[i]-ou[i],0 );
}
if( ou[i] > in[i] ){
MFMC::add( i,TT,ou[i]-in[i],0 );
}
}
int maxflow,mincost;
MFMC::mfmc( SS,TT,maxflow,mincost );
ans += mincost;
printf("%d\n",ans);
return 0;
}