noip2017提高组复赛解题报告

noip2017提高组复赛解题报告定期推送帐号信息学新闻，竞赛自主招生，信息学专业知识，信息学疑难解答，融科教育信息学竞赛培训等诸多优质内容的微信平台，欢迎分享文章给你的朋友或者朋友圈！以下解题思路及代码未经官方评测，仅供参考，复赛成绩以官方（CCF）评测结果为准。

Day11.小凯的疑惑(math.cpp/c/pas)【问题描述】小凯手中有两种面值的金币，两种面值均为正整数且彼此互素。

每种金币小凯都有无数个。

在不找零的情况下，仅凭这两种金币，有些物品他是无法准确支付的。

现在小凯想知道在无法准确支付的物品中，最贵的价值是多少金币?注意:输入数据保证存在小凯无法准确支付的商品。

【输入格式】输入文件名为math.in。

输入数据仅一行，包含两个正整数a 和b，它们之间用一个空格隔开，表示小凯手中金币的面值。

【输出格式】输出文件名为math.out。

输出文件仅一行，一个正整数N，表示不找零的情况下，小凯用手中的金币不能准确支付的最贵的物品的价值。

【输入输出样例1】math.in3 7 math.out11【数据规模与约定】对于30%的数据: 1 ≤a，b ≤50。

对于60%的数据: 1 ≤a，b ≤10,000。

对于100%的数据:1 ≤a，b ≤1,000,000,000。

数学太差只找规律吧。

设：其中一个数为2则：2、3=>1；2、5=>3；2、7=>5；2、11=>9得：2、n=>n-2设：其中一个数为3则：3、5=>7；3、7=>11；3、11=>19；3、13=>23得：3、n=>2n-3设：其中一个数为5则：5、7=>23；5、11=>39；5、13=>47；5、17=>63得：5、n=>4n-5所以：m、n=>(m-1)n-m #includeusing namespace std;int main(){ long long a,m,n;scanf('%lld %lld',&m,&n); a=(m-1)*n-m;printf('%lld',a); return 0;} 2.时间复杂度(complexity.cpp/c/pas)【问题描述】小明正在学习一种新的编程语言A++，刚学会循环语句的他激动地写了好多程序并给出了他自己算出的时间复杂度，可他的编程老师实在不想一个一个检查小明的程序，于是你的机会来啦!下面请你编写程序来判断小明对他的每个程序给出的时间复杂度是否正确。

A++语言的循环结构如下:其中“F i x y”表示新建变量(i 变量i 不可与未被销毁的变量重名)并初始化为x，然后判断i 和y 的大小关系，若i 小于等于y 则进入循环，否则不进入。

每次循环结束后i都会被修改成i +1，一旦i 大于y 终止循环。

x和y 可以是正整数(x 和y 的大小关系不定)或变量n。

n 是一个表示数据规模的变量，在时间复杂度计算中需保留该变量而不能将其视为常数，该数远大于100。

“E”表示循环体结束。

循环体结束时，这个循环体新建的变量也被销毁。

注:本题中为了书写方便，在描述复杂度时，使用大写英文字母“O”表示通常意义下“Θ”的概念。

【输入格式】输入文件名为complexity.in。

输入文件第一行一个正整数t，表示有t(t ≤10)个程序需要计算时间复杂度。

每个程序我们只需抽取其中“F i x y”和“E”即可计算时间复杂度。

注意:循环结构允许嵌套。

接下来每个程序的第一行包含一个正整数L 和一个字符串，L 代表程序行数，字符串表示这个程序的复杂度，“O(1)”表示常数复杂度，“O(n^w)”表示复杂度为n^w，其中w 是一个小于100 的正整数(输入中不包含引号)，输入保证复杂度只有O(1)和O(n^w) 两种类型。

接下来L 行代表程序中循环结构中的“F i x y”或者“E”。

程序行若以“F”开头，表示进入一个循环，之后有空格分离的三个字符(串)i x y，其中i 是一个小写字母(保证不为“n”)，表示新建的变量名，x 和y 可能是正整数或n ，已知若为正整数则一定小于100。

程序行若以“E”开头，则表示循环体结束。

【输出格式】输出文件名为complexity.out。

输出文件共t 行，对应输入的t 个程序，每行输出“Yes”或“No”或者“ERR”(输出中不包含引号)，若程序实际复杂度与输入给出的复杂度一致则输出“Yes”，不一致则输出“No”，若程序有语法错误(其中语法错误只有: F 和E 不匹配；新建的变量与已经存在但未被销毁的变量重复两种情况)，则输出“ERR”。

注意:即使在程序不会执行的循环体中出现了语法错误也会编译错误，要输出“ERR”。

【输入输出样例1】complexity.in82 O(1)F i 1 1E2O(n^1)F x 1 nE1 O(1)F x 1 n4 O(n^2)F x 5 nFy 10 nEE4 O(n^2)F x 9 nEF y 2 nE4 O(n^1)F x 9 nF y n 4EE4 O(1)F y n 4F x 9 nEE4 O(n^2)F x 1 nF x 110EEcomplexity.outYesYesERRYesNoYesYesERR【数据规模与约定】对于30%的数据:不存在语法错误，数据保证小明给出的每个程序的前L/2 行一定为以 F 开头的语句，第L/2+1 行至第L 行一定为以E 开头的语句，L，若x、y 均为整数，x 一定小于y，且只有y 有可能为n。

对于50%的数据:不存在语法错误，L，且若x、y均为整数，x 一定小于y，且只有y 有可能为n。

对于70%的数据:不存在语法错误，L。

对于100%的数据:L。

用STL stack模拟，对于用代码量堆出来的OIer来说，这就是信心倍增的大力模拟题。

代码就不上了。

3. 逛公园(park.cpp/c/pas)【问题描述】策策同学特别喜欢逛公园。

公园可以看成一张N个点M条边构成的有向图，且没有自环和重边。

其中1号点是公园的入口，N号点是公园的出口，每条边有一个非负权值，代表策策经过这条边所要花的时间。

策策每天都会去逛公园，他总是从1号点进去，从N号点出来。

策策喜欢新鲜的事物，他不希望有两天逛公园的路线完全一样，同时策策还是一个特别热爱学习的好孩子，他不希望每天在逛公园这件事上花费太多的时间。

如果1号点到N号点的最短路长为d，那么策策只会喜欢长度不超过d+K的路线。

策策同学想知道总共有多少条满足条件的路线，你能帮帮他吗?为避免输出过大，答案对P取模。

如果有无穷多条合法的路线，请输出?1。

【输入格式】输入文件名为park.in。

第一行包含一个整数T, 代表数据组数。

接下来T组数据，对于每组数据:第一行包含四个整数N, M, K, P，每两个整数之间用一个空格隔开。

接下来M行，每行三个整数ai , bi , ci，代表编号为ai, bi的点之间有一条权值为ci的有向边，每两个整数之间用一个空格隔开。

【输出格式】输出文件名为park.out。

输出文件包含T行，每行一个整数代表答案。

【输入输出样例1】park.in25 7 2 101 2 12 4 04 5 22 3 23 4 13 5 21 5 32 2 0 101 2 02 1 0park.out3-1【数据规模与约定】对于不同的测试点，我们约定各种参数的规模不会超过如下：对于100%的数据, 1≤P≤109,1≤ai,bi≤??,0≤ci ≤1000。

数据保证:至少存在一条合法的路线。

最短路+拓扑排序+DP先跑最短路。

发现K只有50，所以一定是要从K入手。

所以考虑DP，令f[i][j]表示走到i，多走的长度是j的方案数。

（多走指的是比最短路多的部分的长度）。

但是发现这个DP方程是存在环的，因为最短路径图上的边以及零边都是可以同行转移的。

将最短路径图上的边以及零边都拿出来跑拓扑排序，然后让这些边在转移时必须沿着拓扑序转移即可。

特别地，如果一个零环位于一条从1到n长度的路径上，则输出-1即可。

#include#include#include#include#include#include#define mp(A,B) make_pair(A,B)using namespace std;const intmaxn=100010;const int maxm=200010;intn,m,cnt,K,P,tot,ans;intto[maxm],next[maxm],val[maxm],head[maxn],dis[maxn],vis[m axn],d[maxn],q[maxnintto2[maxm],next2[maxm],val2[maxm],head2[maxn],dis2[maxn]; priority_queue> pq;inline int rd(){ int ret=0; chargc=getchar(); while(gc'9') gc=getchar();while(gc>='0'&&gc return ret;}inline void add(int a,int b,int c){ to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt;to2[cnt]=a,val2[cnt]=c,next2[cnt]=head2[b],head2[b]=cnt++;}in line void upd(int &x,int y){ x=(x+y)%P;}voidwork(){ n=rd(),m=rd(),K=rd(),P=rd(); int i,j,k,a,b,c,u; memset(head,-1,sizeof(head)),memset(head2,-1,sizeof(head2)),c nt=tot=ans=0; for(i=1;i //1memset(vis,0,sizeof(vis)),memset(dis,0x3f,sizeof(dis)),memset( d,0,sizeof(d)); pq.push(mp(0,1)),dis[1]=0;while(!pq.empty()) { u=pq.top().second,pq.pop(); if(vis[u]) continue; vis[u]=1;for(i=head[u];i!=-1;i=next[i])if(dis[to[i]]>dis[u]+val[i])dis[to[i]]=dis[u]+val[i],pq.push(mp(-d is[to[i]],to[i])); } //2memset(dis2,0x3f,sizeof(dis2)),memset(vis,0,sizeof(vis));pq.push(mp(0,n)),dis2[n]=0; while(!pq.empty()){ u=pq.top().second,pq.pop(); if(vis[u]) continue; vis[u]=1;for(i=head2[u];i!=-1;i=next2[i])if(dis2[to2[i]]>dis2[u]+val2[i])dis2[to2[i]]=dis2[u]+val2[i],pq.p ush(mp(-dis2[to2[i]],to2[i])); } //3 for(i=1;ifor(i=1;i for(j=1;j { u=q[j];for(i=head[u];i!=-1;i=next[i]) if(dis[u]+val[i]==dis[to[i]]){ d[to[i]]--; if(!d[to[i]])q[++tot]=to[i]; } } for(i=1;i{ printf('-1\n'); return ; } //DP memset(f,0,sizeof(f)); f[0][1]=1; for(k=0;k{ for(i=1;i{ upd(f[k][to[j]],f[k][u]); }for(i=1;i{ upd(f[k+dis[i]+val[j]-dis[to[j]]][to[j]],f[k][i]);} } for(i=0;i printf('%d\n',ans);}intmain(){ int T=rd(); while(T--) work(); return 0;} Day2 1.奶酪(cheese.cpp/c/pas)【问题描述】现有一块大奶酪，它的高度为h，它的长度和宽度我们可以认为是无限大的，奶酪中间有许多半径相同的球形空洞。