BZOJ3836: [Poi2014]Tourism

因为任意两点间不存在点数超过10的简单路径，所以dfs搜索树的层数不超过10
对每个联通块建dfs搜索树，因为所有非树边都是返祖边，可以顺着dfs序dp，三进制状态压每一层的状态，0表示未被覆盖，1表示建了旅游站，2表示没建被覆盖

其实感觉这个复杂度很虚qwq

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void up(int &x,const int &y){if(x<y)x=y;}
inline void down(int &x,const int &y){if(x>y)x=y;}
const int maxn = 51000;
const int maxv = 210000;

int pw[12];
int n,m;
int c[maxn];
struct edge{int y,nex;}a[maxn<<1]; int len,fir[maxn];
inline void ins(const int x,const int y){a[++len]=(edge){y,fir[x]};fir[x]=len;}

int f[2][maxv],now;
int dfn[maxn],id,dep[maxn],mxd;
int get(const int i,const int d){return i/pw[d]%3;};
void dfs(const int x)
{
    dfn[x]=++id; up(mxd,dep[x]);
    for(int i=0;i<pw[mxd+1];i++) if(f[now][i]!=INT_MAX)
    {
        int tmp=f[now][i];
        int kk=get(i,dep[x]); int ii=i-kk*pw[dep[x]];
        down(f[!now][ii],tmp);
        int i1=-1,i2=ii+pw[dep[x]];
        for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(dfn[y])
        {
            kk=get(ii,dep[y]);
            if(kk==0) i2+=pw[dep[y]]<<1;
            if(kk==1) i1=ii;
        }
        if(i1!=-1) i1+=pw[dep[x]]<<1,down(f[!now][i1],tmp);
        down(f[!now][i2],tmp+c[x]);
    }
    for(int i=0;i<pw[mxd+1];i++) f[now][i]=INT_MAX; 
    now=!now;

    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(!dfn[y])
        dep[y]=dep[x]+1,dfs(y);
    for(int i1=0;i1<pw[dep[x]];i1++)
        for(int i2=0;i2<pw[mxd-dep[x]];i2++)
            f[now][i1+i2*pw[dep[x]+1]]=INT_MAX;
    int ok=1;
}

int re;

int main()
{
    pw[0]=1; for(int i=1;i<12;i++) pw[i]=pw[i-1]*3;

    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&c[i]);
    for(int i=1;i<=m;i++)
    {
        int x,y; scanf("%d%d",&x,&y);
        ins(x,y); ins(y,x);
    }

    for(int i=1;i<=n;i++) if(!dfn[i])
    {
        mxd=0;
        for(int k=0;k<2;k++) for(int j=0;j<pw[11];j++) f[k][j]=INT_MAX;
        now=0; f[now][0]=0;
        dfs(i);
        int tmp=INT_MAX;
        for(int j=0;j<pw[mxd+1];j++) down(tmp,f[now][j]);
        re+=tmp;
    }
    printf("%d\n",re);

    return 0;
}