24.3.10周报( 线性筛,区间dp, 重载运算符,二维全局vector）

让我康米

已于 2024-05-10 21:48:10 修改

阅读量956

点赞数 25

文章标签：算法图论数据结构

于 2024-03-10 21:06:54 首次发布

本文链接：https://blog.csdn.net/breath_of_yuan/article/details/136462629

版权

星期一：

cf round930 div2 C题交互：cf传送门

思路：先找最大值，再找和最大值异或最大的值，异或值最大，二进制中即为0对应1，按位或的结果即为1，所以也应是最大，按位或结果相同时取更小值

代码如下：

ll n;
int m,ma,mi;
int ask(int a,int b,int c,int d){
	cout << "? " << a << " " << b << " " << c << " " << d << endl;
	char x; cin >> x;
	if(x=='<') return 0;
	if(x=='=') return 1;
	return 2;
}
void solve(){
	cin >> n;
	for(int i=0;i<n;i++)
		if(!ask(ma,ma,i,i)) ma=i;
	for(int i=0;i<n;i++){
		int tmp=ask(ma,mi,ma,i);
		if(!tmp) mi=i;
		else if(tmp==1)
			if(!ask(i,i,mi,mi)) mi=i;
	}
	cout << "! " << ma << " " << mi << endl;
	ma=mi=0;
}

再敲遍线性筛：复杂度O(n）

const int N=1e5+10;
ll n;
int p[N],cnt;
bool vi[N];
void getp(int x){
	for(int i=2;i<=x;i++){
		if(!vi[i]) p[++cnt]=i;
		for(int j=1;1ll*i*p[j]<=x;j++){
			vi[i*p[j]]=1;
			if(i%p[j]==0) break;
		}
	}
}

区间dp板子石子合并

代码如下：

ll n;
int a[310],sum[310];
ll dp[310][310];
void solve(){
	cin >> n;
	memset(dp,0x3f,sizeof dp);
	for(int i=1;i<=n;i++){
		cin >> a[i];
		sum[i]=sum[i-1]+a[i];
		dp[i][i]=0;
	}
	for(int i=2;i<=n;i++){
		for(int j=1;j+i-1<=n;j++){
			int r=j+i-1;
			for(int k=j;k<r;k++){
				dp[j][r]=min(dp[j][k]+dp[k+1][r]+sum[r]-sum[j-1],dp[j][r]);
			}
		}
	}
	cout << dp[1][n];
}

星期二：

天梯选拔赛，手速可以，但上限是卡在那了

补一题：

思路：dp，左区间从1到n，右区间从n到 1，要注意不合法的区间状态（我就wa在这了），比如n为奇数，dp[1][2]即不合法，因为每次删去偶数个，所以遇到不合法区间continue

代码如下：

ll n;
int a[2020];
ll dp[2020][2020];
ll ans;
void solve(){
	cin >> n;
	for(int i=1;i<=n;i++) cin >> a[i];
	for(int i=1;i<=n;i++){
		for(int j=n;j>i;j--){
			if(((j-i)&1)==(n&1)) continue;
			ll r1=0,r2=0,r3=0;
			if(i>2) r1=dp[i-2][j]+(a[i-1]+a[i-2])%3;
			if(j<n-1) r2=dp[i][j+2]+(a[j+1]+a[j+2])%3;
			if(i>1 && j<n) r3=dp[i-1][j+1]+(a[i-1]+a[j+1])%3;
			dp[i][j]=max({r1,r2,r3,dp[i][j]});
			if(i==j-1) ans=max(dp[i][j]+(a[i]+a[j])%3,ans);
			ans=max(dp[i][j],ans);
		}
	}
	cout << ans;
}

晚cf C没写出来，且手速慢了，掉一点分，但离上蓝还会远吗

星期三：

补选拔赛二题：

思路：分解质因数，把n分解为p1^a1 * p2^a2 * p3^a3...答案为(a1+1)*(a2+1)*(a3+1)...，n<=1e8，所以线性筛筛到根号1e8即1e4多一点(1e4+10)就够了

代码如下：

const int N=1e6+10;
const int mod=998244353;
ll n;
int p[N],cnt;
bool vi[N];
vector<int>ve;
bool isp(int x){
	for(int i=2;i<=x/i;i++)
		if(x%i==0) return 0;
	return 1;
}
void solve(){
	for(int i=2;i<=1e4+10;i++){
		if(!vi[i]) p[++cnt]=i;
		for(int j=1;1ll*i*p[j]<=1e4+10;j++){
			vi[1ll*i*p[j]]=1;
			if(i%p[j]==0) break;
		}
	}
	int t; cin >> t;
	while(t--){
		cin >> n;
		for(int i=1;p[i]<=n/p[i] && n!=1;i++){
			int tmp=0;
			while(n%p[i]==0){
				n/=p[i];
				tmp++;
			}
			if(tmp) ve.push_back(tmp+1);
		}
		if(n!=1) ve.push_back(2);
		ll ans=1;
		for(auto i:ve) ans*=i;
		cout << ans << "\n";
		ve.clear();
	}
}

cf round 932 div2 C题：

思路：按b排序，然后跑一遍n^2，b的权值能一次算出，范围内的a塞优先队列里，时间超了再从优先队列里排出

代码如下：

ll n;
int k,t;
pair<ll,ll> ba[2020];
ll ans,aa,bb,minn=INT_MAX;
void solve(){
	cin >> n >> t;
	for(int i=1;i<=n;i++){
		cin >> ba[i].second >> ba[i].first;
		minn=min(ba[i].second,minn);
	}
	if(t>=minn) ans=1;
	sort(ba+1,ba+n+1);
	for(int i=1;i<=n;i++){
		priority_queue<int>pq;                    //大根堆
		aa=0;
		for(int j=i;j<=n;j++){
			bb=ba[j].first-ba[i].first;
			if(bb>=t) break;
			pq.push(ba[j].second);
			aa+=ba[j].second;
			while(aa+bb>t){
				aa-=pq.top();
				pq.pop();
			}
			if(pq.size()>ans) ans=pq.size();
		} 
	}
	cout << ans << "\n";
	ans=0;minn=INT_MAX;
}

学了下kmp，没学明白，板子先放这 O(n+m) 洛谷传送门

const int N=1e6+10;
const int mod=998244353;
ll n;
int m;
string s,p;
int ne[N];
void solve(){
	cin >> s >> p;
	n=p.size(),m=s.size();
	s=" "+s,p=" "+p;
	ne[1]=0;
	for(int i=2,j=0;i<=n;i++){
		while(j && p[i]!=p[j+1]) j=ne[j];
		if(p[i]==p[j+1]) j++;
		ne[i]=j;
	}
	for(int i=1,j=0;i<=m;i++){
		while(j && s[i]!=p[j+1]) j=ne[j];
		if(s[i]==p[j+1]) j++;
		if(j==n) cout << i-n+1 << "\n";
	}
//	for(int i=1;i<=n;i++) cout << ne[i] << " ";
}

补选拔赛三题（最STL的一集）：

思路：将冰和鞋分别按光滑度和防滑度从小到大排序,set维护鞋能走的冰,multiset维护任意两冰面间的最小距离,判断鞋能否通过,STL操作较多,对其熟练度要求比较高

代码如下：

const int N=1e6+10;
const int mod=998244353;
ll n;
int m;
PII f[N];
struct nod{
	int k,s,idx;
	bool if1;
}b[N];
set<int>se;
multiset<int>ms;
bool cmp1(nod a,nod b){
	return a.k<b.k;
}
bool cmp2(nod a,nod b){
	return a.idx<b.idx;
}
void solve(){
	cin >> n >> m;
	for(int i=1;i<=n;i++){
		cin >> f[i].first;
		f[i].second=i;
	}
	sort(f+1,f+n+1);
	for(int i=1;i<=m;i++){
		cin >> b[i].k >> b[i].s;
		b[i].idx=i;
	}
	sort(b+1,b+m+1,cmp1);
	int cnt=1;
	while(cnt<=n && !f[cnt].first){
		se.insert(f[cnt].second);
		cnt++;
	}
	for(auto i=se.begin();next(i)!=se.end();i=next(i)){
		ms.insert(*(next(i))-*(i));
	}
	for(int i=1;i<=m;i++){
		while(f[cnt].first<=b[i].k && cnt<=n){
			int l=*prev(se.upper_bound(f[cnt].second));
			int r=*(se.upper_bound(f[cnt].second));
			se.insert(f[cnt].second);
			ms.erase(ms.find(r-l));
			ms.insert(f[cnt].second-l),ms.insert(r-f[cnt].second);
			cnt++;
		}
		if(b[i].s>=*(ms.rbegin())) b[i].if1=1;
		else b[i].if1=0;
	}
	sort(b+1,b+m+1,cmp2);
	for(int i=1;i<=m;i++) cout << b[i].if1 << "\n";
}

星期四：

补cf round932 D

思路：答案为总数-(y+x)在s的对数-(y-x)在s的对数+(y+x)和(y-x)都在s的对数，最后一个需要把奇偶分开算

代码如下：

const int N=1e6+10;
ll n;
ll c,s[N];
int num1,num2;
void solve(){
	cin >> n >> c;
	for(int i=1;i<=n;i++){
		cin >> s[i];
		if(s[i]&1) num1++;
		else num2++;
	}
	ll ans=(c+2)*(c+1)/2;
	for(int i=1;i<=n;i++){
		ans-=s[i]/2+1;
		ans-=c-s[i]+1;
	}
	ans+=1ll*num1*(num1+1)/2,ans+=1ll*num2*(num2+1)/2;
	cout << ans << "\n";
	num1=num2=0;
}

接着看了下E，思路懂了，具体的代码实现感觉还是欠缺了点功夫

牛客练习赛122 C题马走日牛客传送门

思路：3之内的数据范围特判下，更大就是n*m，wa在n或m为2时的特判没判对，粗糙了

D题 ⚪ :

思路：拆⚪为链，逆向，区间dp，buff叠满了这题，代码是敲了，但不能说懂了

把⚪拆为两倍长的链，链中每个长为n的区间都可视为一个环，求删除的最小代价即为不重合的最大价值，然后区间dp先枚举终点，再枚举起点，再枚举起点的边，感觉有点玄

代码如下:

const int N=2e3+10;
const int mod=1e9+7;
ll n;
int m;
vector<PII>ve[2*N];
ll dp[2*N][2*N],sum;
void solve(){
	cin >> n >> m;
	for(int i=1;i<=m;i++){
		int x,y,w;
		cin >> x >> y >> w;
		sum+=w;
		if(x>y) swap(x,y);
		ve[x].push_back({y,w});
		ve[y].push_back({x+n,w});
	}
	for(int i=2;i<=2*n;i++){
		for(int j=i-1;j;j--){
			if(i-j+1>n) break;
			dp[j][i]=dp[j+1][i];
			for(auto [v,w]:ve[j]){
				if(v>i) continue;
				ll sum=w;
				if(v>j+2) sum+=dp[j+1][v-1];
				if(i>v+1) sum+=dp[v+1][i];
				dp[j][i]=max(sum,dp[j][i]);
			}
		}
	}
	ll cnt=0;
	for(int i=1;i<=n+1;i++) cnt=max(dp[i][i+n-1],cnt);
	cout << sum-cnt;
}

星期五：

牛客小白月赛，有道题赛时没做出来很无语

我不是大富翁（牛客传送门

题不难，坑有点多，比较吃细节，和cf风格不太一样，平时多打点牛客也挺好的

多重映射（牛客传送门

思路：并查集，但又没全并，倒着并，并且用不着find找祖宗

代码如下：

const int N=2e6+10;
const int mod=1e9+7;
ll n;
int m;
int a[N];
map<int,int>fa;
PII op[N];
void solve(){
	cin >> n >> m;
	for(int i=1;i<=n;i++){
		cin >> a[i];
		fa[a[i]]=a[i];
	}
	for(int i=1;i<=m;i++){
		cin >> op[i].first >> op[i].second;
		fa[op[i].first]=op[i].first;
		fa[op[i].second]=op[i].second;
	}
	for(int i=m;i;i--){
		fa[op[i].first]=fa[op[i].second];
	}
	for(int i=1;i<=n;i++){
		cout << fa[a[i]] << " ";
	}
	cout << "\n";
	fa.clear();
}

星期六：

天体选拔赛 round2 打拉了，但我要嘴硬一下这题不太好

精灵融合

思路：一道普普通通树形dfs，但我算经验时漏了等级为1所需经验值为0的情况，把经验值默认为1了，真太粗糙了，反思，woyouzui

幸运号码

思路：对每个数枚举前缀作为Si和后缀作为Sj的情况，不过这里枚举Si时j跑完，但枚举为Sj时j不到0的原因没懂，是处理和自己匹配的情况吗？

代码如下：

const int N=2e6+10;
const int mod=1e9+7;
ll n;
int sum[N];
int cnt[100][2];
string s[N];
ll ans;
int ask(string s){
	int res=0;
	for(int i=0,j=s.size();i<j;i++) res+=s[i]-'0';
	return res;
}
void solve(){
	cin >> n;
	for(int i=1;i<=n;i++){
		cin >> s[i];
		sum[i]=ask(s[i]);
		cnt[sum[i]][s[i].size()&1]++;
	}
	for(int i=1;i<=n;i++){
		int l=0;
		for(int j=0;j<s[i].size();j++){
			l+=s[i][j]-'0';
			int r=sum[i]-l;
			if(l<=r) continue;
			ans+=cnt[l-r][s[i].size()&1];
		}
	}
	for(int i=1;i<=n;i++){
		int r=0;
		for(int j=s[i].size()-1;j>0;j--){
			r+=s[i][j]-'0';
			int l=sum[i]-r;
			if(r<=l) continue;
			ans+=cnt[r-l][s[i].size()&1];
		}
	}
	cout << ans;
}

周日：

昨晚初试atcoder ABC344，补道D ac传送门

思路：比较新的字符串dp，dp[i,j]表示考虑到前i个包能组成T前j个字符的最小花费，这里用find函数来continue被坑了

代码如下：

ll n;
vector<string>ve[110];
string t;
ll dp[110][110];
void solve(){
	cin >> t >> n;
	for(int i=1;i<=n;i++){
		int a; cin >> a;
		while(a--){
			string s; cin >> s;
			ve[i].push_back(s);
		}
	}
	for(int i=0;i<=n+1;i++){
		for(int j=1;j<=(int)t.size();j++) dp[i][j]=1e9;
	}
	for(int i=1;i<=n;i++){
		for(auto j:ve[i]){
//			if(t.find(j)==string::npos) continue;
			int len=j.size();
			for(int k=1;k<=(int)t.size();k++){
				dp[i][k]=min(dp[i-1][k],dp[i][k]);
				if(k<len) continue;
				if(t.substr(k-len,len)!=j) continue;
				dp[i][k]=min(dp[i-1][k-len]+1,dp[i][k]);
			}
		}
	}
	if(dp[n][t.size()]==1e9) cout << "-1";
	else cout << dp[n][t.size()];
}

开二维不定长数组，防止爆内存

注意后面一定要用a=，命名在括号前会成局部数组

vector<vector<int>>a;                 //全局声明

a=vector<vector<int>>(n+1,vector<int>(m+1,0));

思路：贪心bfs，贪心用优先队列实现

代码如下：

ll n;
int m,x,y,t;
int dx[]={0,0,-1,1};
int dy[]={1,-1,0,0};
ll ans;
struct nod{
	int x,y,w;
	bool operator<(const nod &b) const{   //重载运算符来排序
		return w>b.w;
	}
};
vector<vector<int>>a;            //不定长数组的全局声明
vector<vector<bool>>vi;
vector<vector<bool>>if1;
void bfs(){
	priority_queue<nod>pq;       //优先队列 instead of 队列
	pq.push({x,y,a[x][y]});
	vi[x][y]=1; bool if2=0;
	while(!pq.empty()){
		nod tp=pq.top(); pq.pop();
		if(ans>=tp.w || !if2) ans+=tp.w,if2=1;
		else break;
		for(int i=0;i<4;i++){
			int xx=tp.x+dx[i],yy=tp.y+dy[i];
			if(xx<1 || xx>n || yy<1 || yy>m || vi[xx][yy]) continue;
			vi[xx][yy]=1;
			if(if1[xx][yy]) ans+=a[xx][yy],pq.push({xx,yy,0});
			else pq.push({xx,yy,a[xx][yy]});
		}
	}
	
}
void solve(){
	cin >> n >> m >> x >> y >> t;
	a=vector<vector<int>>(n+1,vector<int>(m+1,0));      //定义要用等号，不然会成局部数组
	vi=vector<vector<bool>>(n+1,vector<bool>(m+1,0));
	if1=vector<vector<bool>>(n+1,vector<bool>(m+1,0));
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++)
			cin >> a[i][j];
	}
	for(int i=1;i<=t;i++){
		int x,y;
		cin >> x >> y;
		if1[x][y]=1;
	}
	bfs();
	cout << ans;
}

这里的重载运算符的逻辑捋捋

因为优先队列默认大根堆，所以在排序时a<b代表着b放a前，所以我们要将<符号重载为>符号的含义，这样优先队列自动排序时，得出a<b，实际上a是大于b的，这样在b放a前面后，就成了小根堆

晚上打打牛客和ac，牛客最后一道不会，ac第一道就不会

让我康米

关注

25
点赞
踩
22

收藏

觉得还不错? 一键收藏
0
评论
24.3.10周报( 线性筛,区间dp, 重载运算符,二维全局vector）

星期一：cf round930 div2 C题交互：cf传送门思路：先找最大值，再找和最大值异或最大的值，异或值最大，二进制中即为0对应1，按位或的结果即为1，所以也应是最大，按位或结果相同时取更小值代码如下：ll n;int m,ma,mi;int ask(int a,int b,int c,int d){ cout << "? " << a << " " << b << " " << c << "
复制链接

扫一扫