[BZOJ]1063 道路设计(Noi2008)

　　省选一试后的第一篇blog！

Description

　　Z国坐落于遥远而又神奇的东方半岛上，在小Z的统治时代，公路成为这里主要的交通手段。Z国共有n座城市，一些城市之间由双向的公路所链接。很是神奇的是Z国的每一个城市所处的经度都不相同，而且最多只和一个位于它东边的城市直接经过公路相连。Z国的首都是Z国政治经济文化旅游的中心，天天都有成千上万的人从Z国的其余城市涌向首都。为了使Z国的交通更加便利顺畅，小Z决定在Z国的公路系统中肯定若干条规划路线，将其中的公路所有改建为铁路。咱们定义每条规划路线为一个长度大于1的城市序列，每一个城市在该序列中最多出现一次，序列中相邻的城市之间由公路直接相连(待改建为铁路)。而且，每一个城市最多只能出如今一条规划路线中，也就是说，任意两条规划路线不能有公共部分。固然在通常状况下是不可能将全部的公路修建为铁路的，所以从有些城市出发去往首都依然须要经过乘坐长途汽车，而长途汽车只往返于公路链接的相邻的城市之间，所以从某个城市出发可能须要不断地换乘长途汽车和火车才能到达首都。咱们定义一个城市的“不便利值”为从它出发到首都须要乘坐的长途汽车的次数，而Z国的交通系统的“不便利值”为全部城市的不便利值的最大值，很明显首都的“不便利值”为0。小Z想知道如何肯定规划路线修建铁路使得Z国的交通系统的“不便利值”最小，以及有多少种不一样的规划路线的选择方案使得“不便利值”达到最小。固然方案总数可能很是大，小Z只关心这个天文数字modQ后的值。注意：规划路线1-2-3和规划路线3-2-1是等价的，即将一条规划路线翻转依然认为是等价的。两个方案不一样当且仅当其中一个方案中存在一条规划路线不属于另外一个方案。

Input

　　第一行包含三个正整数N、M、Q，其中N表示城市个数,M表示公路总数,N个城市从1~N编号,其中编号为1的是首都。Q表示上文提到的设计路线的方法总数的模数。接下来M行，每行两个不一样的正数ai、bi表示有一条公路链接城市ai和城市bi。输入数据保证一条公路只出现一次。

Output

　　包含两行。第一行为一个整数，表示最小的“不便利值”。第二行为一个整数，表示使“不便利值”达到最小时不一样的设计路线的方法总数modQ的值。若是某个城市没法到达首都，则输出两行-1。

Sample Input

　　5 4 100
　　1 2
　　4 5
　　1 3
　　4 1

Sample Output

　　1
　　10

HINT

　　对于100%的数据，知足1≤N,M≤100000，1≤Q≤120000000，1≤ai,bi≤N。

 

Solution

　　题面描述得至关别扭啊，小C也是看过discuss后才知道它给的是一座森林。

　　正常来讲学过树链剖分的人只要看到这一条件，题目就已经作完了。

　　题目实际上要咱们求的是，求一棵树的树链剖分方案，使得全部点到根的路径上所走的轻边数量的最大值最小。

　　学过树链剖分的人都应该知道，若是按照重链剖分(即重边连向结点更多的子树)来划分一棵树，那么从某个点出发走到根节点最多只要走logN条重链。

　　树剖复杂度的证实在此提一下，其实也很简单，设x个结点组成的树中，全部点到根路径上轻边数量最大值的最大值为f(x)。

　　那么就有，因为f(x)是递增的，因此f(x)一定是从f((x-1)/2)+1转移获得的。

　　因此转移式化简就是：，很显然f(N)是logN级别的。

　　所以答案最大不会超过logN，并且这样形态的树是一棵二叉树。

　　但实际上这道题所说的链剖不是真正意义上的链剖，划分出来的路径是能够拐弯的，所以答案会比logN更小一些，实际上在树的形态是三叉树的状况下答案达到最大，大概是log_3(N)级别。

　　而后知道了以上这些结论咱们能够很是容易地写出DP的状态和转移方程：

　　设f[i][j][0/1]表示在结点i的子树中，答案为j，且该子树可否继续向上连边的链剖方案数。转移花log^2随便转移一下就行。

　　最后判断答案的时候简单地根据该处的DP值是否为0来取最小值是不可行的，由于方案数极可能对Q取模等于0，因此这里有一个小技巧，就是方案数为Q的倍数的时候，将对Q取模的值设为Q便可。

　　时间复杂度。

#include <cstdio>
#include <algorithm>
#include <cstring>
#define ll long long
#define MN 100005
#define MS 12
using namespace std;
struct edge{int nex,to;}e[MN<<1];
int n,m,mod,pin;
int hr[MN],f[MN][2][MS];

inline int read()
{
    int n=0,f=1; char c=getchar();
    while (c<'0' || c>'9') {if(c=='-')f=-1; c=getchar();}
    while (c>='0' && c<='9') {n=n*10+c-'0'; c=getchar();}
    return n*f;
}

inline void rw(int& x,int y) {x+=y; if(x>mod)x-=mod;}
inline int modu(ll x) {return !x?0:(x-1)%mod+1;}
inline void ins(int x,int y) {e[++pin]=(edge){hr[x],y}; hr[x]=pin;}

void dfs(int x,int fat)
{
    register int i,j,k;
    int g[3][MS],h[3][MS+5];
    memset(g,0,sizeof(g));
    memset(h,0,sizeof(h));
    g[0][0]=1;
    for (i=hr[x];i;i=e[i].nex)
    {
        if (e[i].to==fat) continue;
        dfs(e[i].to,x);
        for (j=0;j<MS;++j)
            for (k=0;k<MS;++k)
            {
                rw(h[2][max(j+1,k)],modu(1LL*g[2][k]*(f[e[i].to][0][j]+f[e[i].to][1][j])));
                rw(h[2][max(j,k)],modu(1LL*g[1][k]*f[e[i].to][0][j]));
                rw(h[1][max(j+1,k)],modu(1LL*g[1][k]*(f[e[i].to][0][j]+f[e[i].to][1][j])));
                rw(h[1][max(j,k)],modu(1LL*g[0][k]*f[e[i].to][0][j]));
                rw(h[0][max(j+1,k)],modu(1LL*g[0][k]*(f[e[i].to][0][j]+f[e[i].to][1][j])));
            }
        for (j=0;j<MS;++j)
            g[0][j]=h[0][j],h[0][j]=0,
            g[1][j]=h[1][j],h[1][j]=0,
            g[2][j]=h[2][j],h[2][j]=0;
    }
    for (i=0;i<MS;++i)
        f[x][0][i]=g[0][i],rw(f[x][0][i],g[1][i]),f[x][1][i]=g[2][i];
}

int main()
{
    register int i,x,y;
    n=read(); m=read(); mod=read(); 
    if (m!=n-1) return 0*printf("-1\n-1");
    for (i=1;i<=m;++i)
        x=read(),y=read(),
        ins(x,y),ins(y,x);
    dfs(1,0);
    for (i=0;i<MS;++i)
        if (f[1][0][i]||f[1][1][i])
            return 0*printf("%d\n%d",i,(f[1][0][i]+f[1][1][i])%mod);
}

 

Last Word

　　心血来潮写了一下小D的mint，交上去居然T了，好感度--，表示这模板不是很想再用了。