NOIP 2012普及组组题解

最新推荐文章于 2024-05-08 22:04:44 发布

_ducati

最新推荐文章于 2024-05-08 22:04:44 发布

阅读量259

点赞数 1

文章标签：算法

本文链接：https://blog.csdn.net/Cherrt/article/details/105440005

版权

话说我为什么要写NOIP普及组的题解?

毕竟~~全做出来了~~ ，要纪念一下。

D1T1 质因数分解

思路:
直接模拟。注意枚举小的质因数，而不是大的质因数，否则会超时。

评价:
①难度: 入门
②做题时间: 1min
③得分: 100

代码:

#include <bits/stdc++.h>
using namespace std;

bool moid(int n)
{
	if (n==0||n==1)  return false;
	for (int i=2;i*i<=n;i++)
	{
		if (n%i==0)  return false;
	}
	return true;
}

int main()
{
	int n;
	cin>>n;
	for (int i=2;i*i<=n;i++)
	{
		if (n%i==0&&moid(i))
		{
			int k=n/i;
		    cout<<max(k,i)<<endl;
		    return 0;
		}
	}
	return 0;
}

D1T2 寻宝

思路:
直接模拟即可。调试时间有些长，但是心态绝对不能崩!

评价:
①难度: 普及-
②时间: 42min(话说我做一道普及-的题目做个42分钟是不是太弱了)
③得分: 100

代码: 略(由于调试过度，代码丑的要命)

D1T3 摆花

思路:

首先，我翻译了一下题目；他就是需要你构造一个由n个数组成的序列，使得第i个数不大于ai；并且这n个数的和需要为m。求出这样的不同序列的个数，并将答案对1000007取模。

接着，通过我的经验(普及组通常，第一题水题，第二题模拟，第三题动规，第四题普通数据结构(什么高级数据结构，如平衡树等考不到))，得到这是一道dp题。

然后开始推公式，~~一下就推出来了~~ 。

状态设计为dp[i][j]表示目前看到第i个数，并且前i个数的和为j的方案数。

状态转移比较难表达(看来我还是太弱了)。事实上就是这样的: dp[i][j]可以为下一行(dp[i+1])作出贡献，即dp[i+1][j]至dp[i][min(m,j+a[i+1])]都可以通过这里的dp[i][j]得到。简略地证明一下，这里得到的和与下一个数能够拼成新的和，且下一个数一定在0至a[i+1]之间，而且和不能超过m，就表达为dp[i+1][j]至dp[i][min(m,j+a[i+1])]均加上dp[i][j]。

但是这种思路的时间复杂度为O( $nm^2$ )，会超时。

接着发现，上面说到了区间加法(dp[i+1][j]至dp[i][min(m,j+a[i+1])]均加上dp[i][j)，那么我们就可以用差分来加速这一过程(树状数组与线段树也可以，但是绝对不推荐，因为这样会加上调试时间)。

啥也不用说了，写完程序测完样例对拍完毕，提交!

结果……WA了! o(╥﹏╥)o qaq

然后我仔细地检查了代码，发现并无错误；然后又看了题面，还是觉得自己并没有理解错。于是，在我放弃准备点开题解的一瞬间，我突然发现:

请输出方案数对1000007取模的结果。

本蒟蒻竟然看成了 $10^9+7$ !

于是，我把mod给改成了1000007，最后就AC了。

所以说，大家一定要看清楚模数啊，不然到时候后悔可TM来不及啊! 要记住到时候普及组只有一次提交机会啦!

评价:
①难度: 普及/提高-
②时间: 17(看题想思+编码)+13(颓废,以为模数是1e9+7)+1(AC)=31

代码:

#include <bits/stdc++.h>
#define int long long
using namespace std;

int n,m,mod=1000007;
int a[105],dp[105][10005];

inline void get_sum(int ii,int l,int r,int c)
{
	dp[ii][l]=(dp[ii][l]+c)%mod;
	dp[ii][r+1]=((dp[ii][r+1]-c)%mod+mod)%mod;
}

signed main()
{
	cin>>n>>m;
	for (int i=1;i<=n;i++)  cin>>a[i];
	for (int i=0;i<=a[1];i++)  dp[1][i]=1;
	
	for (int i=1;i<=n-1;i++)
	{
		int sumv=0;
		for (int j=0;j<=m;j++)  get_sum(i+1,j,min(j+a[i+1],m),dp[i][j]);
		for (int j=0;j<=m;j++)
		{
			sumv=(sumv+dp[i+1][j])%mod;
			dp[i+1][j]=sumv;
		}
	}
	cout<<dp[n][m]<<endl;
	
	return 0;
}

D1T4 文化之旅

思路:
这是一道错题。但是玄学做法有很多，这里我采用的做法是Floyd求最短路+简单模拟。我们在松弛的时候，仅仅需要判断用于松弛的那个点是否与另两个点的文化冲突或相同；如果是就不能松弛，否则就直接松弛。

当然这种做法还是有错误的(从一个点到另一点，不能仅仅判断这两个点的文化是否冲突或相同来决定能否互相到达，而要判断这条路之间所有的国家的文化是否存在冲突或系相同)。~~但是谁叫数据水呢~~

后来我为了巩固Dijkstra就又写了一个Dijkstra的版本(不喜欢SPFA, 现在出题人的一大爱好就是卡SPFA)。那么，两个代码就都放吧~

评价:
①难度: 普及+/提高
②时间: 20(Floyd)+21(Dijkstra)=41min
③得分: 100

代码:
①Dijkstra

#include <bits/stdc++.h>
#define inf 2000000000
#define int long long
using namespace std;

int n,k,m,s,t,cnt=0;
int li[105],head[105],dis[105],visited[105]={0};

int ag[105][105],GA[105][105];

struct edge
{
	int next;
	int to;
	int dis;
}e[10005];

struct node
{
	int dis;
	int pos;
	bool operator < (const node &x) const
	{
		return x.dis<dis;
	}
};

std::priority_queue<node> q;

inline void add_edge(int u,int v,int w)
{
	cnt++;
	e[cnt].to=v;
	e[cnt].dis=w;
	e[cnt].next=head[u];
	head[u]=cnt;
}

inline void dijkstra()
{
	dis[s]=0;
	q.push((node){0,s});
	while (!q.empty())
	{
		node tmp=q.top();
		q.pop();
		int x=tmp.pos;
		if (visited[x])  continue;
		visited[x]=1;
		for (int i=head[x];i;i=e[i].next)
		{
			int y=e[i].to;//i优化sy
			int xx=li[s],yy=li[y];
			if (xx==yy||(ag[xx][yy]==1&&ag[yy][xx]==1))  continue;
			if (dis[y]>dis[x]+e[i].dis)
			{
				dis[y]=dis[x]+e[i].dis;
				if (!visited[y])  q.push((node){dis[y],y});
			}
		}
	}
}

inline void init()
{
	cin>>n>>k>>m>>s>>t;
	for (int i=1;i<=n;i++)  cin>>li[i];
	for (int i=1;i<=k;i++)
	{
		for (int j=1;j<=k;j++)  cin>>ag[i][j];
	}
	for (int i=1;i<=m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		if (li[u]==li[v])  continue;
		if (ag[li[u]][li[v]]==0&&ag[li[v]][li[u]]==0)
		{
			add_edge(u,v,w);
			add_edge(v,u,w);
		}
		else if (ag[li[u]][li[v]]==0&&ag[li[v]][li[u]]==1)  add_edge(v,u,w);
		else if (ag[li[u]][li[v]]==1&&ag[li[v]][li[u]]==0)  add_edge(u,v,w);
		else continue;
	}
}

signed main()
{
	init();
	
	for (int i=1;i<=n;i++)  dis[i]=inf;
	
	dijkstra();
	
	if (dis[t]!=inf)  cout<<dis[t]<<endl;
	else cout<<-1<<endl;
	
	return 0;	
}

②Floyd

#include <bits/stdc++.h>
#define inf 2000000000
#define int long long
using namespace std;

int n,k,m,s,t,li[105];
int ag[105][105],GA[105][105];

signed main()
{
	cin>>n>>k>>m>>s>>t;
	for (int i=1;i<=n;i++)  cin>>li[i];
	for (int i=1;i<=k;i++)
	{
		for (int j=1;j<=k;j++)  cin>>ag[i][j];
	}
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=n;j++)
		{
			if (i!=j)  GA[i][j]=inf;
			else GA[i][j]=0;
		}
	}
	for (int i=1;i<=m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		if (li[u]==li[v])  continue;
		if (ag[li[u]][li[v]]==0&&ag[li[v]][li[u]]==0)  GA[u][v]=GA[v][u]=w;
		else if (ag[li[u]][li[v]]==0&&ag[li[v]][li[u]]==1)  GA[v][u]=w;
		else if (ag[li[u]][li[v]]==1&&ag[li[v]][li[u]]==0)  GA[u][v]=w;
		else continue;
	}
    
    if (li[s]==li[t]||(ag[t][s]==1))
    {
    	cout<<-1<<endl;
    	return 0;
	}
    
	for (int k=1;k<=n;k++)
	{
		for (int i=1;i<=n;i++)
		{
			for (int j=1;j<=n;j++)
			{
				int x=li[i],y=li[j];
				if (ag[x][y]==1||x==y)  continue;
				if (GA[i][k]+GA[k][j]<GA[i][j])
				{
					GA[i][j]=GA[i][k]+GA[k][j];
					GA[j][i]=GA[i][k]+GA[k][j];
				}
			}
		}
	}
	if (GA[s][t]!=inf)  cout<<GA[s][t]<<endl;
	else cout<<-1<<endl;
	return 0;
}

_ducati

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
NOIP 2012普及组组题解

话说我为什么要写NOIP普及组的题解?毕竟全做出来了，要纪念一下。D1T1 质因数分解思路:直接模拟。注意枚举小的质因数，而不是大的质因数，否则会超时。评价:①难度: 入门②做题时间: 1min③得分: 100代码:#include <bits/stdc++.h>using namespace std;bool moid(int n){ if (n==0...
复制链接

扫一扫