GXOI2019 题解

最新推荐文章于 2020-06-09 16:27:08 发布
最新推荐文章于 2020-06-09 16:27:08 发布
阅读量299 收藏
点赞数 1
分类专栏： 省选，NOI题解
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sslz_fsy/article/details/102864545
版权

[GXOI/GZOI2019]与或和
按位处理，只有全部是 1 1 1 那么与出来才是 1 1 1，于是问题转换为求全部是 1 的子矩阵个数
对每个点统计它作为左下角的矩阵个数就可以不重不漏，单调栈即可
考虑或，显然有一个 1 就是 1，那么统计全部是 0 的个数减一下即可

#include<bits/stdc++.h>
#define cs const
using namespace std;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch - '0'), ch = getchar();
	return cnt * f;
}
typedef long long ll;
cs int N = 1e3 + 5;
int a[N][N], b[N][N];
int n, h[N], sta[N];
cs int Mod = 1e9 + 7;
int mul(int a, int b){ return 1ll * a * b % Mod; }
void Add(ll &a, ll b){ a += b; if(a >= Mod) a -= Mod; }
ll calc(int typ){
	memset(h, 0, sizeof(h));
	ll sum = 0;
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= n; j++){
			if(b[i][j] ^ typ) h[j] = 0;
			else ++h[j]; 
		}
		int top = 0; ll ret = 0;
		for(int j = 1; j <= n; j++){
			while(top && h[j] < h[sta[top]]){
				Add(ret, Mod - mul(h[sta[top]] - h[j], sta[top] - sta[top - 1]));
				top--;
			} sta[++top] = j; ret += h[j]; Add(sum, ret);
		}
	} return sum;
}
int main(){
	n = read();
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n; j++)
			a[i][j] = read();
	ll ret = n * (n + 1) / 2;
	ll Ans1 = 0, Ans2 = 0;
	for(int k = 0; k < 32; k++){
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++)
				b[i][j] = a[i][j] >> k & 1;
		ll Or = (ret * ret - calc(0) + Mod) % Mod;
		ll And = calc(1);
		Add(Ans1, mul(And, 1 << k));
		Add(Ans2, mul(Or, 1 << k));
	} cout << Ans1 << " " << Ans2;
	return 0;
}

[GXOI/GZOI2019]宝牌一大堆
对于国士无双：
13 种牌暴力枚举哪一个选两个即可
对于小七对：将所有大于等于2张的牌按贡献排序取前 7 个即可
对于 3 ∗ 4 + 2 3*4 + 2 3∗4+2，显然不能直接统计，考虑 d p dp dp
类似麻将 d p dp dp 的套路，维护 i , i + 1 , i + 2 i,i+1,i+2 i,i+1,i+2 用了的个数
i + 1 , i + 2 i+1,i+2 i+1,i+2 可以只枚举到 2，因为 3 张变刻子
d p [ i ] [ x ] [ y ] [ z ] [ k ] [ 0 / 1 ] dp[i][x][y][z][k][0/1] dp[i][x][y][z][k][0/1] 表示 选到 i i i， i , i + 1 , i + 2 i,i+1,i+2 i,i+1,i+2 的出现次数，面子的个数，有无雀头的方案数
暴力枚举转移即可

#include<bits/stdc++.h>
#define cs const
using namespace std;
typedef long long ll;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch - '0'), ch = getchar();
	return cnt * f;
}
cs int N = 40;
int a[N]; // 剩余牌数 
int val[N];  // 宝牌
ll ans;
ll dp[N][5][5][3][3][2];
int T;
int c[5][5] = 
{ 
1, 0, 0, 0, 0,
1, 1, 0, 0, 0,
1, 2, 1, 0, 0, 
1, 3, 3, 1, 0, 
1, 4, 6, 4, 1
};

void Mx(ll &a, ll b){ if(a < b) a = b;}
int ps[N] = {1, 9, 10, 18, 19, 27, 28, 29, 30, 31, 32, 33, 34};
void Yolanda(){
	// 枚举哪一个牌选两张 
	for(int i = 0; i < 13; i++){
		ll sum = 1;
		for(int j = 0; j < 13; j++){
			if(j == i) sum *= c[a[ps[j]]][2] * val[ps[j]] * val[ps[j]];
			else sum *= a[ps[j]] * val[ps[j]];
		} ans = max(ans, sum * 13); 
	} 
}

int value[N]; // 七连对   
void FSY(){
	int ret = 0;
	for(int i = 1; i <= 34; i++) if(a[i] >= 2) value[++ret] = c[a[i]][2] * val[i] * val[i];
	if(ret >= 7){
		ll sum = 1;
		sort(value + 1, value + ret + 1);
		for(int i = ret; i > ret - 7; i--) sum = sum * value[i];
		Mx(ans, sum * 7);
	}
}

bool tail[N];
void DP(){
	dp[1][0][0][0][0][0] = 1;
	for(int i = 1; i <= 34; i++){
		for(int j = 0; j <= 4; j++){ // j 个面子 
			for(int k = 0; k <= 4; k++){ // i 用了多少  
				for(int l = 0; l <= 2; l++){
					for(int r = 0; r <= 2; r++){
						ll *x = dp[i][j][k][l][r];
						if(!x[0] && !x[1]) continue;
						
						if(a[i] - k >= 2){ // 雀头  
							Mx(dp[i][j][k + 2][l][r][1], x[0] / c[a[i]][k] * c[a[i]][k + 2] * val[i] * val[i]) ;
						}
						
						if(j < 4){
						if(a[i] - k >= 3){ // 刻子  
							for(int t = 0; t < 2; t++)
							Mx(dp[i][j + 1][k + 3][l][r][t], x[t] / c[a[i]][k] * c[a[i]][k + 3] * (val[i]*val[i]*val[i]));
						}
						if((!tail[i]) && a[i] - k > 0 && a[i + 1] - l > 0 && a[i + 2] - r > 0 && (l^2) && (r^2)){
							// 顺子  
							for(int t = 0; t < 2; t++){
								Mx(dp[i][j + 1][k + 1][l + 1][r + 1][t], 
								x[t] / c[a[i]][k] * c[a[i]][k+1] * val[i] 
								/ c[a[i + 1]][l] * c[a[i + 1]][l + 1] * val[i + 1] 
								/ c[a[i + 2]][r] * c[a[i + 2]][r + 1] * val[i + 2] 
								);
							}
						}
						}
						
						Mx(dp[i + 1][j][l][r][0][0], x[0]);
						Mx(dp[i + 1][j][l][r][0][1], x[1]);
						if(j == 4) Mx(ans, x[1]); 
					}
				} 
			} 
		}
	}
} 
void Solve(){
	tail[8] = tail[9] = tail[17] = tail[18] = tail[26] = tail[27] = 
	tail[28] = tail[29] = tail[30] = tail[31] = tail[32] = tail[33] = tail[34] = 1;
	while(1){
		char c[3]; scanf("%s", c);
		if(c[0] == '0') break;
		if(c[1] == 'm') a[c[0] - '0']--;
		if(c[1] == 's') a[c[0] - '0' + 9]--;
		if(c[1] == 'p') a[c[0] - '0' + 18]--;
		if(c[0] == 'E') a[28]--;
		if(c[0] == 'S') a[29]--;
		if(c[0] == 'W') a[30]--;
		if(c[0] == 'N') a[31]--;
		if(c[0] == 'Z') a[32]--;
		if(c[0] == 'B') a[33]--;
		if(c[0] == 'F') a[34]--; 
	}
	while(1){
		char c[3]; scanf("%s", c);
		if(c[0] == '0') break;
		if(c[1] == 'm') val[c[0] - '0'] = 2;
		if(c[1] == 's') val[c[0] - '0' + 9] = 2;
		if(c[1] == 'p') val[c[0] - '0' + 18] = 2;
		if(c[0] == 'E') val[28] = 2;
		if(c[0] == 'S') val[29] = 2;
		if(c[0] == 'W') val[30] = 2;
		if(c[0] == 'N') val[31] = 2;
		if(c[0] == 'Z') val[32] = 2;
		if(c[0] == 'B') val[33] = 2;
		if(c[0] == 'F') val[34] = 2; 
	}
	Yolanda(); // 国士无双  
	FSY(); // 七连对  
	DP(); // 3 * 4 + 2
	cout << ans << '\n';
}
void Init(){
	ans = 0;
	for(int i = 1; i <= 34; i++) a[i] = 4, val[i] = 1;
	memset(dp, 0, sizeof(dp));
} 
int main(){
	T = read();
	while(T--) Init(), Solve();
	return 0;
}

[GXOI/GZOI2019]特技飞行
组合题，由于保证交点数 ≤ 5 e 5 \le 5e5 ≤5e5 ，暴力求出，离散化，坐标选择，扫描线
考虑第一个问，显然最大最小取到的时候 a , b a,b a,b 的个数取到端点
最大的时候，发现可以全部交换，次数为逆序对个数
最小的时候，发现每一个置换必须换 l e n − 1 len-1 len−1 次，个数是 n − 置 换 数 n-置换数 n−置换数

#include<bits/stdc++.h>
#define cs const
using namespace std;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch - '0'), ch = getchar();
	return cnt * f;
}
typedef long long ll;
typedef long double ld;
cs int N = 1e5 + 5, M = 5e5 + 5, K = 7e5 + 5;
int n, k; ll a, b, c;
int Xs, Xe;
int y[N][2], tp[N];
bool cmp(int a, int b){ return y[a][1] < y[b][1];}
struct P{
	ld x, y; 
	P(ld _x = 0, ld _y = 0){ x = _x, y = _y; }
	P operator + (cs P &a){ return P(x+a.x, y+a.y); }
	P operator - (cs P &a){ return P(x-a.x, y-a.y); }
	P operator * (ld k){ return P(x * k, y * k);}
}p[M], q[N]; int ret;
struct Line{ P a, b; }L[N];
P Getpos(Line u, Line v){
	ld d1 = u.a.y - v.a.y, d2 = u.b.y - v.b.y;
	if(d1 < 0) d1 = -d1; if(d2 < 0) d2 = -d2;
	ld k = 1.0 * d1 / (d1 + d2);
	return u.a + ((u.b - u.a) * k);
} 
struct node{ 
	int id, h; 
	node(int i = 0, int _h = 0){ id = i; h = _h; }
	bool operator < (cs node &a) cs{ return h < a.h;}
};
multiset<node> S;
typedef multiset<node>::iterator Int;
ll ans1, ans2;
int fa[N];
int find(int x){ return x == fa[x] ? x : fa[x] = find(fa[x]);}
int r[N];
ld Dx[K], Dy[K]; int sx, sy;
struct que{ int l, r, op; };
vector<que> v[K]; 
vector<int> u[K]; 
#define pb push_back
struct BIT{
	int c[K]; 
	void add(int x,int v){ for(;x<=sy;x+=x&-x) c[x] += v; }
	int ask(int x){ int ans = 0; for(;x;x-=x&-x) ans += c[x]; return ans; }
}bit;
int Solve(){
	ll sum = 0;
	for(int i = 1; i <= ret; i++){
		ld x = p[i].x, y = p[i].y;
		p[i].x = x + y; p[i].y = x - y;
		Dx[++sx] = p[i].x;
		Dy[++sy] = p[i].y;
	}
	for(int i = 1; i <= k; i++){
		Dx[++sx] = q[i].x - r[i];
		Dx[++sx] = q[i].x + r[i];
		Dy[++sy] = q[i].y - r[i];
		Dy[++sy] = q[i].y + r[i];
	} 
	sort(Dx + 1, Dx + sx + 1);
	sort(Dy + 1, Dy + sy + 1); 
	sx = unique(Dx + 1, Dx + sx + 1) - (Dx + 1);
	sy = unique(Dy + 1, Dy + sy + 1) - (Dy + 1);
	for(int i = 1; i <= ret; i++){
		int x = lower_bound(Dx + 1, Dx + sx + 1, p[i].x) - Dx;
		int y = lower_bound(Dy + 1, Dy + sy + 1, p[i].y) - Dy;
		u[x].pb(y); 
	}
	for(int i = 1; i <= k; i++){
		int lx = lower_bound(Dx + 1, Dx + sx + 1, q[i].x - r[i]) - Dx;
		int rx = lower_bound(Dx + 1, Dx + sx + 1, q[i].x + r[i]) - Dx;
		int ly = lower_bound(Dy + 1, Dy + sy + 1, q[i].y - r[i]) - Dy;
		int ry = lower_bound(Dy + 1, Dy + sy + 1, q[i].y + r[i]) - Dy;
		v[lx].pb((que){ly, ry, 1});
		v[rx + 1].pb((que){ly, ry, -1});
	}
	for(int i = 1; i <= sx; i++){
		for(int j = 0; j < v[i].size(); j++){
			que nx = v[i][j];
			bit.add(nx.l, nx.op);
			bit.add(nx.r + 1, - nx.op);
		}
		for(int j = 0; j < u[i].size(); j++) sum += bit.ask(u[i][j]) > 0;
	} return sum;
}

int main(){
	n = read(), a = read(), b = read(), c = read(), Xs = read(), Xe = read();
	for(int i = 1; i <= n; i++) fa[i] = i;
	for(int i = 1; i <= n; i++) y[i][0] = read();
	for(int i = 1; i <= n; i++) y[i][1] = read(), tp[i] = i;
	sort(tp + 1, tp + n + 1, cmp);
	for(int i = 1; i <= n; i++){
		int fx = find(i), fy = find(tp[i]);
		if(fx ^ fy) fa[fx] = fy;
	}
	for(int i = 1; i <= n; i++) L[i] = (Line){P(Xs, y[i][0]), P(Xe, y[i][1])};
	for(int i = 1; i <= n; i++){
		Int it = S.lower_bound(node(i, y[i][1]));
		for(;it != S.end(); it++){
			p[++ret] = Getpos(L[it->id], L[i]);
		} S.insert(node(i, y[i][1]));
	}
	ans1 = ret * a;
	int cnt = 0;
	for(int i = 1; i <= n; i++) cnt += (find(i) == i);
	ans2 = (n - cnt) * a + (ret - (n - cnt)) * b;
	k = read();
	for(int i = 1; i <= k; i++){
		int x = read(), y = read(); r[i] = read();
		q[i] = P(x + y, x - y);
	} ll sum = 1ll * Solve() * c; 
	if(ans1 > ans2) swap(ans1, ans2);
	cout << ans1 + sum << " " << ans2 + sum;
	return 0;
}

[GXOI/GZOI2019]逼死强迫症
有一种做法是暴力枚举两个1块的起始和结束位置，显然中间只有 2 种可能，两边是斐波那契拼起来的
考虑递推，令 f i f_i fi​ 为斐波那契数列， F i F_i Fi​ 为长度为 i i i 的答案
首先当前点可以由前两个，放两个横或者放一个竖转移过来
F i = F i − 1 + F i − 2 F_i=F_{i-1}+F_{i-2} Fi​=Fi−1​+Fi−2​
考虑当前点作为一个1的右端点的新增贡献
枚举左端点
F i = 2 ∗ ∑ j = 1 i − 3 f j F_i=2*\sum_{j=1}^{i-3}f_j Fi​=2∗∑j=1i−3​fj​
令 s i = ∑ j = 1 i f j s_i = \sum_{j=1}^{i}f_j si​=∑j=1i​fj​
那么有 F i = F i − 1 + F i − 2 + s i − 3 ∗ 2 = F i − 1 + F i − 2 + 2 ∗ f i − 1 − 2 F_i=F_{i-1}+F_{i-2}+s_{i-3}*2=F_{i-1}+F_{i-2}+2*f_{i-1}-2 Fi​=Fi−1​+Fi−2​+si−3​∗2=Fi−1​+Fi−2​+2∗fi−1​−2
矩阵乘法即可

#include<bits/stdc++.h>
#define cs const
using namespace std;
cs int Mod = 1000000007;
int add(int a, int b){ return a + b >= Mod ? a + b - Mod : a + b; }
int mul(int a, int b){ return 1ll * a * b % Mod; }
int T, n;
cs int N = 5;
struct matrix{
	int a[N][N];
	matrix(){ memset(a, 0, sizeof(a)); }
	matrix operator * (cs matrix &A){
		matrix B; for(int i = 0; i < N; i++)
			for(int j = 0; j < N; j++) for(int k = 0; k < N; k++)
				B.a[i][j] = add(B.a[i][j], mul(a[i][k], A.a[k][j]));
		return B;
	}
}I, TR, ST;
matrix power(matrix A, int b){
	matrix ans = I; for(;b;b>>=1, A=A*A) if(b&1) ans = ans*A; return ans;
}
int main(){
	scanf("%d", &T);
	TR.a[0][0] = 1; TR.a[0][1] = 1;
	TR.a[1][0] = 1;
	TR.a[2][0] = 2; TR.a[2][2] = 1; TR.a[2][3] = 1;
	TR.a[3][2] = 1;
	TR.a[4][0] = Mod-2; TR.a[4][4] = 1;
	ST.a[0][0] = 2; ST.a[0][1] = 0; ST.a[0][2] = 3; ST.a[0][3] = 2; ST.a[0][4] = 1;
	for(int i = 0; i < N; i++) I.a[i][i] = 1;
	while(T--){
		scanf("%d", &n);
		if(n <= 2){ puts("0"); continue; }
		matrix A = TR, B = ST;
		B = B * power(A, n - 3); 
		cout << B.a[0][0] << '\n';
	} return 0;
}

[GXOI/GZOI2019]旅行者
D i j s k t r a Dijsktra Dijsktra 染色即可，枚举连接两个染色块的边更新答案

#include<bits/stdc++.h>
#define N 100050
#define M 500050
#define LL long long
#define inf 10000000000000000
using namespace std;
int read(){
	int cnt = 0; char ch = 0;
	while(!isdigit(ch)) ch = getchar();
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt;
}
int first[N], nxt[M], to[M], w[M], tot;
void add(int x, int y, int z){
	nxt[++tot] = first[x], first[x] = tot;
	to[tot] = y, w[tot] = z;
}
int T, n, m, k, a[N];
struct Node{ int x, y, z;} E[M];
LL dis[N][2]; int col[N][2];
void Dijsktra(int id){
	priority_queue<pair<int, int> > q;
	for(int i=1; i<=n; i++) dis[i][id] = inf, col[i][id] = 0;
	for(int i=1; i<=k; i++){
		col[a[i]][id] = a[i], dis[a[i]][id] = 0;
		q.push(make_pair(-dis[a[i]][id], a[i])); 
	} 
	while(!q.empty()){
		int x = q.top().second; q.pop();
		for(int i=first[x];i;i=nxt[i]){
			int t = to[i]; 
			if(dis[t][id] > dis[x][id] + w[i]){
				dis[t][id] = dis[x][id] + w[i];
				col[t][id] = col[x][id];
				q.push(make_pair(-dis[t][id], t));
			}
		}
	}
	
}
int main(){
	T = read();
	while(T--){
		n = read(), m = read(), k = read();
		for(int i=1; i<=n; i++) first[i] = 0; tot = 0;
		for(int i=1; i<=m; i++){
			int x = read(), y = read(), z = read();
			if(x ^ y) add(x, y, z); E[i] = (Node){x, y, z};
		}
		for(int i=1; i<=k; i++) a[i] = read();
		Dijsktra(0);
		for(int i=1; i<=n; i++) first[i] = 0; tot = 0;
		for(int i=1; i<=m; i++){
			if(E[i].y ^ E[i].x) add(E[i].y, E[i].x, E[i].z);
		}
		Dijsktra(1);
		LL ans = inf;
		for(int i=1; i<=m; i++){
			int x = E[i].x, y = E[i].y, z = E[i].z;
			if(col[x][0] && col[y][1] && col[x][0] != col[y][1]) 
				ans = min(ans, dis[x][0] + dis[y][1] + (LL)z);
		} printf("%lld\n",ans);
	} return 0;
}

[GXOI/GZOI2019]旧词
更 L C A LCA LCA 那道题一样
考虑到 d e p [ l c a ( x , y ) ] k = ( d e p [ l c a ( x , y ) ] k − ( d e p [ l c a ( x , y ) ] − 1 ) k ) + ( ( d e p [ l c a ( x , y ) ] − 1 ) k − ( d e p [ l c a ( x , y ) ] − 2 ) k ) + . . . dep[lca(x,y)]^k=(dep[lca(x,y)]^k-(dep[lca(x,y)]-1)^k)+((dep[lca(x,y)]-1)^k-(dep[lca(x,y)]-2)^k)+... dep[lca(x,y)]k=(dep[lca(x,y)]k−(dep[lca(x,y)]−1)k)+((dep[lca(x,y)]−1)k−(dep[lca(x,y)]−2)k)+...
而一个点恒定加的一个常数 a x a_x ax​，树剖，线段树，预处理 s u m ( a x ) sum(a_x) sum(ax​) 就可以 p u s h d o w n pushdown pushdown 了

#include<bits/stdc++.h>
#define cs const
using namespace std;
cs int Mod = 998244353;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch - '0'), ch = getchar();
	return cnt * f;
}
cs int N = 5e4 + 5;
int n, m, k;
int add(int a, int b){ return a + b >= Mod ? a + b - Mod : a + b; }
int mul(int a, int b){ return 1ll * a * b % Mod; }
void Add(int &a, int b){ a = add(a, b); }
int ksm(int a, int b){ int ans = 1; for(;b;b>>=1, a=mul(a,a)) if(b&1) ans=mul(ans,a); return ans; }
vector<int> v[N];
int fa[N], dep[N], top[N], siz[N], son[N], in[N], sign;
int a[N], pw[N];
typedef pair<int, int> pi;
#define mp make_pair
int ans[N];
vector<pi> q[N];

namespace seg{
	int sum[N << 2], val[N << 2], tg[N << 2];
	#define mid ((l+r)>>1)
	void ins(int x, int l, int r, int p, int v){
		Add(sum[x],v); if(l == r) return; 
		if(p <= mid) ins(x<<1, l, mid, p, v); else ins(x<<1|1, mid+1, r, p, v);
	}
	void pushup(int x){ val[x] = add(val[x<<1], val[x<<1|1]);}
	void pushnow(int x, int v){ Add(val[x], mul(v, sum[x])); Add(tg[x], v); }
	void pushdown(int x){ if(tg[x]) pushnow(x<<1, tg[x]), pushnow(x<<1|1, tg[x]), tg[x] = 0; }
	void modify(int x, int l, int r, int L, int R, int v){
		if(L<=l && r<=R){ pushnow(x, v); return; }
		pushdown(x);
		if(L<=mid) modify(x<<1, l, mid, L, R, v);
		if(R>mid) modify(x<<1|1, mid+1, r, L, R, v);
		pushup(x);
	}
	int query(int x, int l, int r, int L, int R){
		if(L<=l && r<=R) return val[x]; pushdown(x); int ans = 0;
		if(L<=mid) Add(ans, query(x<<1, l, mid, L, R));
		if(R>mid) Add(ans, query(x<<1|1, mid+1, r, L, R));
		return ans;
	}
}
void dfs(int u){
	siz[u] = 1; 
	for(int i = 0; i < v[u].size(); i++){
		int t = v[u][i]; 
		dep[t] = dep[u] + 1; fa[t] = u;
		a[t] = add(pw[dep[t]], Mod - pw[dep[u]]);
		dfs(t); siz[u] += siz[t];
		if(siz[son[u]] > siz[t]) son[u] = t;
	}
}
void dfs2(int u, int tp){
	top[u] = tp; in[u] = ++sign; seg::ins(1,1,n,in[u],a[u]);
	if(son[u]) dfs2(son[u], tp);
	for(int i = 0; i < v[u].size(); i++) if(v[u][i]^son[u]) dfs2(v[u][i], v[u][i]);
}
void modify(int x){ while(x){ seg::modify(1,1,n,in[top[x]],in[x],1); x=fa[top[x]]; } }
int ask(int x){ int ans=0; while(x){ Add(ans,seg::query(1,1,n,in[top[x]],in[x])); x=fa[top[x]];} return ans; }
int main(){
	n = read(), m = read(), k = read();
	for(int i = 1; i <= n; i++) pw[i] = ksm(i, k);
	for(int i = 2; i <= n; i++){
		int x = read(); v[x].push_back(i);
	} dep[1] = a[1] = 1; dfs(1); dfs2(1, 1);
	for(int i = 1; i <= m; i++){
		int x = read(), y = read();
		q[x].push_back(mp(y, i));
	}
	for(int i = 1; i <= n; i++){
		modify(i); 
		for(int j = 0; j < q[i].size(); j++){
			int x = q[i][j].first; 
			ans[q[i][j].second] = ask(x);
		}
	} 
	for(int i = 1; i <= m; i++) cout << ans[i] << '\n';
	return 0;
}
FSYo
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
2019正睿Day1题解.pdf
10-15
正睿2019Day1题解 T1序列seq T2灯泡bulb T3比赛match
蓝桥杯-蓝桥杯2019VIP题目+题解.zip
最新发布
03-04
蓝桥杯 蓝桥杯_蓝桥杯2019VIP题目+题解
GZOI/GXOI2019
weixin_30667301的博客
05-03 143
陆陆续续做完了…… 与或和（单调栈） 这是一道一眼题…… 看到位运算，按位考虑贡献。对于每一位，将矩阵中的元素变为“当前元素的这一位是否为\(1\)”，那么原矩阵变为\(01\)矩阵。在\(01\)矩阵中能够对\(AND\)产生贡献的是全\(1\)的矩阵，能够对\(OR\)产生贡献的是存在\(1\)的矩阵，那么我们需要求全\(1\)和全\(0\)矩阵，单调栈做一下就可以了。 代码 宝牌一大堆（DP...
[GXOI/GZOI2019]宝牌一大堆(dp)
weixin_33762321的博客
06-18 183
luogu bzoj 这个麻将题还算挺友善的，比隔壁zjoi的要好得多。。。 比较正常的做法是五维dp 但事实上六维dp也是完全不会被卡的 七对子选权值最高的七个，国士无双直接$13^2$暴力 $dp[i][j][0/1][k][l][m]$表示枚举到了第i张牌，已经凑了j个面子，有无雀头，第i张牌已经用了k张，第i+1张牌用了l张，第i+2张牌用了m张，直接暴力转移。。。 然...
GXOI/GZOI2019题解
weixin_30394633的博客
04-25 243
GXOI/GZOI2019题解 P5300 [GXOI/GZOI2019]与或和 一眼题。。 显然枚举每个二进制位，答案就变成了全1子矩阵数量。 这个xjb推推，单调栈一下就行了。 #include<bits/stdc++.h> #define il inline #define vd void #define mod 1000000007 typedef long long ll;...
[GXOI/GZOI2019]与或和
psc233的博客
04-22 180
[GXOI/GZOI2019]与或和 解析： 首先我们分位考虑，对于二进制第i位，AND所产生的贡献为2^i*该位所构成的01矩阵中的全为1的子矩阵的个数,OR所产生的贡献为2^i*该位所构成的01矩阵中有1的子矩阵的个数 对于上述化简后的问题，我们只需要单调栈就行了 #include<bits/stdc++.h> #define ll long long using nam...
「GXOI / GZOI2019」与或和
这里RevolIA，_(:з」∠)_
05-17 326
问题 C: 与或和 时间限制:4 Sec内存限制:512 MB 提交:278解决:78 [提交] [状态] [命题人:admin] 题目描述 Freda 学习了位运算和矩阵以后，决定对这种简洁而优美的运算，以及蕴含深邃空间的结构进行更加深入的研究。 对于一个由非负整数构成的矩阵，她定义矩阵的AND值为矩阵中所有数二进制AND（&）的运算结果；定义矩阵的OR值为矩阵中所...
2019正睿Day2题解.pdf
10-15
2019正睿Day2题解 T1石子stone T2内存memory T3子集subset
ZJOI2019题面以及简要题解
03-27
这是ZJOI的题目以及简要题解。参加不了的OIer可以在家里做做看哦！
PAT甲级所有题解代码
02-03
这里包含了179道PAT甲级的代码AC题解，供有需要通过PAT考试的同学下载学习使用。
GX Simulator 7.16（三菱PLC仿真器可仿真Fx-3U CPU）
09-03
GX Simulator 7.16可在三菱PLC编程软件GX-Developer中仿真Fx-3U系列的CPU，而其它低版本的仿真器无法仿真Fx-3U类型的CPU。
【GXOI/GZOI2019】宝牌一大堆（麻将DP）（贪心）
zxyoi_dreamer的博客（不定期诈尸）
09-06 343
传送门 其实还有第三种麻将——四川麻将，规则和其他麻将都不太像，甚至在四川不同地区的规则也不一样。 而关于麻将胡牌的题，能够利用DP解决的，我们可以将这种技巧成为麻将DP 题解： 这篇题解里面用的是川麻的语言习惯。 首先注意到杠子是废的，我们永远不可能考虑用一个杠子代替一个刻子，因为它权值永远小于刻子。（没注意到没关系，反正写到DP里面杠子的转移就一行） 国士无双可以直接枚举重复的那个是什么O(...
[GXOI/GZOI2019]与或和 题解
GKxx
04-28 272
传送门 题意：给出一个n×nn\times nn×n的矩阵，求所有子矩阵的and\text{and}and值之和和or\text{or}or值之和。 显然可以把每一位分开来求，那么对于某一位而言矩阵中的元素不是000就是111。一个仅由000和111构成的子矩阵的and\text{and}and值是111当且仅当这个子矩阵中全是111，否则是000；类似地，or\text{or}or值是000当且...
「GXOI / GZOI2019」简要题解
04-16 218
「GXOI / GZOI2019」简要题解 LOJ#3083. 「GXOI / GZOI2019」与或和 https://loj.ac/problem/3083 题意:求一个矩阵的所有子矩阵的与和 和 或和。 分析: 或和与是一个东西，只要把所有数都异或上\((1<<31)-1\)然后再从总答案中减掉就能互相转化，考虑求与。 枚举每一位，转化成算有多少个全\(1\)子矩形，...
【矩阵快速幂】LOJ3086 「GXOI / GZOI2019」逼死强迫症
MEET YOU FIRST TIME
04-22 378
【题目】 LOJ 有一条2×n2\times n2×n的路，要用n−1n-1n−1个1×21\times 21×2的相同砖块和两个1×11\times 11×1的相同砖块来铺路，要求两个小砖块没有邻边，求方案数。有TTT组数据。 n≤2×109,T≤500n\leq 2\times 10^9,T\leq 500n≤2×109,T≤500 【解题思路】 首先如果没有那两个1×11\times 11×...
【GXOI/GZOI2019】旅行者
weixin_30362233的博客
07-31 117
题面 https://www.luogu.org/problem/P5304 题解 可能是“最短路”标签下最难的一道题了吧，还是要补说几句的。 对于任意一对点，他们的最短路有两种情况： 一条边直达 经过别的点中转 第一种情况直接特判掉。 对于第二种情况 跑多源$dijkstra$，正跑一次，再反跑一次。每次记录这个点的最短路是由哪个点过来的。 本来是这样想的...
[GXOI/GZOI2019]旧词 题解
GKxx
04-27 170
传送门 题意：给一棵树，qqq次询问，每次询问给出x,yx,yx,y，求∑i=1xdepk(lca(i,y))\sum\limits_{i=1}^x\mathrm{dep}^k(\mathrm{lca}(i,y))i=1∑x​depk(lca(i,y)) 当k=1k=1k=1的时候它几乎就是[LNOI2014]LCA。 k=1k=1k=1的做法： lca(i,y)\mathrm{lca}(i,y)...
「GXOI / GZOI 2019」逼死强迫症
wly_sh的博客
06-09 229
题目传送门luogu 题目传送门loj 题目描述 ITX351 要铺一条 2×N2 \times N2×N的路，为此他购买了 N块 2×12 \times 12×1的方砖。可是其中一块砖在运送的过程中从中间裂开了，变成了两块 1×11 \times 11×1的砖块！ ITX351 由此产生了一个邪恶的想法：他想要在这条路上故意把两块 1×11 \times 11×1的砖块分开铺，不让两块砖有相邻的边，其他砖块可以随意铺，直到整条路铺满。这样一定可以逼死自身强迫症 sea5！ 也许下面的剧情你已经猜到了——他
[GXOI/GZOI2019]逼死强迫症 题解
GKxx
04-28 200
传送门 题意：用n−1n-1n−1块大小为1×21\times21×2的砖和两块大小为1×11\times11×1的砖铺满2×n2\times n2×n的路，要求两块小砖不能有边相邻，求方案数。 这题简直是弟中弟。 首先设fif_ifi​表示用1×21\times21×2的砖铺满2×n2\times n2×n的路的方案数，显然fif_ifi​是斐波那契数列。 并注意到斐波那契数列的一个重要性质：∑...
2019北京市赛小学组t2题解
10-12
2019北京市赛小学组t2题解 该题目给出了一个长为n的字符串str，要求使用以下两种操作将字符串转换成相同的字符： 1. 将字符串str的第i个字符改为其相邻字符的ASCII码值的平均数 2. 将字符串str的第i个字符改为其...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值