【比赛报告】NOIP2016真题 NOIP练习赛卷三十一

玩具谜题 模拟

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


按照题意模拟就好

#include<cstdio>
struct node{
	int fa;
	char name[20];
}q[100100];
int main(){
	freopen("testdata.in","r",stdin);
	int n,m,i,f,num,h=0;
	scanf("%d%d",&n,&m);
	for(i=0;i<n;i++)scanf("%d%s",&q[i].fa,q[i].name);
	//i从1开始 PAC 
	for(i=0;i<m;i++){
		scanf("%d%d",&f,&num);
		if(q[h].fa==f)h=(h-num+n)%n;
		else h=(h+num)%n;
	}
	printf("%s",q[h].name);
	return 0;
}

总结

简单模拟


天天爱跑步 LCA

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


学习了大佬题解,讲的非常清晰。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
#define PB(v) push_back(v)
template<class T>inline void read(T&x)
{
	x=0;int f=0;char ch=getchar();
	while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
	while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	if(f)x=-x;
}
const int N=3e5+10;
const int M=6e5+10;
int n,m,hd[N],tot,w[N],bucket[N+M],ans[N],deep[N],f[N][25],dis[N],spn[M];
bool vis[N];
vector<int>v1[M],v2[M],v3[M];
struct Edge{
	int v,nx;
}e[M];
struct node{
	int u,v,dis,lca;
}p[N];
inline void add(int u,int v)
{
	e[++tot].v=v;
	e[tot].nx=hd[u];
	hd[u]=tot;
}
inline void dfs(int now,int dep)
{
	vis[now]=true;
	deep[now]=dep;
	for(int i=1;i<=22;i++)
	    f[now][i]=f[f[now][i-1]][i-1];
	for(int i=hd[now];i;i=e[i].nx)
	    if(!vis[e[i].v])
	    {
	    	f[e[i].v][0]=now;
	    	dis[e[i].v]=dis[now]+1;
	    	dfs(e[i].v,dep+1);
		}
	vis[now]=false;
}
inline int LCA(int u,int v)
{
	if(deep[u]<deep[v])swap(u,v);
	for(int i=0,div=deep[u]-deep[v];i<=22;i++)
	    if(div&(1<<i))u=f[u][i];
	for(int i=22;i>=0;i--)
	    if(f[u][i]!=f[v][i])u=f[u][i],v=f[v][i];
	return u==v?u:f[u][0];
}
//u->lca(u,v)的路径,当且仅当deep[i]+w[i]==deep[u]时,u节点对i节点是有贡献的 
//只要符合deep[i]+w[i]的全部是玩家起点的点,就能对i点产生贡献 
inline void dfs1(int now)
{
	vis[now]=true;
	int prev=bucket[deep[now]+w[now]+N];//全部+N避免负数 
	for(int i=hd[now];i;i=e[i].nx)
	    if(!vis[e[i].v])dfs1(e[i].v);
	bucket[deep[now]+N]+=spn[now];
	ans[now]+=bucket[deep[now]+w[now]+N]-prev;
	int len=v1[now].size();
	for(int k=0;k<len;k++)
	    --bucket[deep[v1[now][k]]+N];
	vis[now]=false;
}
//lca(u,v)->v的路径,当且仅当dis(u,v)-deep[v]=w[i]-deep[i]成立时有贡献 
inline void dfs2(int now)
{
	vis[now]=true;
	int prev=bucket[w[now]-deep[now]+N];
	for(int i=hd[now];i;i=e[i].nx)
	    if(!vis[e[i].v])dfs2(e[i].v);
	int len=v2[now].size();
	for(int k=0;k<len;k++)
	    ++bucket[v2[now][k]+N];
	ans[now]+=bucket[w[now]-deep[now]+N]-prev;
	len=v3[now].size();
	for(int k=0;k<len;k++)
	    --bucket[v3[now][k]+N];
	vis[now]=false;
}
int main()
{
	//freopen("in.txt","r",stdin);
	read(n);read(m);
	for(int i=1,u,v;i<n;i++)
	    read(u),read(v),add(u,v),add(v,u);
    for(int i=1;i<=n;i++)
        read(w[i]);
    f[1][0]=1;dfs(1,0);
    for(int i=1,u,v;i<=m;i++)
    {
    	read(u);read(v);
    	p[i].u=u;p[i].v=v;
    	p[i].lca=LCA(u,v);
    	p[i].dis=dis[u]+dis[v]-dis[p[i].lca]*2;
    	spn[u]++;//以x为路径起点的路径条数 
    	v1[p[i].lca].PB(u);//以x为lca的路径的起点的集合 
    	v2[v].PB(p[i].dis-deep[p[i].v]);//以x为终点的路径的起点集合 
    	v3[p[i].lca].PB(p[i].dis-deep[p[i].v]);//以x为LCA的路径的终点的集合 
	}
	dfs1(1);//从下至上
	dfs2(1);//从上至下 
	for(int i=1;i<=m;i++)//如果lca能观察到这条路径上的人,需要去重 
	    if(deep[p[i].u]==deep[p[i].lca]+w[p[i].lca])ans[p[i].lca]--;
	for(int i=1;i<=n;i++)
	    printf("%d ",ans[i]);
	printf("\n");
	return 0;
}

总结

难度挺大的题


换教室 期望DP+最短路

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


好久以前做的了

#include<cstdio>
#include<cstring>
#include<algorithm>
#define INF 0x3f3f3f3f
using namespace std;
const int N1=2010;
const int N2=310;
double dp[N1][N1][2];
//dp[i][j][0]表示前i个时间段用j次申请(第i个不申请)所需要的最小体力值
//dp[i][j][1]表示申请
double k[N1];//通过概率
int c[N1];//安排教室 
int d[N1];//另外教室 
int dis[N1][N1];//边权 
int main()
{
	//freopen("in.txt","r",stdin);
	int n,m,v,e,a,b;
	scanf("%d%d%d%d",&n,&m,&v,&e);//时间,申请,教室,道路 
	int i,j,l,w;//边权 
	for(i=1;i<=n;i++)
	scanf("%d",&c[i]);//安排 
	for(i=1;i<=n;i++)
	scanf("%d",&d[i]);//另外
	for(i=1;i<=n;i++)
	scanf("%lf",&k[i]);//概率
	memset(dis,127/2,sizeof(dis));
	for(i=1;i<=e;i++)
	{
		scanf("%d%d%d",&a,&b,&w);
		dis[a][b]=min(dis[a][b],w);
		dis[b][a]=min(dis[b][a],w);
	}
	for(i=1;i<=v;i++)dis[i][i]=0;
	for(l=1;l<=v;l++)
	for(i=1;i<=v;i++)
	for(j=1;j<=v;j++)
	if(dis[i][j]>dis[i][l]+dis[l][j])//处理最短路 
	dis[i][j]=dis[i][l]+dis[l][j];
	for(i=1;i<=v;i++)dis[0][i]=0;
	for(i=0;i<=2009;i++)
	for(j=0;j<=2009;j++)
	for(l=0;l<=1;l++)
	dp[i][j][l]=INF;
	//memset(dp,0x3f,sizeof(dp));
	dp[0][0][0]=0;
	for(i=0;i<=n-1;i++)
	for(j=0;j<=m;j++)
	{
		dp[i+1][j][0]=min(dp[i][j][0]+dis[c[i]][c[i+1]],dp[i][j][1]+dis[d[i]][c[i+1]]*k[i]+dis[c[i]][c[i+1]]*(1.0-k[i]));
		if(j!=m)dp[i+1][j+1][1]=min(dp[i][j][0]+dis[c[i]][d[i+1]]*k[i+1]+dis[c[i]][c[i+1]]*(1.0-k[i+1]),dp[i][j][1]+dis[c[i]][c[i+1]]*(1.0-k[i+1])*(1.0-k[i])+dis[d[i]][c[i+1]]*(1.0-k[i+1])*k[i]+dis[c[i]][d[i+1]]*k[i+1]*(1.0-k[i])+dis[d[i]][d[i+1]]*k[i]*k[i+1]);
	}
	double ans=INF;
	for(j=0;j<=m;j++)
	{
		ans=min(ans,min(dp[n][j][0],dp[n][j][1]));
	}
	printf("%.2lf\n",ans);
	return 0;
}

总结


组合数问题 组合计数+前缀和

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


线性递推组合数,然后求一个前缀和。
对着大佬题解改改改一直没对,然后发现m可以大于n……

#include<cstdio>
const int N=2e3+10;
int t,k,n,m;
long long C[N][N],sum[N][N];
int main()
{
    //freopen("in.txt","r",stdin);
    scanf("%d%d",&t,&k);
    for(int i=0;i<=2000;i++)C[i][0]=1;
    C[1][1]=1;
    for(int i=2;i<=2000;i++)
        for(int j=1;j<=i;j++)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%k;
    for(int i=0;i<=2000;i++)
    {
        for(int j=1;j<=i;j++)
        {
        	sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
        	if(C[i][j]==0)sum[i][j]++;
        }
        sum[i][i+1]=sum[i][i];
    }
    while(t--)
    {
        scanf("%d%d",&n,&m);if(m>n)m=n;
        printf("%lld\n",sum[n][m]);
    }
    return 0;
}

总结

考察组合数的递推


蚯蚓 优先队列+模拟

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


最开始敲了个85的模拟……想想也没那么简单。
学习了大佬题解,其中单调性很关键。

#include<cstdio>
#include<queue>
#include<cmath>
#include<algorithm>
#define re register 
using namespace std;
template<class T>inline void read(T&x)
{
	x=0;re int f=0;re char ch=getchar();
	while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
	while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	if(f)x=-x;
}
template<class M>inline void write(M x)
{
	re int buf[40],p=0;
	if(x<0)putchar('-'),x=-x;
	do{
		buf[p++]=x%10;x/=10;
	}while(x);
	for(re int i=p-1;i>=0;i--)putchar(buf[i]+48);
	putchar(' ');
}
const int N=7e6+10;
int n,m,q,u,v,t,sigma,t0,t1,t2,h0,h1,h2,cut1[N],cut2[N],now[N];
double p;
priority_queue<int>ans;
inline bool cmp(const int &a,const int &b){return a>b;}
int main()
{
	//freopen("in.txt","r",stdin);
    read(n);read(m);read(q);read(u);read(v);read(t);
    p=(double)u/v;
    for(re int i=1;i<=n;i++)
        read(now[i]);
    t0=n;t1=t2=0;h0=h1=h2=1;
    sort(now+1,now+n+1,cmp);
    re int top;
    for(re int i=1,a1,a2;i<=m;i++)
    {
    	if(h0>t0)
    	{
    		
    		if(cut1[h1]>cut2[h2])top=cut1[h1++];
    		else top=cut2[h2++];
		}
		else if(now[h0]>=cut1[h1]&&now[h0]>=cut2[h2])top=now[h0++];
		else if(now[h0]<=cut1[h1]&&cut2[h2]<=cut1[h1])top=cut1[h1++];
		else top=cut2[h2++];
		top+=sigma;
		a1=floor(p*(double)top);a2=top-a1;
		a1-=sigma+q;a2-=sigma+q;
		cut1[++t1]=a1;cut2[++t2]=a2;
		if(i%t==0)write(top);sigma+=q;
	}
	putchar('\n');
	for(re int i=h0;i<=t0;i++)ans.push(now[i]);
	for(re int i=h1;i<=t1;i++)ans.push(cut1[i]);
	for(re int i=h2;i<=t2;i++)ans.push(cut2[i]);
	for(re int i=1;ans.size();i++)
	{
		if(i%t==0)write(ans.top()+sigma);
		ans.pop();
	}
	return 0;
}

总结

挺好的一个大模拟。


愤怒的小鸟 状压DP

题目链接
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


dp[i]表示i状态时所需要的最少的小鸟数,state[i]表示第i条抛物线所打掉的小猪状态,
dp[i|state[j]]=min(dp[i|state[j]],dp[i]+1)

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define eps 1e-6
int n,m;
const int N=18;
int dp[1<<N];//dp[i]表示i状态时所需最少小鸟数 
double x[N],y[N];
int state[510];//记录每条抛物线 
int cnt; 
int main()
{
	//freopen("in.txt","r",stdin);
	int t;
	double a,b;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&m);
		cnt=0;
		memset(dp,0x3f,sizeof(dp));
		dp[0]=0;
		for(int i=0;i<n;i++)
		scanf("%lf%lf",&x[i],&y[i]);
		for(int i=0;i<n;i++)
		{
			state[cnt++]=(1<<i);
			for(int j=i+1,k=0;j<n;j++)
			{
				if((k>>j)&1)continue;
				a=(y[j]*x[i]-y[i]*x[j])/(x[i]*x[j]*(x[j]-x[i]));
				b=(y[i]*x[j]*x[j]-y[j]*x[i]*x[i])/(x[i]*x[j]*(x[j]-x[i]));
				if(a>=-eps)continue;
				state[cnt]=(1<<i);
				for(int o=j;o<n;o++)
				if(fabs(a*x[o]*x[o]+b*x[o]-y[o])<=eps)
				{
					k|=(1<<o);
					state[cnt]|=(1<<o);
				}
				cnt++;
			}
		}
		for(int i=0;i<(1<<n);i++)
		for(int j=0;j<cnt;j++)
		dp[i|state[j]]=min(dp[i|state[j]],dp[i]+1);
		printf("%d\n",dp[(1<<n)-1]);
	}
	return 0;
} 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值