BZOJ 4712 洪水 动态dp(LCT+矩阵乘法）

30min写完，细节调了几乎两天~ 

这玩意太tm恶心了~    

#include 
#define N 200005
#define ll long long 
#define inf 10000000005 
#define lson p[x].ch[0]
#define rson p[x].ch[1]  
#define setIO(s) freopen(s".in","r",stdin)  , freopen(s".out","w",stdout)  
using namespace std;
int n,edges;
ll f[N],g[N],val[N], son[N], sta[N];     
int hd[N],to[N<<1],nex[N<<1];  
void add(int u,int v)
{
    nex[++edges]=hd[u],hd[u]=edges,to[edges]=v;  
} 
struct Node
{
    int ch[2],f,rev;
}p[N];    
struct Matrix
{
    ll a[2][2];
    ll*operator[](int x) { return a[x]; }
}t[N],tmp[N];  
Matrix operator*(Matrix a,Matrix b)
{
    Matrix c; 
    c[0][0]=min(a[0][0]+b[0][0],a[0][1]+b[1][0]); 
    c[0][1]=min(a[0][0]+b[0][1],a[0][1]+b[1][1]);   
    c[1][0]=min(a[1][0]+b[0][0],a[1][1]+b[1][0]);       
    c[1][1]=min(a[1][0]+b[0][1],a[1][1]+b[1][1]);  
    return c; 
} 
int get(int x)
{
    return p[p[x].f].ch[1]==x;         
}
int isrt(int x)
{
    return !(p[p[x].f].ch[0]==x||p[p[x].f].ch[1]==x); 
}
void pushup(int x)
{ 
    t[x]=tmp[x]; 
    if(lson) t[x]=t[lson]*t[x];     
    if(rson) t[x]=t[x]*t[rson];                
}  
void rotate(int x)
{
    int old=p[x].f,fold=p[old].f,which=get(x);  
    if(!isrt(old)) p[fold].ch[p[fold].ch[1]==old]=x; 
    p[old].ch[which]=p[x].ch[which^1],p[p[old].ch[which]].f=old; 
    p[x].ch[which^1]=old,p[old].f=x,p[x].f=fold; 
    pushup(old),pushup(x);  
}
void splay(int x)
{
    int v=0,u=x,fa; 
    for(u;!isrt(u);u=p[u].f);              
    for(u=p[u].f;(fa=p[x].f)!=u;rotate(x))
        if(p[fa].f!=u)
            rotate(get(fa)==get(x)?fa:x); 
}
void Access(int x)
{
    for(int y=0;x;y=x,x=p[x].f)
    {
        splay(x);         
        if(rson)      
        {
            tmp[x][0][0]+=min(t[rson][0][0],t[rson][0][1]);                 
        }
        if(y)
        { 
            tmp[x][0][0]-=min(t[y][0][0],t[y][0][1]);                                     
        }
        rson=y; 
        pushup(x);                   
    }
}   
void dfs(int x,int ff)
{
    p[x].f=ff;
    ll sum=0;                            
    for(int i=hd[x];i;i=nex[i])
    {
        int v=to[i];
        if(v!=ff) dfs(v,x), sum+=f[v], ++son[x];        
    } 
    if(sum!=0) f[x]=min(sum, val[x]);   
    else f[x]=val[x];                                  
    tmp[x][0][0]=(sum==0?inf:sum);      
    tmp[x][0][1]=val[x];  
    tmp[x][1][0]=tmp[x][1][1]=0;                                                          
    pushup(x);   
}        
int main()
{
    int i,j,k,m;
    // setIO("input");
    scanf("%d",&n);   
    for(i=1;i<=n;++i) scanf("%lld",&val[i]); 
    for(i=1;i