BZOJ1603 NOI2008 設計路線 樹形dp

題目大意:給定一個根節點爲1的樹,要求在樹中找到一些不相交的鏈,使得每個節點的不便利值得最大值最小,並求出滿足條件的方案個數。一個節點的不便利值就是從該節點到根的路徑上經過的非鏈邊的條數。


題目分析:對於無解的情況是很好處理的,對於題目意思,只要題目所給的是一個森林便可以直接判斷無解,即m<n-1。而這題的關鍵點在於對狀態的分析和設計。事實上因爲n<=100000,所以第一問答案不超過十種(通過滿三叉樹證明)。而利用樹剖也可證明是<logN的,這也已經可以實現了。

考慮f[i][j][k]爲節點i,不便利值<j,與子節點連了k條邊{k=0,1,2}

借將狼踩盡的圖- 這是由狀態直接得到的轉移;


但是顯然會T所以化簡


細節上注意模Q==0時應當保留Q否則會被當做沒方案從而WA掉。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
#include<cassert>
#include<climits>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define RepD(i,n) for(int i=n;i>=0;i--)
#define MEM(a) memset(a,0,sizeof(a))
#define MEMI(a) memset(a,127,sizeof(a))
#define MEMi(a) memset(a,128,sizeof(a))
#define INF (2139062143)
#define phiF (1000000006)
#define MAXN (1000000+10)
typedef long long ll;
struct info{
	int to,next;
}e[200005];
int first[100005],n,m,Q,tot,x,y;ll f[100005][11][3];
void add(int x,int y){
	e[++tot].to=y;
	e[tot].next=first[x];
	first[x]=tot;
}
int get(ll t){return !t?0:t%Q?t%Q:Q;}
void dfs(int u,int fa){
	int cnt=0;
	for (int p=first[u];p;p=e[p].next){
		if (e[p].to!=fa) {
		dfs(e[p].to,u);
		cnt++;}
	}
	Fork(i,0,10) f[u][i][0]=1;
	if (!cnt) return ;
	for (int p=first[u];p;p=e[p].next){
		int v=e[p].to;
		if(v!=fa){
			Fork(j,0,10){
				ll f1;
				if (j){
					f1=f[v][j-1][0]+f[v][j-1][1]+f[v][j-1][2];}else f1=0;
				ll	f2=f[v][j][1]+f[v][j][0];
				ll t;
				t=(ll)f[u][j][2]*f1+(ll)f[u][j][1]*f2;f[u][j][2]=get(t);
				t=(ll)f[u][j][1]*f1+(ll)f[u][j][0]*f2;f[u][j][1]=get(t);
				t=(ll)f[u][j][0]*f1;f[u][j][0]=get(t);

				
			} 
			
		}
		
	}
	
	
}


int main(){
	scanf("%d%d%d",&n,&m,&Q);
	if (m<n-1){
		printf("-1\n-1");
		return 0;
	}
	For (i,m){
		scanf("%d%d",&x,&y);
		add(x,y);
		add(y,x);
	}
	
	dfs(1,0);
	
	Fork(i,0,10)
	if (f[1][i][0]+f[1][i][1]+f[1][i][2]>0) {
		printf("%d\n",i);
		printf("%d",(f[1][i][0]+f[1][i][1]+f[1][i][2])%Q) ;
		return 0;
	}
} 



發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章