NOI2008 图形DP DOIday01
Z 国坐落于遥远而又神奇的东方半岛上,在小 Z 的统治时代公路成为这里主要的交通手段。Z 国共有 N 座城市,一些城市之间由双向的公路所连接。非常神奇的是 Z 国的每个城市所处的经度都不相同,并且最多只和一个位于它东边的城市直接通过公路相连。Z 国的首都是 Z 国政治经济文化旅游的中心,每天都有成千上万的人从 Z 国的其他城市涌向首都。
为了使 Z 国的交通更加便利顺畅,小 Z 决定在 Z 国的公路系统中确定若干条规划路线,将其中的公路全部改建为铁路。
我们定义每条规划路线为一个长度大于 1 的城市序列,每个城市在该序列中最多出现一次,序列中相邻的城市之间由公路直接相连(待改建为铁路)。并且,每个城市最多只能出现在一条规划路线中,也就是说,任意两条规划路线不能有公共部分。
当然在一般情况下是不可能将所有的公路修建为铁路的,因此从有些城市出发去往首都依然需要通过乘坐长途汽车,而长途汽车只往返于公路连接的相邻的城市之间,因此从某个城市出发可能需要不断地换乘长途汽车和火车才能到达首都。
我们定义一个城市的“不便利值”为从它出发到首都需要乘坐的长途汽车的次数,而 Z 国的交通系统的“不便利值”为所有城市的不便利值的最大值,很明显首都的“不便利值”为 0。小 Z 想知道如何确定规划路线修建铁路使得 Z 国的交通系统的“不便利值”最小,以及有多少种不同的规划路线的选择方案使得“不便利值”达到最小。当然方案总数可能非常大,小 Z 只关心这个天文数字 mod Q 后的值。
注意:规划路线 1−2−3 和规划路线 3−2−1 是等价的,即将一条规划路线翻转,依然认为是等价的。两个方案不同当且仅当其中一个方案中存在一条规划路线不属于另一个方案。
输入描述
第一行包含三个正整数 N,M,Q,其中 N 表示城市个数, M 表示公路总数,N 个城市从 1∼N 编号,其中编号为 1 的是首都。Q 表示上文提到的设计路线的方法总数的模数。
接下来 M 行,每行两个不同的整数 ai,bi(1≤ai,bi≤N),表示有一条公路连接城市 ai 和城市 bi。
输入数据保证一条公路只出现一次。
输出描述
应包含两行。
第 1 行为一个整数,表示最小的“不便利值”。
第 2 行为一个整数,表示使“不便利值”达到最小时不同的设计路线的方法总数 mod Q 的值。
如果某个城市无法到达首都,则输出两行 −1。
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL maxn=1e6, lim=10;
struct node{
LL to,next;
}dis[maxn];
LL n,m,p,num;
LL head[maxn],dp[maxn][lim+9][5];
inline void Add(LL u,LL v){
dis[++num]=(node){v,head[u]}, head[u]=num;
}
inline LL Mod(LL x){
return (x-1)%p+1;
}
void Dfs(LL u,LL fa){
for(LL i=head[u];i;i=dis[i].next){
LL v(dis[i].to);
if(v==fa) continue;
Dfs(v,u);
}
for(LL i=0;i<=lim;++i) dp[u][i][0]=1;
for(LL i=head[u];i;i=dis[i].next){
LL v(dis[i].to);
if(v==fa) continue;
dp[u][0][2]=Mod( dp[u][0][1]*(dp[v][0][0]+dp[v][0][1]) ),
dp[u][0][1]=Mod( dp[u][0][0]*(dp[v][0][0]+dp[v][0][1]) ),
dp[u][0][0]=0;
for(LL j=1;j<=lim;++j){
dp[u][j][2]=Mod( Mod(dp[u][j][2]*(dp[v][j-1][0]+dp[v][j-1][1]+dp[v][j-1][2])) + Mod(dp[u][j][1]*(dp[v][j][0]+dp[v][j][1])) ),
dp[u][j][1]=Mod( Mod(dp[u][j][1]*(dp[v][j-1][0]+dp[v][j-1][1]+dp[v][j-1][2])) + Mod(dp[u][j][0]*(dp[v][j][0]+dp[v][j][1])) ),
dp[u][j][0]=Mod( dp[u][j][0]*(dp[v][j-1][0]+dp[v][j-1][1]+dp[v][j-1][2]) );
}
}
}
int main(){
scanf("%lld%lld%lld",&n,&m,&p);
if(m!=n-1){
printf("-1\n-1");
return 0;
}
for(LL i=1;i<=m;++i){
LL u,v; scanf("%lld%lld",&u,&v);
Add(u,v), Add(v,u);
}
Dfs(1,0);
for(LL i=0;i<=lim;++i)
if(dp[1][i][0]+dp[1][i][1]+dp[1][i][2]){
printf("%lld\n%lld",i,(dp[1][i][0]+dp[1][i][1]+dp[1][i][2])%p);
return 0;
}
return 0;
}