JZOJ5442. 【NOIP2017提高A组冲刺11.1】荒诞

Description

我的灵魂与我之间的距离如此遥远，而我的存在却如此真实。
——加缪《局外人》
我醒来的时候，发现满天星斗照在我的脸上。田野上的声音一直传到我的耳畔。夜的气味，土地的气味，海盐的气味，使我的两鬓感到清凉。这沉睡的夏夜的奇妙安静，像潮水一般浸透我的全身。这时，长夜将尽，汽笛叫了起来。它宣告有些人踏上旅途，要去一个从此和我无关痛痒的世界。
这时我在想一个问题：我有一个n个点，m条边的无向图，第i个点建立一个旅游站点的费用是c_i。特别地，这张图中的任意两点间不存在节点数超过10的简单路径。
为了把一切都做得完善，为了使我感到不那么孤独，我想要建造一些旅游站点使得每个点要么建立了旅游站点，要么与它有边直接相连的点里至少有一个点建立了旅游站点。我还希望这个建造方案总花费尽量少。
请求出这个花费。

Input

文件第一行有两个正整数n和m。
第二行包含n个整数，其中第i个数为c_i，表示在第i个点建立旅游站点的费用。
接下来m行，每行两个正整数u,v，表示一条边(u,v)，保证没有重边。

Output

输出只有一行，表示最小的总花费。

Sample Input

6 6
3 8 5 6 2 2
1 2
2 3
1 3
3 4
4 5
4 6

Sample Output

7

Data Constraint

对于前10%的测试点，满足所有的c_i相等。
对于前30%的测试点，满足1<=n<=20,0<=m<=50。
对于另外15%的测试点，满足每个连通块都是一棵树。
对于100%的测试点，满足1<=n<=2*10^4,0<=m<=2.5*10^4,0<=c_i<=10^4。

题解

注意到题目上“这张图中的任意两点间不存在节点数超过10的简单路径。”
这说明在这个图上面建dfs树的深度是不会大于10的。

这是图，所以即使位于下面的点也会影响到它的祖先，也就是说还要记录一个祖先的状态。
设$f_{x,s}$表示当前节点x，它祖先的三进制状态，0/1/2，表示没选且周围的都没有选/没有选但周围有至少一个选了/选了

一个点首先先要继承它父亲的状态，然后将状态下传给自己的儿子，最后又从儿子的状态更新。

数组要开滚动数组。

code

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string.h>
#include <algorithm>
#define ll long long
#define N 25003
#define M 10000000 
using namespace std;
char ch;
void read(int& n)
{
    n=0;
    ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while('0'<=ch && ch<='9')n=(n<<1)+(n<<3)+ch-'0',ch=getchar();
}
void write(int x)
{
    if(x>9)write(x/10);
    putchar(x%10+48);
}
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}

int nxt[N*2],to[N*2],b[N],cost[N],ans,n,m,x,y,tot,root;
int bz[N],f[3][60000],z[13],deep[N];

void ins(int x,int y)
{
    nxt[++tot]=b[x];
    to[tot]=y;
    b[x]=tot;
}

void dfs(int x,int fa)
{
    bz[x]=root;
    int p=deep[x]%2;
    bool myfather[13];
    memset(f[p],127,sizeof(f[p]));
    memset(myfather,0,sizeof(myfather));
    for(int i=b[x];i;i=nxt[i])
        if(deep[to[i]]<deep[x] && bz[to[i]]==root)myfather[deep[to[i]]]=1;
    for(int s=0;s<z[deep[x]-1];s++)
    {
        if(f[p^1][s]>M)continue;
        int ss=s+z[deep[x]-1];
        bool pd=0;
        for(int de=1;de<deep[x];de++)
            if(myfather[de])
            {
                if(!(ss/z[de-1]%3))ss+=z[de-1]<<1;
                if(s/z[de-1]%3==1)pd=1;
            }
        f[p][ss]=min(f[p][ss],f[p^1][s]+cost[x]);
        f[p][s+(z[deep[x]-1]<<1)*pd]=f[p^1][s];
    }

    for(int i=b[x];i;i=nxt[i])
        if(bz[to[i]]==0)
        {
            deep[to[i]]=deep[x]+1;
            dfs(to[i],x);
            for(int s=0;s<z[deep[x]];s++)
                f[p][s]=min(f[p^1][s+z[deep[x]]],f[p^1][s+(z[deep[x]]<<1)]);
        }

}

int main()
{
    freopen("absurdity.in","r",stdin);
    freopen("absurdity.out","w",stdout);
    read(n);read(m);
    z[0]=1;
    for(int i=1;i<=10;i++)
        z[i]=z[i-1]*3;
    for(int i=1;i<=n;i++)
        read(cost[i]);
    for(int i=1;i<=m;i++)
        read(x),read(y),ins(x,y),ins(y,x);

    memset(bz,0,sizeof(bz));


    ans=0;
    for(int i=1;i<=n;i++)
        if(bz[i]==0)memset(f,127,sizeof(f)),deep[root=i]=1,f[0][0]=0,dfs(i,0),
            ans+=min(f[1][1],f[1][2]);

    write(ans);

    return 0;
}