P2488 [SDOI2011] 工作安排

前言

这是卡 spfa？？？

建图

首先我们可以想到的就是题意建模，也就是每个人向能做的物品连边。（这种连边方法是错的就不详细说了）

但是我们会发现一点：每人所造的物品在图中会变得不已知。

所以在这里要使用建模中经常会遇到的方法：拆点。

把代表人的点拆成成入点和出点，在出点流向的点的流量设为 inf，源点向入点连流量为 inf 边。让后再让入点连向出点，边数为 $S_i+1$，每条边的流量就为 $T_{i,j}-T_{i,j-1}$，费用就为 $w_{i,j}$（只有这种边有费用其余均没有）。（因为在数据中愤怒值递增不用怕 spfa 跑错边的顺序）

最后将所有代表物品的点连向汇点流量为所需个数，这个题就做完了。

但真的做完了吗？

测评结果告诉我们：No。

这样只能得 85 分（如果常数差可能为 75 分）因为边数过多。

应该怎么解决？

我们都知道，一次费用流的 dinic 只会增广源点到汇点的最短路，所以我们只需要让他们只跑这一条路，就不会 TLE。

代码

#include<bits/stdc++.h>
#define int long long 
#define inf 21474836471
#define N 50121
using namespace std ;
int m , n , head[N] , now[N] , dep[N] , s , t , cnt=0 , money=0 , S[N] , to[N] , go[N] ;
bool f[N] ;
struct node{
	int to , next , w , c , I ;
}e[1123456];
void use(int u,int v,int w,int c){
	e[cnt].I = u ;
	e[cnt].to = v ;
	e[cnt].w = w ;
	e[cnt].next = head[u] ;
	e[cnt].c = c ;
	head[u] = cnt++ ;
}
void newnet(int u,int v,int w,int c){use(u,v,w,c) , use(v,u,0,-c) ;}
bool bfs(){
	for(int i=1;i<=t;i++) dep[i] = inf , now[i] = head[i] , f[i] = 0 ;
	queue<int> q ;
	q.push(s) ;
	dep[s] = 0 ;
	while(!q.empty()){
		int u=q.front() ;
		f[u] = 0 ;
		q.pop() ;
		for(int i=head[u];~i;i=e[i].next){
			int x=e[i].to ;
			if(e[i].w&&dep[u]+e[i].c<dep[x]){
				to[x] = i ;
				dep[x] = dep[u]+e[i].c ;
				if(!f[x]){
					f[x] = 1 ;
					q.push(x) ;
				}
			}
		}
	}
	return dep[t]!=inf ;
}
int dfs(int u,int sum){
	if(u==t) return sum ;
	f[u] = 1 ;
	int use=0 ;
	while(1){
		int i = go[u] ;
		int x=e[i].to ;
		if(dep[u]+e[i].c!=dep[x]||!e[i].w||f[x]) continue ;
		int temp=dfs(x,min(sum,e[i].w)) ;
		if(!temp) continue ;
		e[i].w -= temp ;
		e[i^1].w += temp ;
		money += temp*e[i].c ;
		sum -= temp ;
		use += temp ;
		break ;
	}
	return use ;
}
void dinic(){
	while(bfs()){
		int temp=t ;
		while(temp!=s){
			go[e[to[temp]].I] = to[temp] ;
			temp = e[to[temp]].I ;
		}
		dfs(s,inf) ;
	}
	cout << money ;
}
int read(){
	int x=0 , f=1 ;
	char a=getchar() ;
	while(!(a>='0'&&a<='9')){
		if(a=='-') f = -f ; 
		a = getchar() ;
	}
	while(a>='0'&&a<='9'){
		x  *= 10 ;
		x += a-'0' ;
		a = getchar() ;
	}
	return x*f ;
}
signed main(){
	memset(head,-1,sizeof head) ;
	cin >> m >> n ;
	s = m*2+n+1 ;
	t = s+1 ;
	for(int i=1;i<=n;i++){
		int w=read() ;
		newnet(i+m*2,t,w,0) ;
	}
	for(int i=1;i<=m;i++){
		newnet(s,i,inf,0) ;
		for(int j=1;j<=n;j++){
			int temp=read() ;
			if(temp) newnet(i+m,j+m*2,inf,0) ;
		}
	}
	for(int i=1;i<=m;i++){
		int temp1=read() ;
		S[temp1+1] = inf ;
		for(int j=1;j<=temp1;j++)
			S[j] = read() ;
		for(int j=1;j<=temp1+1;j++){
			int temp=read() ;
			newnet(i,i+m,S[j]-S[j-1],temp) ;
		}
	}
	dinic() ;
	return 0 ;
}