AtCoder Petrozavodsk Contest 001 F - XOR Tree 状压dp

题意

给出一棵树,边有边权,每次操作可以将一条链上的所有边异或上一个相同的数。问最少多少次操作后每条边的边权都变为0。
n<=100000,边权<=15

分析

我们可以把一个点的点权设为所有与它相连的边的权值异或和。
那么所有点的点权为0和所有边的边权为0互为充要条件,这样每次操作就变成了修改两个点的点权。
状压dp一下即可。

代码

#include
#include
#include
#include
#include
#include
#define MIN(x,y) x=min(x,y)
using namespace std;

const int N=100005;
const int inf=1000000000;

int n,a[N],bin[20],cnt[70005],f[70005],w[20];
vector<int> vec[20];

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

void dp(int S)
{
    for (int i=0;i16];i++) cnt[i]=cnt[i>>1]+(i&1),vec[cnt[i]].push_back(i),f[i]=inf;
    f[S]=0;
    for (int i=16;i>0;i--)
        for (int j=0;jint s=vec[i][j];
            if (f[s]==inf) continue;
            for (int x=0;x<16;x++)
                for (int y=0;y<16;y++)
                {
                    if (!(s&bin[x])||!(s&bin[y])||x==y) continue;
                    MIN(f[s^bin[x]^bin[x^y]^bin[y]],f[s]+1+((s&bin[x^y])>0));
                }
        }
}

int main()
{
    bin[0]=1;
    for (int i=1;i<=16;i++) bin[i]=bin[i-1]*2;
    n=read();
    for (int i=1;iint x=read(),y=read(),z=read();
        a[x+1]^=z;a[y+1]^=z;
    }
    for (int i=1;i<=n;i++) w[a[i]]++;
    int S=0,ans=0;
    for (int i=1;i<=15;i++) ans+=w[i]/2,w[i]%=2,S+=bin[i]*w[i];
    dp(S);
    printf("%d",ans+f[0]);
    return 0;
}

你可能感兴趣的:(状压dp)