NOIP2012提高组解析

本文解析了NOIP2012提高组的竞赛题目,涉及维吉尼亚解密、国王游戏、开车旅行等算法问题。通过参考程序展示了如何运用数学和逻辑解决这些挑战,包括高精度计算、预处理、模拟、同余方程、二分搜索、最短路径及LCA(最近公共祖先)等算法思想。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目见此



day1

维吉尼亚解密:

从图中可以看出,就是将明文的字母,向前移动对应密钥字母在26字母中的排名数个格子就是原文。


参考程序:

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cctype>
using namespace std;
char key[110];
char dark[1100];
int main(){
	freopen("vigenere.in","r",stdin);
	freopen("vigenere.out","w",stdout);
	scanf("%s%s",&key,&dark);
	int j=0;
	int l=strlen(key);
	for (int i=0;i<strlen(dark);i++){
		int ch=tolower(key[j])-'a';
		char cmd=dark[i]>'Z'?'a':'A';
		int x=(dark[i]-cmd-ch+26)%26;
		printf("%c",cmd+x);
		j=(j+1)%l;
	}
	return 0;
}


国王游戏:

其实就是按照大臣左右手乘积的大小排序加上一个高精度。

证明要点:从最后一个人开始推,假设其与前面的一人交换,则导致结果变大。


参考程序:

#include<cstdio>
#include<algorithm>
#include<cstring>
#define maxn 11000
using namespace std;
struct King{
	int l,r;
	bool operator < (const King &b)const{
		return l*r<b.l*b.r;
	}
}a[maxn];
int t[maxn];
int ans[maxn];
int f[maxn];
int len=0,alen=0,n;
void Div(int x){
	int y=0;
	memset(t,0,sizeof(t));
	for (int i=len-1;i>=0;i--){
		y=y*10+f[i];
		t[i]=y/x;
		y%=x;
	}
	int j=len;
	while (!t[j-1])j--;
	if (alen<j){
		alen=j;
		for (int i=0;i<j;i++)
			ans[i]=t[i];
	}else if (alen==j){
		bool flag=true;
		for (int i=0;i<j;i++)
			if (ans[i]<t[i]){flag=false;break;}
		if (!flag)for (int i=0;i<j;i++)ans[i]=t[i];
	}
}
void mult(int x){
	for (int i=0;i<len;i++)
		f[i]*=x;
	for (int i=0;i<len;i++){
		f[i+1]+=f[i]/10;
		f[i]%=10;
	}
	while (f[len]>0){
		f[len+1]=f[len]/10;
		f[len]=f[len]%10;
		len++;
	}
}
int main(){
	freopen("game.in","r",stdin);
	freopen("game.out","w",stdout);
	scanf("%d",&n);
	for (int i=0;i<=n;i++)
		scanf("%d%d",&a[i].l,&a[i].r);
	sort(a+1,a+n+1);
	f[0]=a[0].l;
	len=1;
	for (int i=1;i<=n;i++){
		Div(a[i].r);
		mult(a[i].l);
	}
	for (int i=alen-1;i>=0;i--)
		printf("%d",ans[i]);
	return 0;
}


开车旅行;

预处理出每个城市向后走2^j次后到达哪个城市,然后就是模拟。


参考程序:

#include<cstdio>
#include<algorithm>
#include<set>
#define maxn 110000
using namespace std;
struct City{
	int h,id;
	bool operator < (const City &b) const{
		return h<b.h;
	}
}c[maxn];
struct Temp{
	int pos,dis;
	Temp(int pos,int dis):pos(pos),dis(dis){}
	Temp(){}
	bool operator < (const Temp &b) const {
		if (dis != b.dis)return dis<b.dis;
		return c[pos].h<c[b.pos].h;
	}
}tmp[5];
int M,n;
int f[maxn][20];
long long fa[maxn][20];
long long fb[maxn][20];
set<City> s;
int nexta[maxn];
int nextb[maxn];
void prepare(int i){
	set<City> :: iterator it=s.find(c[i]); 
	int add=0;
	if (it!=s.begin()){
		--it;tmp[++add]=Temp(it->id,abs(c[i].h-it->h));
		if (it!=s.begin()){
			--it;tmp[++add]=Temp(it->id,abs(c[i].h-it->h));++it;
		}
		++it;
	}
	if ((++it)!=s.end()){
		tmp[++add]=Temp(it->id,abs(c[i].h-it->h));
		if ((++it)!=s.end()){
			tmp[++add]=Temp(it->id,abs(c[i].h-it->h));
		}
	}
	sort(tmp+1,tmp+add+1);
	nextb[i]=tmp[1].pos;
	if (add!=1)nexta[i]=tmp[2].pos;
}
void solve(int begin,int x,long long &ta,long long &tb){
	for (int j=19;j>=0;j--)
		if (f[begin][j] && fa[begin][j]+fb[begin][j]<=x){
			ta+=fa[begin][j];tb+=fb[begin][j];
			x-=fa[begin][j]+fb[begin][j];begin=f[begin][j];
		}
	int p1=nexta[begin];
	if (!p1)return;
	int dis=abs(c[p1].h-c[begin].h);
	if (dis<=x)ta+=dis;
}
int main(){
	freopen("drive.in","r",stdin);
	freopen("drive.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;i++){
		scanf("%d",&c[i].h);
		c[i].id=i;
	}
	for (int i=n;i>0;i--){
		s.insert(c[i]);
		if (i!=n)prepare(i);
	}
	for (int i=1;i<=n;i++){
		int p1=nexta[i];
		int p2=nextb[nexta[i]];
		fa[i][0]=p1?abs(c[i].h-c[p1].h):0;
		fb[i][0]=p2?abs(c[p1].h-c[p2].h):0;
		f[i][0]=p2;
	}
	for (int j=1;j<20;j++)
		for (int i=1;i<=n;i++){
			f[i][j]=f[f[i][j-1]][j-1];
			fa[i][j]=fa[i][j-1]+fa[f[i][j-1]][j-1];
			fb[i][j]=fb[i][j-1]+fb[f[i][j-1]][j-1];
		}
	int x0;
	scanf("%d",&x0);
	long long ansa=1e15,ansb=0ll;
	int num=0;
	for (int i=1;i<=n;i++){
		long long sa=0ll,sb=0ll;
		solve(i,x0,sa,sb);
		if (sb && (!num || sa*ansb<sb*ansa)){
			ansa=sa;ansb=sb;
			num=i;
		}
	}
	printf("%d\n",num);
	scanf("%d",&M);
	while (M--){
		int begin,x0;
		scanf("%d%d",&begin,&x0);
		long long sa=0ll,sb=0ll;
		solve(begin,x0,sa,sb);
		printf("%lld %lld\n",sa,sb);
	}
	return 0;
}



day2:

同余方程:

模的原理


参考程序:

#include<cstdio>
#include<algorithm>
using namespace std;
int n,m;
int exgcd(int n,int m,int &x,int &y){
	if (!m){
		x=1;y=0;
		return n;
	}
	int k=exgcd(m,n%m,y,x);
	y-=n/m*x;
	return k;
}
int main(){
	freopen("mod.in","r",stdin);
	freopen("mod.out","w",stdout);
	scanf("%d%d",&n,&m);
	int b=m,x,y;
	int d=exgcd(n,m,x,y);
	printf("%d",(x+b)%b);
	return 0;
}



借教室:

IMOS算法,二分即可。当然线段树也可以。


参考程序:

#include<cstdio>
#include<algorithm>
#define maxn 1100000
using namespace std;
int x[maxn],d[maxn],s[maxn],t[maxn];
int a[maxn];
int n,m;
int main(){
	freopen("classroom.in","r",stdin);
	freopen("classroom.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		scanf("%d",&x[i]);
	for (int i=1;i<=m;i++)
		scanf("%d%d%d",&d[i],&s[i],&t[i]);
	int l=0,r=m+1;
	int last=0;bool flag;
	int mid=0;
	while (l<r){
		mid=(l+r)>>1;
		if (last<mid){
			for (int i=last+1;i<=mid;i++){
				a[s[i]]=a[s[i]]+d[i];
				a[t[i]+1]=a[t[i]+1]-d[i];
			}
		}else{
			for (int i=mid+1;i<=last;i++){
				a[s[i]]=a[s[i]]-d[i];
				a[t[i]+1]=a[t[i]+1]+d[i];
			}
		}
		last=mid;
		flag=true;
		int sum=0;
		for (int i=1;i<=n;i++){
			sum=sum+a[i];
			if (sum>x[i]){
				flag=false;
				break;
			}
		}
		if (flag)l=mid+1;else r=mid;
	}
	if (mid==m && flag)printf("0");else{
		printf("-1\n%d",l);
	}
	return 0;
}

补上线段树的:

#include<cstdio>
#include<algorithm>
#define maxn 1110000
using namespace std;
int s[4*maxn];
long long sign[4*maxn];
int d[maxn],b[maxn],t[maxn];
int a[maxn];
int n,m;
void update(int no){
	s[no]-=sign[no];
	sign[2*no]+=sign[no];
	sign[2*no+1]+=sign[no];
	sign[no]=0;
}
void add(int no,int l,int r,int ql,int qr,int x){
	update(no);
	if (l>qr || r<ql)return;
	if (l>=ql && qr>=r){
		sign[no]+=x;
		update(no);
		return;
	}
	int mid=(l+r)>>1;
	add(2*no,l,mid,ql,qr,x);
	add(2*no+1,mid+1,r,ql,qr,x);
	s[no]=min(s[2*no],s[2*no+1]);	
}
void build(int no,int l,int r){
	if (l==r)s[no]=a[l];
	else{
		int mid=(l+r)>>1;
		if (l<=mid)build(2*no,l,mid);
		if (mid<r)build(2*no+1,mid+1,r);
		s[no]=min(s[2*no],s[2*no+1]);
	}
}
int main(){
	freopen("classroom.in","r",stdin);
	freopen("classroom.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	build(1,1,n);
	for (int i=1;i<=m;i++)
		scanf("%d%d%d",&d[i],&b[i],&t[i]);
	for (int i=1;i<=m;i++){
		add(1,1,n,b[i],t[i],d[i]);
		if (s[1]<0){
			printf("-1\n%d",i);
			return 0;
		}
	}
	printf("0");
	printf("0");
	return 0;
}


疫情控制:

又是LCA,看来NOIP是把倍增思想翻来覆去的考。

二分结果,将所有的军队都向上走,如果能到根节点那是更好,如果不能就驻守在那个点,打上标记,然后统计根节点的哪些儿子没有驻守,用贪心策略按照剩下时间的多少派剩下的军队去对应所需时间多少的儿子结点,可以分配则缩小区间,反之同理。


参考程序:

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=211000;
typedef long long LL;
struct Node{
	int j,v,next;
}e[maxn];
struct NODE{
	int n1;
	LL n2;
	bool operator < (const NODE &b)const{
		return n2<b.n2;
	}
}b[maxn],c[maxn];
bool bo[maxn];
int a[maxn];
int fa[maxn][25];
LL len[maxn][25];
int id[maxn];
int n,m,cnt1,cnt2,NodeCnt=0;
void addedge(int u,int v,int w){
	int p=++NodeCnt;
	e[p].j=v;e[p].v=w;e[p].next=a[u];
	a[u]=p;
}
void dfs(int u){
	for (int p=a[u];p;p=e[p].next){
		int j=e[p].j;
		if (j!=fa[u][0]){
			fa[j][0]=u;
			len[j][0]=e[p].v;
			dfs(j);
		}
	}
}
void dfs2(int u){
	if (bo[u])return;
	bool done=true,first=true;
	for (int p=a[u];p;p=e[p].next){
		int j=e[p].j;
		if (j!=fa[u][0]){
			dfs2(j);
			done=false;
			if (!bo[j])first=false;
		}
	}
	if (u!=1 && first && !done)bo[u]=true;
}
bool check(long long time){
	memset(bo,0,sizeof(bo));cnt1=cnt2=0;
	memset(c,0,sizeof(c));
	memset(b,0,sizeof(b));
	for (int i=1;i<=m;i++){
		int up=id[i];
		LL tt=time;
		for (int j=19;j>=0;j--)
			if (fa[up][j] && tt>=len[up][j]){
				tt-=len[up][j];
				up=fa[up][j];
			}
		if (up!=1)bo[up]=true;
			else {
				c[++cnt1].n2=tt;
				up=id[i];
				for (int j=19;j>=0;j--)
					if (fa[up][j]>1)up=fa[up][j];
				c[cnt1].n1=up;
			}
	}
	dfs2(1);
	for (int p=a[1];p;p=e[p].next){
		int j=e[p].j;
		if (!bo[j])b[++cnt2].n1=j,b[cnt2].n2=len[j][0];
	}
	if (cnt1<cnt2)return false;
	sort(c+1,c+1+cnt1);
	sort(b+1,b+1+cnt2);
	int j=1;
	b[cnt2+1].n2=0x7fffffff;
	for (int i=1;i<=cnt1;i++){
		if (!bo[c[i].n1])bo[c[i].n1]=true;
		else{
			if (c[i].n2>=b[j].n2)bo[b[j].n1]=true,j++;
		}
		while (bo[b[j].n1])j++;
	}
	return j>cnt2;
}
int main(){
	freopen("blockade.in","r",stdin);
	freopen("blockade.out","w",stdout);
	scanf("%d",&n);
	LL r=0;
	for (int i=1;i<n;i++){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		r+=z;
		addedge(x,y,z);
		addedge(y,x,z);
	}
	dfs(1);
	for (int j=1;j<20;j++)
	    for (int i=1;i<=n;i++)
			fa[i][j]=fa[fa[i][j-1]][j-1],len[i][j]=len[fa[i][j-1]][j-1]+len[i][j-1];
	scanf("%d",&m);
	for (int i=1;i<=m;i++)scanf("%d",&id[i]);
	LL l=0;
	while (l<r){
		LL mid=(l+r)>>1;
		if (check(mid))r=mid;
			else l=mid+1;
	}
	printf("%lld",l);
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值