Day3

P4322 [JSOI2016]最佳团体

$01$分数规划与树形$dp$结合

题意：
求$\frac{\sum p_i}{\sum s_i}$的最大值，我们可以二分一个$v$，则我们只需要求$\sum p_i>\sum s_i\ast v==>\sum p_i-\sum s_i\ast v>0$就$ok$

让每一个点的权值等于$p_i-s_i\ast v$，然后做一个树形DP。

树形$dp$就直接用有依赖的背包问题则可以解决

最后一个注意点就是这题有一个虚根，也就是0号节点

Code

#include<bits/stdc++.h>
#define rep(i,a,b) for(int (i)=(a);(i)<=(b);(i)++)
#define ll long long
#define eps 1e-4
using namespace std;
const int maxn=1e5+10;
const int maxm=2500+10;
int n,cnt,m;
double f[maxm][maxm],head[maxn],w[maxn];
int fa[maxn],size[maxn];
double l=0.0,r=10000.0,ans=0;

template <class t> inline void read(t &x) {
	int f=1;x=0;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-') ch=getchar();ch=getchar();}
	while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
	x*=f;
}

struct node{
	double cost,v;
}a[maxn];

struct {
	int v,nex;
}e[maxn];

void add(int u,int v) {
	e[++cnt].v=v;
	e[cnt].nex=head[u];
	head[u]=cnt;
}

void readdata() {
	read(m);read(n);
	m++;n++;
	rep(i,2,n) {
		int fa;
		scanf("%lf%lf",&a[i].cost,&a[i].v);
		read(fa);fa++;
		add(fa,i);
	}
}

void dfs(int u,double mid) {
	f[u][0]=0;f[u][1]=a[u].v-mid*a[u].cost;
	size[u]=1;
	for(int i=head[u];i;i=e[i].nex) {
		int v=e[i].v;
		dfs(v,mid);
		size[u]+=size[v];
		for(int j=size[u];j>=0;j--) {
			for(int k=0;k<=size[v] && k<j;k++) {
				f[u][j]=max(f[u][j],f[u][j-k]+f[v][k]);
			}
		}
	}
}

void work() {
	while(l+eps<=r) {
		double mid=(l+r)/2.0;
		//rep(i,1,n) w[i]=a[i].v-mid*a[i].cost;
		memset(f,-0x3f3f,sizeof(f));
	//	memset(size,false,sizeof(size));
		dfs(1,mid);
		if(f[1][m]>=0) l=mid;
		else r=mid;
	}
	printf("%.3lf",l);
}

int main()  {
	//freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}

P2704 [NOI2001]炮兵阵地

此题和SCOI差不多，只不过要多判几个条件，横竖的相邻两个炮兵之间的距离必须差三，在判断平原和山坡（也就是哪里能放炮兵也需要一定的技巧）

一开始在存图的时候就用二进制数保留下来，山坡为$1$为了和后面方便状态判断

此外状态表示$f[i][s][k]$表示在$k$行时的状态为$s$，上一行的状态为$i$最多放的炮兵的数目

Code

#include<bits/stdc++.h>
#define rep(i,a,b) for(int (i)=(a);(i)<=(b);(i)++)
#define ll long long
using namespace std;
const int maxn=1e5+10;
const int maxm=1<<10;
int sum[maxm],a[maxn];
int f[maxm][maxm][102];
int m,n;

template <class t> inline void read(t &x) {
	int f=1;x=0;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-') ch=getchar();ch=getchar();}
	while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
	x*=f;
}

void pre() {
	rep(i,0,(1<<m)-1) {
		rep(j,0,m-1) {
			if((1<<j)&i) sum[i]++;
		}
	}
}

void readdata() {
	read(n);read(m);
	rep(i,1,n)
	 rep(j,1,m) {
	 	char x;cin>>x;
	 	a[i]<<=1;
	 	if(x=='H') a[i]++;
	 }
	pre();
	rep(s,0,(1<<m)-1) {
		if(((a[1]&s) || ((s<<1)&s) || ((s<<2)&s))) continue;
		f[0][s][1]=sum[s];
	}
}

void work() {
	rep(i,0,(1<<m)-1)
	 rep(j,0,(1<<m)-1) {
	 	if((i&j) || (i&a[1]) || (j&a[2]) || (i&(i<<1)) || (i&(i<<2)) || (j&(j<<1)) || (j&(j<<2))) continue;
	 	f[i][j][2]=sum[i]+sum[j];
	 }
	rep(i,3,n) {
		rep(s,0,(1<<m)-1) {
			if((s&a[i]) || (s&(s<<1)) || (s&(s<<2))) continue;
			rep(j,0,(1<<m)-1) {
				if(((s&j) || (j&(j<<1)) || (j&(j<<2)) || (j&a[i-1]))) continue;
				rep(k,0,(1<<m)-1) {
					if(((k&s) || (k&j) || (k&(k<<1)) || (k&(k<<2)) || (k&a[i-2]))) continue;
					f[j][s][i]=max(f[k][j][i-1]+sum[s],f[j][s][i]);
				}
			}
		}
	}
	int ans=0;
	rep(i,0,(1<<m)-1)
	 rep(j,0,(1<<m)-1) {
	 	ans=max(ans,f[i][j][n]);
	 }
	printf("%d",ans);
}

int main() {
//	freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}

P1273 有线电视网

这是一个树上背包问题
$设f[u][i][j]表示以u为根节点的前i棵子树有j个用户所能获得的最大价值$
由于可以通过背包的方式（是选还是不选）滚动掉一维

for(int k=head[u];k;k=e[k].nex) {
		int v=e[k].v;
		if(v==fa) continue;
		t=dfs(v,u),tot+=t;//t为可选择数，tot为背包容积
		don(j,tot,1) {//由于是01背包的滚动优化，所以倒着来
			rep(i,1,t)
			if(j>=i) f[u][j]=max(f[u][j],f[u][j-i]+f[v][i]-e[k].w);
		}
	}

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(register int (i)=(a);(i)<=(b);(i)++)
#define don(i,a,b) for(register int (i)=(a);(i)>=(b);(i)--)
using namespace std;
const int maxn=1e5+10;
const int maxm=3e3+10;
int m,n,cnt=0;
int head[maxn],w[maxn],sum[maxn],f[maxm][maxm];

template <class t> inline void read(t &x)
{
    x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
    x*=f;
}

struct node{
	int v,nex,w;
}e[maxn<<1];

void add(int u,int v,int w) {
	e[++cnt].v=v;
	e[cnt].w=w;
	e[cnt].nex=head[u];
	head[u]=cnt;
}

void readdata()
{
	read(n),read(m);
	rep(i,1,n-m) {
		int num;
		read(num);
		while(num--) {
			int x,w;
			read(x),read(w);
			add(i,x,w);
			add(x,i,w);
		}
	}
	rep(i,n-m+1,n) {
		read(w[i]);
	}
}

inline int dfs(int u,int fa) {
	if(u>n-m) {//由于题上说明了只有从n-m开始才是用户的编号
		f[u][1]=w[u];//用户单独，只选用户一个就是用户自己
		return 1;
	}
	int t,tot=0;
	for(int k=head[u];k;k=e[k].nex) {
		int v=e[k].v;
		if(v==fa) continue;
		t=dfs(v,u),tot+=t;
		don(j,tot,1) {
			rep(i,1,t)
			if(j>=i) f[u][j]=max(f[u][j],f[u][j-i]+f[v][i]-e[k].w);
		}
	}
	return tot;
}



void work()
{
	int ans=0;
	memset(f,-0x3f,sizeof(f));
	memset(sum,false,sizeof(sum));
	rep(i,1,n) f[i][0]=0;
	dfs(1,0);
	don(i,m,1) {//找最大，如果成立则是最大
		if(f[1][i]>=0) {
			ans=i;
			break;
		}
	}
	printf("%d\n",ans);
}

int main()
{
	//freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}