HDU - 4322

題目鏈接：HDU - 4322

---

因爲不喜歡造成的都是一樣的，所以我們其實是讓喜歡的代價儘量多。

然後每個人向匯點連邊流量爲b/k，費用爲k

但是可能某個人拿糖之後會超過本身的B就變得不優了。但是也可能會繼續拿喜歡的糖果（當糖果剩餘時），所以當b不能整除k的時候，加一個到匯點流量爲1，費用爲b%k的邊。

最後看剩下用1的糖果是否滿足即可。

---

AC代碼：

#pragma GCC optimize("-Ofast","-funroll-all-loops")
#include<bits/stdc++.h>
//#define int long long
using namespace std;
const int inf=0x3f3f3f3f;
const int N=50,M=1e5+10;
int n,m,k,s,t,d[N],st[N],vis[N],b[N],g[N][N],cost,maxflow,ts,sum;
int head[N],nex[M],to[M],w[M],flow[M],tot;
inline void ade(int a,int b,int c,int d){
	to[++tot]=b; nex[tot]=head[a]; w[tot]=d; flow[tot]=c; head[a]=tot;
}
inline void add(int a,int b,int c,int d){ade(a,b,c,d);	ade(b,a,0,-d);}
inline int spfa(){
	queue<int> q; q.push(s); memset(st,0,sizeof st); memset(d,0x3f,sizeof d); d[s]=0;
	while(q.size()){
		int u=q.front();	q.pop();	vis[u]=0;
		for(int i=head[u];i;i=nex[i]){
			if(flow[i]&&d[to[i]]>d[u]+w[i]){
				d[to[i]]=d[u]+w[i];	
				if(!vis[to[i]])	q.push(to[i]),vis[to[i]]=1;
			}
		}
	}
	return d[t]<inf;
}
int dfs(int x,int f){
	if(x==t)	return f;	st[x]=1;	int fl=0;
	for(int i=head[x];i&&f;i=nex[i]){
		if(!st[to[i]]&&flow[i]&&d[to[i]]==d[x]+w[i]){
			int mi=dfs(to[i],min(flow[i],f));
			flow[i]-=mi,flow[i^1]+=mi,fl+=mi,f-=mi;
		}
	}
	return fl;
}
inline void solve(){
	cin>>n>>m>>k; t=n+m+1; tot=1; memset(head,0,sizeof head); cost=maxflow=sum=0;
	for(int i=1;i<=m;i++)	cin>>b[i],sum+=b[i];
	for(int i=1;i<=m;i++)	for(int j=1;j<=n;j++)	cin>>g[i][j];
	for(int i=1;i<=n;i++)	add(s,i,1,0);
	for(int i=1;i<=m;i++)	for(int j=1;j<=n;j++)	if(g[i][j])	add(j,i+n,1,0);
	for(int i=1;i<=m;i++){
		add(i+n,t,b[i]/k,-k);
		if(b[i]%k)	add(i+n,t,1,-b[i]%k);
	}
	while(spfa()){int f=dfs(s,inf);	cost-=f*d[t],maxflow+=f;}
	printf("Case #%d: ",++ts);
	puts((n-maxflow)+cost>=sum?"YES":"NO");
}
signed main(){
	int T; cin>>T; while(T--) solve();
	return 0;
}