话说我为什么要写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( n m 2 nm^2 nm2),会超时。
接着发现,上面说到了区间加法(dp[i+1][j]至dp[i][min(m,j+a[i+1])]均加上dp[i][j),那么我们就可以用差分来加速这一过程(树状数组与线段树也可以,但是绝对不推荐,因为这样会加上调试时间)。
啥也不用说了,写完程序测完样例对拍完毕,提交!
结果……WA了! o(╥﹏╥)o qaq
然后我仔细地检查了代码,发现并无错误;然后又看了题面,还是觉得自己并没有理解错。于是,在我放弃准备点开题解的一瞬间,我突然发现:
请输出方案数对1000007取模的结果。
本蒟蒻竟然看成了 1 0 9 + 7 10^9+7 109+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;
}