bzoj4568 [Scoi2016]幸运数字（倍增，线性基，点分治）

最新推荐文章于 2021-01-09 11:14:27 发布

1035719430

最新推荐文章于 2021-01-09 11:14:27 发布

阅读量175

点赞数

分类专栏：线性基倍增

本文链接：https://blog.csdn.net/a1035719430/article/details/86671830

版权

线性基同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

倍增

5 篇文章 0 订阅

订阅专栏

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 个非负整数，表示这名旅行者可以保留的最大幸运值。
Sample Input
4 2
11 5 7 9
1 2
1 3
1 4
2 3
1 4
Sample Output
14
11
HINT
Source

线性基具有可合并性
一种做法是倍增维护线性基后直接合并
复杂度 $O(60^2nlogn+60^2qlogn)$
跑不满加上是总时限可以卡过去
正解应该是点分治+线性基？
复杂度 $O(60nlogn+60^2q)$
代码是倍增的

#include<bits/stdc++.h>
using namespace std;
#define rep(i,j,k) for(int i = j;i <= k;++i)
#define repp(i,j,k) for(int i = j;i >= k;--i)
#define rept(i,x) for(int i = linkk[x];i;i = e[i].n)
#define P pair<int,int>
#define Pil pair<int,ll>
#define Pli pair<ll,int>
#define Pll pair<ll,ll>
#define pb push_back 
#define pc putchar
#define mp make_pair
#define file(k) memset(k,0,sizeof(k))
#define ll long long
int rd()
{
    int sum = 0;char c = getchar();bool flag = true;
    while(c < '0' || c > '9') {if(c == '-') flag = false;c = getchar();}
    while(c >= '0' && c <= '9') sum = sum * 10 + c - 48,c = getchar();
    if(flag) return sum;
    else return -sum;
}
const int p = 1e9+7;
int n,k;
int a[4010],B[4010],g[4010];
int f[4010],c[4010];
inline int calc(int a,int b){return (a+b)%p;}
inline int mul(int a,int b){return 1ll*a*b%p;}
inline void mo(int &x){x = (x%p+p)%p;}
void change(int *a,int *b)
{
	rep(i,0,2*k-2) c[i] = 0;
	rep(i,0,k-1) rep(j,0,k-1) c[i+j] = calc(c[i+j],mul(a[i],b[j]));
	repp(i,2*k-2,k)
	if(c[i])
	{
		rep(j,i-k,i-1) c[j] = calc(c[j],mul(c[i],B[j-i+k+1]));
		c[i] = 0;
	}
	rep(i,0,k-1) a[i] = c[i];
}
void work(int n)
{
	for(int i = 1;i <= n;i <<= 1)
	{
		if(n&i) change(f,g);
		change(g,g);
	}
}
int main() 
{
	n = rd();k = rd();
	repp(i,k,1) B[i] = rd();rep(i,0,k-1) a[i] = rd();
	rep(i,1,k) mo(a[i]),mo(B[i]);
    f[0] = 1;g[1] = 1;
	work(n);
	int ans = 0;
	rep(i,0,k-1) ans = calc(ans,mul(f[i],a[i]));
	printf("%d\n",ans);
	return 0;
}