2019.03.16【NOIP提高组】模拟B组 状压枚举+贪心+DP+网络流

0 病毒

D(1<=D<=15)种病毒正袭击着牧场里的奶牛，这给牧场主约翰造成了极大的麻烦，因为按照卫生检疫条例的规定：如果牛奶中的病毒超过K(1<=K<=D)种则不允许被销售，约翰的牧场里共有N(1<=N<=1,000)头奶牛，每头奶牛都带有若干种病毒，约翰挤出的牛奶是混放在一起的，要求约翰最多可以挤多少头奶牛的牛奶，使得混放在一起的牛奶不包含超过K种的病毒。

由病毒的种数可得应用状态压缩
枚举状态，如状态满足病毒种数（即状态中1数量）不大于k，则用满足该种状态的奶牛总数更新答案

#include <cstdio>
#include <algorithm>

using namespace std;

int n,d,k,ans;
int a[1005],c[16];
int f[80000];

bool bz(int x){
	int g=0;
	for (int i=1;i<=d;i++)
		if (x&c[i]) g++;
	if (g>k) return 0;
	return 1;
}

int main(){
	scanf("%d%d%d",&n,&d,&k);
	int maxx=1;
	if (d==k) {printf("%d",n);return 0;}
	c[1]=1;
	for (int i=2;i<=d;i++) {
		c[i]=c[i-1]<<1;
		maxx+=c[i];
	}
	for (int i=1;i<=n;i++){
		int v,d;
		scanf("%d",&v);
		for (int j=1;j<=v;j++){
			scanf("%d",&d);
			a[i]|=c[d];
		}
	}
	for (int i=1;i<=maxx;i++)
	if (bz(i)){
		int res=0;
		for (int j=1;j<=n;j++)
			if ((a[j]&i)==a[j]) 
				res++;
		ans=max(ans,res);
	}
	printf("%d",ans);
}

1 奶酪厂

奶牛买了一个奶酪厂生产奶酪,已知每周生产一单位奶酪的费用为C_i，每周可以生产任意数量的奶酪，现在要为接下来N(1<=N<=10,000)周做生产计划。
　　厂里有一个仓库，存储量无穷大，可以用来存储暂时不用的奶酪，每单位奶酪每周花费S(1<=S<=100)。
　　告诉你每周客户的需求量Y_i(0<=Y_i<=10,000)，请你帮忙用最少的钱满足这些需求。
　　
易得每周的奶酪全都选单位费用最小的来制作，那么可以直接做（或单调栈/队列优化）

#include <cstdio>
#include <algorithm>

using namespace std;

int n,s;
int c[10004],y[10004];
int f[10004][3],h=1,t=1;
long long ans;

int main(){
	scanf("%d%d",&n,&s);
	f[1][0]=1000;
	for (int i=1;i<=n;i++){
		scanf("%d%d",&c[i],&y[i]);
		int k=c[i];
		for (int j=i-1;j>=1;j--)
			k=min(k,c[j]+(i-j)*s);
		ans=(long long)ans+k*y[i];		
	}
	printf("%lld",ans);
}

2 太空电梯

奶牛们想用K(1<=K<=400)中石块制造一个太空电梯去太空旅行，每种石块有自己的高度h_i(1<=h_i<=100)和数量c_i(1<=c_i<=10),为了避免宇宙射线的干扰，每种石块不能超过最高可以达到的高度a_i(1<=a_i<=40000)。
　　帮助奶牛用石块堆积一个最高的太空电梯。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

int k;
struct node{
	int h,a,c;
}s[405];
int f[40005],fz[40005];

bool comp(node a,node b){
	return a.a<b.a;
}

int main(){
	scanf("%d",&k);
	for (int i=1;i<=k;i++){
		scanf("%d%d%d",&s[i].h,&s[i].a,&s[i].c);
	}
	f[0]=1;
	sort(s+1,s+1+k,comp);
	for (int i=1;i<=k;i++){
		memcpy(fz,f,sizeof(int)*40003);
		for (int l=s[i].c;l>=1;l--)
			for (int j=s[i].a;j>=s[i].h*l;j--)
			if (fz[j-l*s[i].h]){			
				f[j]|=fz[j-l*s[i].h];
			}
	}
	for (int i=40000;i>=0;i--){
		if (f[i]) {
			printf("%d",i);
			break;
		}
	}
}

3 躲雨

FJ的奶牛真的很怕雨，所以一下雨他们就急着要找地方躲雨。为了做好防雨准备，他们装了警报器，一有雨就拉响警报器，他们利用天气预报来预测雨的到来，但有时候不准导致拉错警报，为了减少拉错警报的可能，他们想在保证奶牛有足够时间躲起来的前提下尽可能晚地拉响警报。

农场有F(1<=F<=200)个奶牛吃草的区域，有P(1<=P<=1500)条路连接一些区域，路是双向的。有些区域有遮雨棚，每个遮雨棚有自己的容量限制。

现在要你计算所有牛都能找到地方躲雨最少需要多少时间。

二分答案，网络流判定
连边方式：新建区域的个数个点，连向汇点，流量为牛棚数；原来的点连向源点，流量为牛数。Floyd算出点与点间最短路，每次判断都重新建图，原来的点和新点间的最短路小于最大值则连边，流量设为大值。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

int m,p,inf=0x3f3f3f,n;
long long ans=-1;
int tot=1,ls[405],ne[100005],y[100005],f[100005];
int v[100010],d[405],cur[405];
long long a[405][405];
int w[405],o[405];

bool bfs(){
	int h=0,t=1;
	memset(d,0,sizeof d);
	memset(v,0,sizeof v);
	v[1]=0;d[0]=1;
	while (h<t){
		int u=ls[v[++h]];
		while (u){
			int xx=v[h];
			if ((d[y[u]]==0)&&(f[u]>0)){
				d[y[u]]=d[xx]+1;
				v[++t]=y[u];
			}
			u=ne[u];
		}
	}
	if (d[2*m+1]>0) return 1;
	return 0;
}

long long dfs(int x,int maxf){
	if ((x==2*m+1)) return maxf;
	for (int &i=cur[x];i;i=ne[i]){
		if ((d[y[i]]==d[x]+1)&&(f[i]>0)){
			long long ee=dfs(y[i],min(f[i],maxf));
			if (ee>0){
				f[i]-=ee;
				f[i^1]+=ee;
				return ee;
			}
		}
	}
	return 0;
}

int dicnic(){
	long long aad=0;
	while (bfs()){
		for (int i=0;i<=2*m+1;i++) cur[i]=ls[i];
		while (long long cd=dfs(0,inf))
			aad+=cd;
	}
	return aad;
}

void set(int u,int v,int z){
	ne[++tot]=ls[u];ls[u]=tot;y[tot]=v;f[tot]=z;
	ne[++tot]=ls[v];ls[v]=tot;y[tot]=u;f[tot]=0;
}

void built(long long mid){
	for (int i=1;i<=m;i++){
		if (w[i]) set(0,i,w[i]);
		if (o[i]) set(i+m,2*m+1,o[i]);
	}
	for (int i=1;i<=m;i++){
		if (w[i]||o[i])
			for (int j=1;j<=m;j++)
			if ((w[j]||o[j])&&(a[i][j]<=mid)){
				set(i,j+m,inf);
		}
	}	
}

int main(){
	scanf("%d%d",&m,&p);
	for (int i=0;i<=m;i++)
		for (int j=0;j<=m;j++) 
			a[i][j]=1e12+1;
	for (int i=1;i<=m;i++){
		scanf("%d%d",&w[i],&o[i]);
		a[i][i]=a[i+m][i+m]=0;
		n+=w[i];
	}	
	for (int i=1;i<=p;i++){
		int u,v,z;
		scanf("%d%d%d",&u,&v,&z);
		if (z<a[u][v]) a[u][v]=z,a[v][u]=z;
	}
	for (int i=1;i<=m;i++)
		for (int j=1;j<=m;j++)
			for (int k=1;k<=m;k++)
				if (a[i][k]+a[k][j]<a[i][j])
					a[i][j]=a[j][i]=a[i][k]+a[k][j];
	long long l=0,r=1e12;
	while (l<=r){
		long long mid=(l+r)>>1;
		tot=1;
		memset(ls,0,sizeof ls);
		built(mid);	
		if (dicnic()==n) 	{ans=mid;r=mid-1;}
				else l=mid+1;
	}
	printf("%lld",ans);
}

好好学习天天向上