BZOJ4346 [POI2016]Nadajniki

我们考虑树形dp，对每个子树求出覆盖这个子树需要的最小代价

首先发现一个点上放多于2个路由器没卵子意义

那么我们考虑一下状态，我们需要用f[i][j][k][l]表示以i为根的子树，i这个点放了j个路由器，i的儿子里有k个路由器（如果大于2个则视为2个），i的父亲至少需要有l个路由器的最小代价，然后挨个子树大力转移即可，不过这样的话复杂度是n*（3^6）的，爆炸了……

然后我们就需要压缩一下状态……首先发现当i不等于0的时候k和l是没有意义的，然后l是一定不会超过2-k的

这样的话有意义的状态就只剩下8种，复杂度变成n*(8^2)，就可以过了

然后尝试着在代码层面对状态进行一些优化，结果还是TLE了……

生气了，手写转移矩阵，加读入优化，define min，成功get rank1

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<bitset>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define MAXN 200010
#define MAXM 1010
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
#define min(a,b) (a)<(b)?(a):(b)
char xch,xB[1<<15],*xS=xB,*xTT=xB;
#define getc() (xS==xTT&&(xTT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xTT)?0:*xS++)
inline int read()
{
	int x=0,f=1;char ch=getc();
	while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getc();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();}
	return x*f;
}
struct vec{
	int to;
	int fro;
};
int n;
vec mp[MAXN*2];
int tai[MAXN],cnt;
const int D=8;
int f[MAXN][D];
int t[D];
//                	 0  1  2  3  4  5  6  7  
int trans[D][D]={	{0 ,0 ,0 ,0 ,0 ,0 ,-1,0 },//0
					{1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 },//1
					{3 ,4 ,6 ,5 ,2 ,-1,-1,-1},//2
					{4 ,4 ,7 ,3 ,3 ,-1,-1,-1},//3
					{4 ,4 ,4 ,4 ,4 ,-1,-1,-1},//4
					{3 ,4 ,6 ,5 ,5 ,-1,-1,-1},//5
					{7 ,4 ,6 ,6 ,6 ,-1,-1,-1},//6
					{4 ,4 ,7 ,7 ,7 ,-1,-1,-1}//7
				};
inline void be(int x,int y){
	mp[++cnt].to=y;
	mp[cnt].fro=tai[x];
	tai[x]=cnt;
}
inline void bde(int x,int y){
	be(x,y);
	be(y,x);
}
void dfs(int x,int F){
	int i,j,k,I,J,K,ii,y;
	memset(f[x],0x3f,sizeof(f[x]));
	f[x][0]=1;
	f[x][1]=2;
	f[x][2]=0;
	for(ii=tai[x];ii;ii=mp[ii].fro){
		y=mp[ii].to;
		if(y!=F){
			dfs(y,x);
			memset(t,0x3f,sizeof(t));
			for(i=0;i<D;i++){
				for(j=0;j<D;j++){
					if(trans[i][j]!=-1){
						t[trans[i][j]]=min(t[trans[i][j]],f[x][i]+f[y][j]);
					}
				}
			}
			memcpy(f[x],t,sizeof(t));
		}
	}
}
int main(){
	int i,x,y;
	n=read();
	for(i=1;i<n;i++){
		x=read();
		y=read();
		bde(x,y);
	}
	dfs(1,0);
	int ans=INF;
	for(i=0;i<=4;i++){
		ans=min(ans,f[1][i]);
	}
	printf("%d\n",ans);
	return 0;
}

/*
7
1 2
2 3
4 3
5 4
6 3
7 6

*/