NOI2022联合省选 day2代码

已于 2022-08-30 01:34:25 修改
阅读量1.2k 收藏 1
点赞数 1
分类专栏： 题解 代码 省选系列 文章标签： 算法 c++ NOI 动态规划
于 2022-04-17 20:59:53 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/liuzhangfeiabc/article/details/124236331
版权
题解 同时被 3 个专栏收录
10 篇文章 0 订阅
订阅专栏
代码
3 篇文章 0 订阅
订阅专栏
省选系列
1 篇文章 0 订阅
订阅专栏

注：由于官方数据还没公布，所以代码不一定保熟
本来想在上一篇题解里放代码，但是这样文章就太长了，所以新开一篇文章放代码

d2t1：

#include<bits/stdc++.h>
#define gc getchar()
#define pc putchar
#define li long long
using namespace std;
inline li read(){
	li x = 0;
	int y = 0,c = gc;
	while(c < '0' || c > '9') y = c,c = gc;
	while(c >= '0' && c <= '9') x = x * 10 + c - '0',c = gc;
	return y == '-' ? -x : x;
}
inline void prt(li x){
	if(x >= 10) prt(x / 10);
	pc(x % 10 + '0');
}
inline void print(li x){
	if(x < 0) pc('-'),x = -x;
	prt(x);
}
inline void file(char *s){
	char c[50];
	sprintf(c,"%s.in",s);
	freopen(c,"r",stdin);
	sprintf(c,"%s.out",s);
	freopen(c,"w",stdout);
}
int p[2010],cnt,n,m,dy[2010],wei[10010];
bool notp[2010];
int sl[8200][310],ss[2010],a[310],k,zc[2010],bg[2010],tto[8210];
const int mo = 998244353;
li mi[1000010];
bool vis[2010];
int main(){
	//file("card");
	int i,j,x,l;
	for(i = 1;i < (1 << 13);++i) wei[i] = wei[i >> 1] ^ (i & 1);
	notp[1] = 1;
	for(i = 2;i <= 2000;++i){
		if(!notp[i]) p[++cnt] = i,dy[i] = cnt;
		for(j = 1;j <= cnt && i * p[j] <= 2000;++j){
			notp[i * p[j]] = 1;
			if(i % p[j] == 0) break;
		}
	}
	n = read();
	mi[0] = 1;for(i = 1;i <= n;++i) mi[i] = mi[i - 1] * 2 % mo;
	for(i = 1;i <= n;++i){
		++ss[read()];
	} 
	for(i = 14;i <= cnt;++i) for(j = 1;p[i] * j <= 2000;++j) bg[p[i] * j] = i;
	for(i = 1;i <= 2000;++i) for(j = 1;j <= 13;++j) 
		if(i % p[j] == 0) zc[i] |= (1 << j - 1);
	for(i = 0;i < (1 << 13);++i){//不能出现的因子集合 
		for(j = 1;j <= 2000;++j) if((i & zc[j]) == 0){
			sl[i][bg[j]] += ss[j];
			tto[i] += ss[j];
		}
	}
	m = read();
	li qqq = 0;
	for(i = 1;i <= m;++i){
		k = read();
		int r = 0,q = 0;
		memset(vis,0,sizeof(vis));
		for(j = 1;j <= k;++j){
			x = read();
			if(vis[x]) continue;
			vis[x] = 1;
			if(x <= 41) q |= 1 << (dy[x] - 1);
			else a[++r] = dy[x];
		}
		li ans = 0;
		for(j = q;;j = (j - 1) & q){
			int *g = sl[j];
			li tmp = 1;
			int oo = tto[j];
			for(l = 1;l <= r;++l){
				(tmp *= mi[g[a[l]]] - 1) %= mo;
				oo -= g[a[l]];
			}
			(tmp *= mi[oo]) %= mo;
			if(wei[j]) ans -= tmp;
			else ans += tmp;
			if(!j) break;
		}
		print((ans % mo + mo) % mo);pc('\n');
		
	}
	return 0;
}

d2t2：

#include<bits/stdc++.h>
#define gc getchar()
#define pc putchar
#define li long long
using namespace std;
inline li read(){
	li x = 0;
	int y = 0,c = gc;
	while(c < '0' || c > '9') y = c,c = gc;
	while(c >= '0' && c <= '9') x = x * 10 + c - '0',c = gc;
	return y == '-' ? -x : x;
}
inline void prt(li x){
	if(x >= 10) prt(x / 10);
	pc(x % 10 + '0');
}
inline void print(li x){
	if(x < 0) pc('-'),x = -x;
	prt(x);
}
inline void file(char *s){
	char c[50];
	sprintf(c,"%s.in",s);
	freopen(c,"r",stdin);
	sprintf(c,"%s.out",s);
	freopen(c,"w",stdout);
}
int n,a[500010],x,y;
char c[1000010];
vector<li> p[500010];
multiset<li> ss;
int mx;
li wk(int qs){
	int i,j;
	li qjmx = 0;int mxwz = 0;bool fg = 0;
	li ans = 0,mm,nn,cd;
	int sz,sl = ss.size();
	for(i = qs;i <= mx;++i){//找全局最大值，值为qjmx，位置为mxwz 
		sz = p[i].size();
		if(sz > 1) fg = 1;
		if(fg){
			if(p[i][sz - 1] >= qjmx){
				qjmx = p[i][sz - 1];
				mxwz = i;
			}
		}
	} 
	for(i = qs;i <= n;++i){
		//每一层用最小的把剩下的放进去，最后用最大的把最小的放进去
		int sz = i <= mx ? p[i].size() : 0;
		sl += sz;
		if(sl <= 1){sl = 0;continue;} 
		for(j = 0;j < p[i].size();++j) ss.insert(p[i][j]);
		mm = *ss.begin();//当前最小值 
		nn = *ss.rbegin();//当前最大值 
		if(i < mxwz){//还没遇到最大值 
			ans += mm * (sl - 2) + nn;//这一层的权值是mn*(sl-2)+mx
			--sl;ss.erase(ss.find(nn));//这一层留下最大值 
		} 
		else{//已经遇到了最大值 
			cd = *(++ss.rbegin());//当前次大值 
			ans += mm * (sl - 2) + cd;//这一层的权值是mn*(sl-2)+cd
			--sl;ss.erase(ss.find(cd));//这一层留下次大值 
		} 
	}
	ss.clear();
	return ans;
}
int main(){
	//file("bracket");
	int i,j;
	n = read();x = read();y = read();
	scanf("%s",c + 1);
	int nw = 0;
	for(i = 1,j = 1;i <= n;++i,++j){
		while(c[j] == ')') --nw,++j;
		p[++nw].push_back(read());
		mx = max(mx,nw);
	}
	for(i = 1;i <= mx;++i) sort(p[i].begin(),p[i].end());
	li ans = 0,mm = 1e9,nn = 0,tot = 0;
	int sl = 0;
	int qs = 0;
	for(qs = 1;qs <= mx;++qs) if(p[qs].size() > 1) break; 
	if(qs == mx + 1){
		pc('0');pc('\n');return 0;
	} 
	if(x == 1 && y == 1){//每层用最小的把别的放进去，在外面留一个最大的 
		for(i = qs;i <= n;++i){
			int sz = i <= mx ? p[i].size() : 0;
			sl += sz;
			if(sl <= 1){sl = 0;continue;} 
			for(j = 0;j < p[i].size();++j){
				tot += p[i][j];
				ss.insert(p[i][j]);
			} 
			mm = *ss.begin();//当前最小值 
			nn = *ss.rbegin();//当前最大值 
			ans += tot + mm * (sl - 2);//这一层的权值是sum+mn*(sl-2) 
			--sl;tot -= nn;ss.erase(ss.find(nn));//外面留下最大值 
		}
	}
	else if(y == 1){//直接把小的扔进去，外面留一个最大值 
		for(i = qs;i <= n;++i){
			int sz = i <= mx ? p[i].size() : 0;
			sl += sz;
			if(sl <= 1){sl = 0;continue;} 
			for(j = 0;j < p[i].size();++j){
				tot += p[i][j];
				ss.insert(p[i][j]);
			} 
			mm = *ss.begin();//当前最小值 
			nn = *ss.rbegin();//当前最大值 
			ans += tot - nn;//这一层的权值是sum-mx
			--sl;tot -= nn;ss.erase(ss.find(nn));//外面留下最大值 
		}
	}
	else if(x == 1){//两种贪心取最优 
		ans = wk(qs); //贪心1：把全局最大值放进最里层
		//贪心2：如果开局是2 1 1...1 k结构，把前面这些层的最小值放进k这一层，后续不变 
		int o;for(o = qs + 1;o <= mx;++o) if(p[o].size() != 1) break;
		if(p[qs].size() == 2){//只有全局最大出现在前半段时可能会有影响 
			tot = 0;mm = 1e9;
			for(i = qs;i < o;++i){
				for(j = 0;j < p[i].size();++j){
					tot += p[i][j];
					mm = min(mm,p[i][j]);
				}
			}
			ss.insert(mm);
			li an = wk(o);
			ans = min(ans,an + tot - mm);
		} 
	}
	print(ans);pc('\n');
	return 0;
}

d2t3：

#include<bits/stdc++.h>
#define gc getchar()
#define pc putchar
#define li long long
#define mp make_pair
#define pli pair<li,int>
#define fi first
#define se second
#define pb push_back
using namespace std;
inline li read(){
	li x = 0;
	int y = 0,c = gc;
	while(c < '0' || c > '9') y = c,c = gc;
	while(c >= '0' && c <= '9') x = x * 10 + c - '0',c = gc;
	return y == '-' ? -x : x;
}
inline void prt(li x){
	if(x >= 10) prt(x / 10);
	pc(x % 10 + '0');
}
inline void print(li x){
	if(x < 0) pc('-'),x = -x;
	prt(x);
}
inline void file(char *s){
	char c[50];
	sprintf(c,"%s.in",s);
	freopen(c,"r",stdin);
	sprintf(c,"%s.out",s);
	freopen(c,"w",stdout);
}
int n;
#define M 5005
li a[M];
int fa[M],ls[M],rs[M];
vector<li> f[M][M];
vector<int> ch[M];
int dpt[M];
bool inzs[M][M];
li vl[M],vr[M],dl[M],dr[M],dvl[M],dvr[M];
void dfs(int x){
	if(!x) return;
	int lc = ls[x],rc = rs[x];
	dfs(lc);dfs(rc);
	dpt[x] = max(dpt[lc],dpt[rc]) + 1;
	int i = 0,j = 0,k,l,o,p,q,d,dt = dpt[x];
	li ax = a[x];
	ch[x].pb(x);
	int szl = ch[lc].size(),szr = ch[rc].size();
	while(i < szl && j < szr){
		if(a[ch[lc][i]] < a[ch[rc][j]]) ch[x].pb(ch[lc][i++]);
		else ch[x].pb(ch[rc][j++]);
	}
	while(i < szl) ch[x].pb(ch[lc][i++]);
	while(j < szr) ch[x].pb(ch[rc][j++]);
	for(i = 1;i < ch[x].size() && ax > a[ch[x][i]];++i) swap(ch[x][i - 1],ch[x][i]);
	for(i = 0;i < ch[x].size();++i) inzs[x][ch[x][i]] = 1;
	
	if(!lc) f[x][x].pb(ax);
	
	if(lc && !rc){
		f[x][x].resize(dt);
		for(d = 0;d < dt;++d) f[x][x][d] = 1e18l;
		for(i = 0;i < szl;++i){
			p = ch[lc][i];
			for(d = 1;d <= f[lc][p].size();++d){
				f[x][x][d] = min(f[x][x][d],f[lc][p][d - 1] + ax);
			}
		}
		for(i = 0;i < szl;++i){
			p = ch[lc][i];
			f[x][p].pb(1e18l);
			for(d = 0;d < f[lc][p].size();++d){
				f[x][p][0] = min(f[x][p][0],f[lc][p][d] + ax * (d + 1) + a[p]);
			}
		}
	}
	
	if(rc){
		memset(vl,0x1f,sizeof(vl));
		memset(vr,0x1f,sizeof(vr));
		memset(dl,0x1f,sizeof(dl));
		memset(dr,0x1f,sizeof(dr));
		memset(dvl,0x1f,sizeof(dvl));
		memset(dvr,0x1f,sizeof(dvr));
		for(i = 0;i < szl;++i){
			p = ch[lc][i];
			for(d = 0;d < f[lc][p].size();++d){
				vl[p] = min(vl[p],f[lc][p][d] + ax * (d + 1));
				dl[d] = min(dl[d],f[lc][p][d]);
			} 
			for(d = 0;d < dpt[rc];++d){
				dvr[d] = min(dvr[d],vl[p] + a[p] * (d + 1));
			}
		}
		for(j = 0;j < szr;++j){
			q = ch[rc][j];
			for(d = 0;d < f[rc][q].size();++d){
				vr[q] = min(vr[q],f[rc][q][d] + ax * (d + 1));
				dr[d] = min(dr[d],f[rc][q][d]);
			}
			for(d = 0;d < dpt[lc];++d){
				dvl[d] = min(dvl[d],vr[q] + a[q] * (d + 1));
			}
		}
		
		
		for(i = 0;i < szl;++i){
			p = ch[lc][i];
			f[x][p].resize(dpt[rc] + 1);
			f[x][p][0] = 1e18l;
			for(d = 1;d <= dpt[rc];++d){
				f[x][p][d] = vl[p] + dr[d - 1] + a[p];
			}
		}
		for(j = 0;j < szr;++j){
			q = ch[rc][j];
			f[x][q].resize(dpt[lc] + 1);
			f[x][q][0] = 1e18l;
			for(d = 1;d <= dpt[lc];++d){
				f[x][q][d] = vr[q] + dl[d - 1] + a[q];
			}
		}
		
		
		f[x][x].resize(dt);
		for(d = 0;d < dt;++d) f[x][x][d] = 1e18l;
		for(i = 0;i < szl;++i){
			p = ch[lc][i];
			li tmp = 1e18l;
			for(o = 0;o < dpt[rc];++o) tmp = min(tmp,dr[o] + a[p] * (o + 1));
			for(d = 0;d < f[lc][p].size();++d){
				f[x][x][d + 1] = min(f[x][x][d + 1],f[lc][p][d] + tmp + a[x]);
				f[x][p][0] = min(f[x][p][0],dvl[d] + f[lc][p][d] + a[p]);
			}
		}
		for(j = 0;j < szr;++j){
			q = ch[rc][j];
			li tmp = 1e18l;
			for(o = 0;o < dpt[lc];++o) tmp = min(tmp,dl[o] + a[q] * (o + 1));
			for(d = 0;d < f[rc][q].size();++d){
				f[x][x][d + 1] = min(f[x][x][d + 1],f[rc][q][d] + tmp + a[x]);
				f[x][q][0] = min(f[x][q][0],dvr[d] + f[rc][q][d] + a[q]);
			}
		}
	}
}
int main(){
	//file("mis");
	int i,k;
	n = read();
	for(i = 1;i <= n;++i) a[i] = read();
	for(i = 2;i <= n;++i){
		k = read();fa[i] = k;
		if(ls[k]) rs[k] = i;
		else ls[k] = i;
	}
	dfs(1);
	li ans = 1e18l;
	for(i = 0;i < dpt[1];++i) ans = min(ans,f[1][1][i] + a[1] * (i - 1));
	print(ans);pc('\n');
	return 0;
}
liuzhangfeiabc
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2018NOIP天津赛区选手代码
03-15
2018NOIP天津赛区选手代码
[省选联考 2022] 预处理器（民间数据）
绿树公司的博客
04-18 208
[省选联考 2022] 预处理器（民间数据） 题目背景 1 s 512 MiB 有问题请找 tiger2005。目前的民间数据放在附件了。 为了达到强力 hack 效果，我们选择直接对 hack 进行追加操作。这将会导致测试点数量多于 10 个，请大家注意。 感谢 IdnadRev 提供 Hack 数据 #11。 题目描述 宏是 C/C++ 语言的一项特性，它根据预先定义的规则进行文本替换（也被称为 “宏展开”），能够实现定义常量、简化代码重复输入等功能。例如： #define PI 3.14159 dou
NOI2022联合省选 题解
liuzhangfeiabc的博客
04-17 2574
NOI2022联合省选 题解 目前只有day2
【省选联考 2022 D2T2】序列变换 题解
Cherrt的博客
04-22 620
文章目录DescriptionSolutionPart 1: 抽象为树上问题与初步解决Part 2: x=0,y=1Part 3: x=1,y=1Part 4: x=1,y=0 Description 传送门 Solution Part 1: 抽象为树上问题与初步解决 首先，我们对括号序列建出一棵树。接着，我们尝试将两种操作转化为树上的操作。 定义位置 iii 的匹配点 pip_ipi​ 为与第 iii 个括号匹配的括号的位置，则考虑每一个满足 si=s_i=si​=( 对应的 [i,pi][i,p_i]
【luogu P8293】[省选联考 2022] 序列变换（贪心）（分类讨论）
欢迎您的光顾awa
05-01 253
给你一个括号序列，每次你可以把 p(A)(B)q 的串变成 P(A()B)q。 你还可以不用花费交换任意两个相邻合法括号序列的位置。 其中 A,B 是合法括号序列，p,q 可以不是。 然后每个左括号有费用，每次边的费用是左边左括号费用和右边左括号费用分别乘 x,y 的和。 x,y 都是 0/1，要你用最小费用使得不能找到可以变换的串。
2016山东省NOI省队集训DAY6 题目题解数据标程 NOIP省选集训
04-18
2016山东省NOI省队集训DAY6 题目题解数据标程 NOIP省选集训
NOI2019省选
04-13
2019年12省（黑龙江、吉林、辽宁、山西、河北、山东、河南、安徽、江苏、上海、湖北、重庆）省选，Day1+Day2（4月6~7日）原题，无测试数据，无样例文件。
NOI 2021 统一省选
04-16
CCF NOI 2021 统一省选（含有 A、B 卷题目）
2016山东省NOI省队集训DAY3 题目题解数据标程 省选模拟赛
04-18
2016山东省NOI省队集训DAY3 题目题解数据标程 省选集训模拟赛
2016山东省NOI省队集训DAY5 题目题解数据标程 NOIP省选集训
04-18
2016山东省NOI省队集训DAY5 题目题解数据标程 NOIP省选集训
NOI2017退役记
commonc的博客
07-24 7400
前文链接：去年NOI2016的游记，对比产生美以后会有一个长篇的OI生涯回忆录，敬请关注DAY-??有一些故事背景的前提下：我："我赌我国赛肯定不能A题"ljss：“那咱俩赌啊，我赌你能A题，赌十块的”我：“赌20吧，我赌我A不了题”ljss：“行啊”我：“那我要是全都爆零了是不是也算我没A题啊”ljss：“那当然算了”于是我们找了yycc和alone_wolf当公证人DAY?ljss：“我赌你这...
NOI Online2022 入门组题解
NOI RP++
04-05 4417
2022-03-26，CCF官方举行了NOI Online能力测试，因为我是个蒟蒻，才考了235分 ，最后一题也是看了题解才写出来 题目： T1王国比赛 T2数学游戏 T3字符串 题解可以看 这个 就给大家展示一下我的满分代码吧 T1 #include<bits/stdc++.h> #define FO(x)\ freopen(#x".in","r",stdin);\ freopen(#x".out","w",stdout); using namespace std; int k[1100][
联合省选2022游记
StaroForgin的博客
04-18 1642
这就是我们的现实。
[省选联考 2022] 序列变换（暂无数据）
绿树公司的博客
04-18 275
[省选联考 2022] 序列变换（暂无数据） 题目背景 1 s 512 MiB 题目描述 你手里有一个长度为 2n2 n2n 的合法括号序列 sss。sss 的每一个左括号有一个权值。 在你眼中，不同的括号序列带来的视觉美感不尽相同。因此，你对具有某一种结构的括号序列特别喜欢，而讨厌具有其他一些结构的括号序列。你希望对 sss 进行一些变换，以消除掉一些自己不喜欢的结构。 具体而言，形如 (A()B)\texttt{(A()B)}(A()B)（其中 A\texttt{A}A、B\texttt{B}B 均为合
NOI2022 题解
liuzhangfeiabc的博客
08-30 2591
这里是liuzhangfeiabc的NOI2022题解。 目前更新了两天的t1和t2，两天的t3实在没心情去研究了…… 注：非官方题解，不保证绝对正确，做法不一定是官方做法。
NOI2020联合省选A题组，中山纪中游记
A1847225889的博客
06-22 773
Day1 饭堂的早餐特意点得丰盛了一点，还买了咖啡，感觉精神饱满。 比赛前以为可以试机，然后迫不及待地进了考场，接着才知道不能试机…… 桌面上放着面包和士力架。 在座位上冥想，想起忘了复习FFT三次变两次的优化怎么打，于是在脑中略微复习了一下。 比赛开始，密码的最后一位是‘y’，监考员写的看起来像‘g’，害得我输入了好几遍。 开始看题。 正式比赛不像平常那样那么浪，于是BFS式想题，先粗略看三题题意，再细看题意，将理解的题意和想到的模型写在草稿纸上。 开始干T1。 略微化了下模型发现就是个线段树二分，O(n
NOIonline 2022
lahlah的博客
04-02 558
降大智了属于是 考后光速补完 A 丹钓战 sb题,首先可以处理出一个L[i]L[i]L[i]表示iii这个点最多可以向前一路pop到那个点 然后就可以用扫描线+树状数组简单维护了 处理L[i]L[i]L[i]可以简单的用个单调栈模拟题目，但是我考场上降智了，去找了个性质，首先可以二分+st表找出iii前面一串比a[i]a[i]a[i]小的区间[l,i][l,i][l,i]，然后找这个区间中最小的L[pos]L[pos]L[pos] 发现a[L[pos]−1]a[L[pos]-1]a[L[pos]−1]一定是
【NOI2022省选挑战赛 Contest10 B】早该砍砍了（笛卡尔树）（DP）
欢迎您的光顾awa
03-21 460
给你一个序列，保证任意两个位置的数都不相同，你每次可以选择一个区间把里面的数都变成其中的最小值，问你操作若干次（可以不操作）之后序列的样子有多少种。
浙江省选day2第二天
Greninja的博客
04-27 772
ZJOI2017day2第二天
[noi2022] 二次整数规划问题
最新发布
09-18
二次整数规划问题是指在给定的约束条件下，要求在整数解中最小化二次目标函数的优化问题。它是数学规划中一类重要的问题。 二次整数规划可以形式化表示为： minimize f(x) = ½ x^T Q x + c^T x ...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值