前言
这是卡 spfa???
建图
首先我们可以想到的就是题意建模,也就是每个人向能做的物品连边。(这种连边方法是错的就不详细说了)
但是我们会发现一点:每人所造的物品在图中会变得不已知。
所以在这里要使用建模中经常会遇到的方法:拆点。
把代表人的点拆成成入点和出点,在出点流向的点的流量设为 inf,源点向入点连流量为 inf 边。让后再让入点连向出点,边数为 S i + 1 S_{i}+1 Si+1,每条边的流量就为 T i , j − T i , j − 1 T_{i,j}-T_{i,j-1} Ti,j−Ti,j−1,费用就为 w i , j w_{i,j} wi,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 ;
}