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);
}