zjoi2015d1题解

闲来无事做了丽洁姐姐的题

t1给一棵树 每个点有点权 每次修改点权 修改后询问每个点到树的带权重心的带权距离是多少 每个点度数不超过20

很显然的一个点分树。。。

我们记一下

每个点的子树中的所有点到该点的带权距离。
每个点的子树中的所有点到该点的父亲的带权距离。
每个点的子树中的所有点的权值和。

大概就可以随便做了（当然说起来简单实际上写的又臭又长。。。）

转移的时候从分治重心开始往分治子树走看哪个更优

#include<bits/stdc++.h>
#define LL long long
const int maxn = ;
using namespace std;
int m,n;
int first[maxn],to[maxn],next[maxn],val[maxn],cnt;
inline void add(int u,int v,int w)
{
    to[++cnt]=v;
    val[cnt]=w;
    next[cnt]=first[u];
    first[u]=cnt;
}
int dep[maxn],dis[maxn],fat[maxn][],s[maxn];
int id[maxn],ip[maxn],top;
void dfs(int u,int fa)
{
    s[++top]=u;
    if(!id[u])id[u]=top;
    dep[top]=dep[ip[fa]]+;ip[u]=top;
    for(int i=first[u];i;i=next[i])
    {
        int v=to[i];
        if(v==fa)continue;
        dis[v]=dis[u]+val[i];
        dfs(v,u);s[++top]=u;dep[top]=dep[ip[fa]+];
   }
}
void make()
{
    for(int i=;i<=top;i++) fat[i][]=i;
    for(int j=;j<=;j++)
        for(int i=;i<=top;i++)
            if(i+(<<j)-<=top)
            {
                int x=fat[i][j-],y=fat[i+(<<j-)][j-];
                if(dep[fat[i][j-]]<dep[fat[i+(<<j-)][j-]]) fat[i][j]=fat[i][j-];
                else fat[i][j]=fat[i+(<<j-)][j-];
            }
}
inline int query(int l,int r)
{
    int len=r-l+,k=;
    for(k=;<<k+<=len;k++);
    if(dep[fat[l][k]]<dep[fat[r-(<<k)+][k]])return fat[l][k];
    else return fat[r-(<<k)+][k];
}
inline int lca(int u,int v)
{
    if(id[u]>id[v]) swap(u,v);
    return s[query(id[u],id[v])];
}
inline int caldis(int u,int v)
{
    int LCA=lca(u,v);
    return dis[u]+dis[v]-*dis[LCA];
}
int rt,sum,f[maxn],size[maxn],vis[maxn];
inline void GetRT(int x,int fa)
{
    size[x]=;f[x]=;
    for(int i=first[x];i;i=next[i])
    {
        if(to[i]==fa || vis[to[i]])continue;
        GetRT(to[i],x);size[x]+=size[to[i]];
        f[x]=max(f[x],size[to[i]]);
    }
    f[x]=max(f[x],sum-size[x]);
    if(f[x]<f[rt])rt=x;
}
int ret,dv[maxn],par[maxn];
LL ans[maxn],anss[maxn],summ[maxn];
inline void work(int x)
{
    vis[x]=;summ[x]=ret;
    for(int i=first[x];i;i=next[i])
    {
        if(vis[to[i]])continue;
        rt=,sum=size[to[i]];
        GetRT(to[i],);
        par[rt]=x;work(rt);
    }
}
LL cal(int u)
{
    LL ret=ans[u];
    for(int i=u;par[i];i=par[i])
    {
        LL delt=caldis(par[i],u);
        ret+=(ans[par[i]]-anss[i]);
        ret+=delt*(summ[par[i]]-summ[i]);
    }
    return ret;
}
LL update(int u,int va)
{
    summ[u]+=va;
    for(int i=u;par[i];i=par[i])
    {
        LL di=caldis(par[i],u);
        summ[par[i]]+=va;
        anss[i]+=va*di;
        ans[par[i]]+=va*di;
    }
}
int last=;
LL query(int u)
{
    LL ka=cal(u);
    for(int i=first[u];i;i=next[i])
    {
        LL tmp=cal(to[i]);
        if(tmp < ka)return query(to[i]);
    }
    last=u;
    return ka;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=;i<n;i++)
    {
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        add(u,v,w),add(v,u,w);
    }top=,dfs(,);
    //cout<<top<<endl;
    make();sum=f[]=n;GetRT(,);
    work(rt);
    for(int i=;i<=m;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        update(a,b);
        printf("%lld\n",query(last));
    }
}

t2图上每个边是0~1的随机实数 求最小生成树上最大边的期望 n<=10

clj：我们可以积分啊

进一步分析可以注意到，考虑一个x，如果<x的边合起来不能使得图联通，<=x的边合起来能够使得图联通，那么这个图的最小瓶颈生成上的最大边就是x。

那么，用WC讲过的同样的方法，我们可以得到一个多项式P(x)，表示<x的边不能使得图联通的概率。

那么注意到，我们只需要对P(x)从0到1求积分就是答案了。为什么呢？因为P(x)也是答案>x的概率，这样相当于一个分部积分。

然后积了两个多小时的分...

后来觉得可以考虑状压dp

可以点这个链接去看一下状压dp做法（白积了半天分QAQ

http://blog.csdn.net/skywalkert/article/details/47792065

#include <cstdio>
const int maxn = , maxm = ;
int n, m, e[maxn], sz[ << maxn], cnt[ << maxn];
long long c[maxm][maxm], f[ << maxn][maxm], g[ << maxn][maxm];
double ans;
int main()
{
    scanf("%d%d", &n, &m);
    for(int i = ; i < m; ++i)
    {
        int u, v;
        scanf("%d%d", &u, &v);
        --u;
        --v;
        e[u] |=  << v;
        e[v] |=  << u;
    }
    c[][] = ;
    for(int i = ; i <= m; ++i)
    {
        c[i][] = c[i][i] = ;
        for(int j = ; j < i; ++j)
            c[i][j] = c[i - ][j - ] + c[i - ][j];
    }
    for(int s = ; s <  << n; ++s)
    {
        sz[s] = sz[s >> ] + (s & );
        if(sz[s] == )
        {
            g[s][] = ;
            continue;
        }
        for(int i = ; i < n; ++i)
            if((s >> i) & )
                cnt[s] += sz[e[i] & s];
        cnt[s] >>= ;
        int lowbit = s & -s;
        for(int t = (s - ) & s; t; t = (t - ) & s)
            if(t & lowbit)
                for(int i = ; i <= cnt[t]; ++i)
                    for(int j = ; j <= cnt[s ^ t]; ++j)
                        f[s][i + j] += g[t][i] * c[cnt[s ^ t]][j];
        for(int i = ; i <= cnt[s]; ++i)
            g[s][i] = c[cnt[s]][i] - f[s][i];
    }
    for(int i = ; i <= m; ++i)
        ans += (double)f[( << n) - ][i] / c[cnt[( << n) - ]][i];
    ans /= m + ;
    printf("%.6f\n", ans);
    return ;
}

状压dp

t3题意不可描述...bzoj3926

广义后缀自动机 解法不是很好描述 可以看代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = ;
int n,c;
int first[maxn],next[maxn],to[maxn],val[maxn],cnt;
int rd[maxn];
ll ans;
inline void add(int u,int v)
{
    to[++cnt]=v;
    next[cnt]=first[u];
    first[u]=cnt;
    to[++cnt]=u;
    next[cnt]=first[v];
    first[v]=cnt;
}
struct SAM
{
    int cnt;
    int fa[maxn],mx[maxn],a[maxn][];
    SAM(){cnt=;}
    inline int extend(int p,int c)
    {
        int np=++cnt;mx[np]=mx[p]+;
        while(!a[p][c]&&p)a[p][c]=np,p=fa[p];
        if(!p)fa[np]=;
        else
        {
            int q=a[p][c];
            if(mx[p]+==mx[q])fa[np]=q;
            else
            {
                int nq=++cnt;mx[nq]=mx[p]+;
                memcpy(a[nq],a[q],sizeof(a[q]));
                fa[nq]=fa[q];
                fa[np]=fa[q]=nq;
                while(a[p][c]==q)a[p][c]=nq,p=fa[p];
            }
        }
        return np;
    }
    inline void calc(){for(int i=;i<=cnt;i++)ans+=mx[i]-mx[fa[i]];}
}sam;
inline void dfs(int u,int fa,int curpos)
{
    int nexpos = sam.extend(curpos,val[u]);
    for(int i=first[u];i;i=next[i])
        if(to[i] != fa)dfs(to[i],u,nexpos);
}
int main()
{
    scanf("%d%d",&n,&c);
    for(int i=;i<=n;i++)scanf("%d",&val[i]);
    for(int i=;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        add(u,v);rd[u]++,rd[v]++;
    }
    for(int i=;i<=n;i++)
        if(rd[i] == )dfs(i,,);
    sam.calc();
    printf("%lld",ans);
}