【bzoj4568】【Scoi2016】幸运数字 （线性基+树上倍增）

Description

A 国共有 n 座城市，这些城市由 n-1 条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。一些旅行者希望游览 A 国。旅行者计划乘飞机降落在 x 号城市，沿着 x 号城市到 y 号城市之间那条唯一的路径游览，最终从 y 城市起飞离开 A 国。在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。然而，幸运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如，游览者拍了 3 张照片，幸运值分别是 5，7，11，那么最终保留在自己身上的幸运值就是 9（5 xor 7 xor 11）。有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择 5 和 11 ，可以保留的幸运值为 14 。现在，一些游览者找到了聪明的你，希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。

Input

第一行包含 2 个正整数 n ，q，分别表示城市的数量和旅行者数量。第二行包含 n 个非负整数，其中第 i 个整数 Gi 表示 i 号城市的幸运值。随后 n-1 行，每行包含两个正整数 x ，y，表示 x 号城市和 y 号城市之间有一条道路相连。随后 q 行，每行包含两个正整数 x ，y，表示这名旅行者的旅行计划是从 x 号城市到 y 号城市。N<=20000,Q<=200000,Gi<=2^60

Output

输出需要包含 q 行，每行包含 1 个非负整数，表示这名旅行者可以保留的最大幸运值。

 

题解：

这首先一看，线性基啊！不过这是树上的，我们可以用树上的数据结构维护线性基。

于是就想到了倍增。

v[i][j] 表示从 i （包括）到其第 2^j 个祖先（不包括）的线性基，每次查询往上跳，每次将线性基合并。

线性基的合并直接暴力，因为线性基只有60的长度，暴力合并不会耗太多时间。

剩下的就是裸的线性基了！！

CODE：

#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdio>
using namespace std;

int n,m,u[250005],v[250005],deg[505];
double a[505][505],g[250005],res;

void gauss(){
    for(int i=1,maxn=i;i<n;maxn=++i){
        for(int j=i+1;j<=n;j++)
            if(fabs(a[j][i])>fabs(a[maxn][i]))maxn=j;
        for(int j=1;j<=n+1;j++)swap(a[i][j],a[maxn][j]);
        for(int j=i+1;j<=n;j++){
            if(fabs(a[j][i])<1e-10)continue;
            double s=a[j][i]/a[i][i];
            for(int k=1;k<=n+1;k++)a[j][k]-=a[i][k]*s;
        }
    }
    for(int i=n;i>=1;i--){
        for(int j=i+1;j<=n;j++)
            a[i][n+1]-=a[i][j]*a[j][n+1];
        a[i][n+1]/=a[i][i];
    }
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        scanf("%d%d",u+i,v+i);
        deg[u[i]]++,deg[v[i]]++;
    }
    a[1][n+1]=-1,a[n][n]=1;
    for(int i=1;i<=m;i++){
        if(u[i]^n)a[u[i]][v[i]]=1.0/deg[v[i]];
        if(v[i]^n)a[v[i]][u[i]]=1.0/deg[u[i]];
    }
    for(int i=1;i<n;i++)a[i][i]=-1;
    gauss();
    for(int i=1;i<=m;i++)
        g[i]=a[u[i]][n+1]/deg[u[i]]+a[v[i]][n+1]/deg[v[i]];
    sort(g+1,g+m+1);
    for(int i=1;i<=m;i++)res+=g[i]*(m-i+1);
    printf("%.3f\n",res);
}