[考试反思]1010csp-s模拟测试67:摸索

嗯。。。所谓RP守恒?

仍然延续着好一场烂一场的规律。

虽说我也想打破这个规律,但是并不想在考烂之后打破这个规律。(因为下一场要考好???)

我也不知道我现在是什么状态,相较于前一阶段有所提升(第一鸡房的buff?)

但是明显还有提升空间。也可以看到离上面的差距有多大。

(在去食堂的路上一顿爆捶捶傻skyh就能进步一名了hhh)

这场考试呢,T1正路是推式子然而我又开始疯狂打表找规律,最后得到了一个和正解没什么关系但是能拿到80分的式子。

也肝了一个小时,然后就停止在80分了。

并不是没有在化式子,而是化的不彻底,没有到最终可做的形式。

然后T2是个板子就写了。没调没对拍,着急进T3。(又不像某些人考前打对拍。。。)

T3想直接奔着最终答案去,于是就先上了个暴力,对着excel看了挺久,找到一些规律但是不能用于做题。

最后20分钟左右的时候想到了n3的思路,但是细节较多来不及打了,时间分配不够。还是给T3分的时间不够多。

最后交了暴力。最后几分钟打了21~25的表但是没有数据,数据是28。如果提前几分钟开始打就多20分了。。。

然后就没有然后了。

还好,仍然是在持续的思考。没挂机没摸鱼。

但是某些考试技巧还是要练的。。。如打表之类。。。

不要尝试找n2的表的规律,很困难而且浪费时间。

分配一些时间给打表交表也是不错的。

T1:神炎皇

考场上的发现是$- n +\sum \limits_{i=1}^{n} \sqrt{i的最大平方因子}$

枚举这个因子,此后就不应再有平方因子,运用莫比乌斯函数的实际意义进行容斥,得到

$= - n + \sum \limits_{d=1}^{\sqrt{n}} d \sum \limits _{k=1}^{\left \lfloor \frac{n}{d^2} \right \rfloor} \mu (k) \times \left \lfloor \frac{n}{d^2 k^2} \right \rfloor$

这样的式子的复杂度外层是$\sqrt{n}$,内层是调和级数,所以总的复杂度是$O( \sqrt{n} ln \sqrt{n} )$,有80分。

而正解是考虑ab的最大公约数,设$gcd(i,j)=g$,$a=gx$,$b=gy$

那么有$xyg^2 mod (x+y)g =0$。因为xy互质,所以xy与x+y互质

所以$g mod (x+y) =0 $,所以$g>=x+y$。那么根据题目的限制$a+b<=n$则$(x+y)g<=n$

有$x+y<=n$。那么就枚举$k=x+y$。求数对$x,y$的数量。

因为$y$与$x$互质,所以$x+y$与$x$互质,所以确定$k$后可行的$x$就有$\varphi (k)$个。

然后可行的$g$的个数就是$\frac{n}{k^2}$个。

那么最终的答案就是$\sum \limits_{k=1}^{\sqrt{n}} \varphi (k) \times \frac{n}{k^2}$

最终复杂度$ O( \sqrt{n} ) $

 1 #include<cstdio>
 2 int phi[10000005],p[10000005],pcnt;char np[10000005];long long n,ans;
 3 int main(){
 4     for(int i=2;i<=10000000;++i){
 5         if(!np[i])p[++pcnt]=i,phi[i]=i-1;
 6         for(int j=1;j<=pcnt&&i*p[j]<=10000000;++j)
 7             if(i%p[j])phi[i*p[j]]=phi[i]*(p[j]-1),np[i*p[j]]=1;
 8             else {phi[i*p[j]]=phi[i]*p[j],np[i*p[j]]=1;break;}
 9     }
10     scanf("%lld",&n);
11     for(long long i=1;i*i<=n;++i)ans+=n/i/i*phi[i];
12     printf("%lld\n",ans);
13 }

T2:降雷皇

一个比较普通的线段树dp,只不过记录了一下方案。

算是个板子吧。

 1 #include<cstdio>
 2 #define mod 123456789
 3 int Mod(int p){return p>=mod?p-mod:p;}
 4 struct state{
 5     int dp,pl;
 6     friend state operator+(state x,state y){
 7         if(x.dp<y.dp)return y;
 8         if(y.dp<x.dp)return x;
 9         return (state){x.dp,Mod(x.pl+y.pl)};
10     }
11 }w[800005],O;
12 int cl[800005],cr[800005],x[800005];
13 void build(int p,int l,int r){
14     cl[p]=l;cr[p]=r;
15     if(l==r){w[p]=(state){0,l?0:1};return;}
16     build(p<<1,l,l+r>>1);build(p<<1|1,(l+r>>1)+1,r);
17     w[p]=w[p<<1]+w[p<<1|1];
18 }
19 void chg(int p,int pos,state nw){
20     if(cl[p]==cr[p]){w[p]=w[p]+nw;return;}
21     chg(pos<=cl[p]+cr[p]>>1?p<<1:p<<1|1,pos,nw);
22     w[p]=w[p<<1]+w[p<<1|1];
23 }
24 state ask(int p,int l,int r){
25     if(l<=cl[p]&&cr[p]<=r)return w[p];
26     return ask(p<<1,l,r)+(r>=cl[p<<1|1]?ask(p<<1|1,l,r):O);
27 }
28 int main(){
29     int n,opt,x;scanf("%d%d",&n,&opt);
30     build(1,0,100000);
31     while(n--){
32         scanf("%d",&x);
33         state bst=ask(1,0,x-1);bst.dp++;
34         chg(1,x,bst);
35     }
36     state ans=ask(1,0,100000);
37     printf("%d\n",ans.dp);
38     if(opt==1)printf("%d\n",ans.pl);
39 }

T3:幻魔皇

斐波那契树当然离不开斐波那契数列。还记得那个兔子的例子?

白点是小兔子,黑点是老兔子。

在相同时间年龄相同的兔子,其子树完全相同。

既然要求树上的距离,那么肯定是要知道最近公共祖先的。

如果两个白点存在祖先关系那么它们的祖先就是老的那一个白点,否则就是一个黑点。

那么考虑枚举每一个点,考虑它为lca时子数的贡献。

如果是个白点,那么其对距离为i的点对的贡献就是它第i代后代中白点的数量,就是fib[i]。

(具体是fib的第几项可能因人而异,因为我的fib[0]=1,fib[1]=0)

然后枚举每一层的全部白点,它们的数量也就是个fib数,n2解决所有祖先关系的白点对的贡献。

如果lca是黑点,那么黑点有一个黑儿子一个白儿子,lca是它的话那么两个白点要分别在两个子数内。

处理层数为1的,答案为0。

不断考虑层数更高的,其新的贡献就是白儿子的新的一层与原有黑儿子的贡献,以及黑儿子新的一层与新的白儿子的贡献。

具体就是一大堆的fib数。处理出f[i][j]表示深度为i的黑根树对答案j的贡献,乘上整棵树里有多少个这样的黑点即可。

答案数组开10000。而不是5000。(我稍蠢,嗯)

复杂度$O(n^2)$

 1 #include<cstdio>
 2 #define int long long
 3 #define mod 123456789
 4 int fib[5005],ans[10005],f[5005][10005];
 5 main(){
 6     int n;scanf("%lld",&n);
 7     fib[0]=1;
 8     for(int i=2;i<=n;++i)fib[i]=(fib[i-1]+fib[i-2])%mod;
 9     for(int i=0;i<n;++i)for(int j=1;j<n-i;++j)(ans[j]+=fib[i]*fib[j])%=mod;
10     for(int i=2;i<n;++i){
11         for(int j=1;j<=n<<1;++j)f[i][j]=f[i-1][j];
12         for(int j=1;j<=i;++j)(f[i][i+j]+=fib[i-1]*fib[j])%=mod;
13         for(int j=1;j<i;++j)(f[i][i+j]+=fib[i]*fib[j-1])%=mod;
14     }
15     for(int i=1;i<n;++i)for(int j=1;j<=n<<1;++j)(ans[j]+=f[n-i][j]*fib[i])%=mod;
16     for(int i=1;i<=n<<1;++i)printf("%lld ",ans[i]);puts("");
17 }

原文地址:https://www.cnblogs.com/hzoi-DeepinC/p/11648052.html

时间: 2024-11-05 18:57:27

[考试反思]1010csp-s模拟测试67:摸索的相关文章

神炎皇(模拟测试67)(数论)

神炎皇: 题意: 对于一个整数对$(a,b)$,若满足$a+b<=n$且$a+b$是$a*b$的因子,则成为神奇的数对.请问这样的数对共有多少个?($N<=10^{14}$) 题解: 已知$a+b<=n\\ (a+b)|ab$. 设$d=\gcd(a,b),x=a/d,y=b/d$. 上式为$(x+y)*d<=n(1)\\ (x+y)|x*y*d(2)$. 因为$\gcd(x+y,x)=\gcd(x+y,y)=\gcd(x,y)=1$. (2)式可化减为$(x+y)|d$. 又由(

[考试反思]0929csp-s模拟测试55:沦陷

菜得过分. 面对T1的大板子不知所措,然后T2的贪心不小心把排序语句删了... T1这种大模板啊...其实我是觉得我能打出来的,然后先用一个小时码了一个2k. 然后做T2想贪心就出来了.十分钟码完T3暴力之后回T1打对拍瞬间爆炸. 于是又重新打了一个2k,WA0.对拍发现. 然后考试就没几分钟了交暴力走了. 不要打完就跑,记得早点对拍改进思路. T1: 的确是挺裸的线段树.离散化或者权值线段树都可以. 但是考场上两个都打出来都死了. 最后用离散化A的. 1 #include<cstdio> 2

[考试反思]1002csp-s模拟测试56:凌乱

放假回来状态回升??(玩够了-但是稍困) T1打的不完全对,但是过掉了.很快的想到了二分吧喇叭啦.. 然后T2也挺快想出来了但是挂细节没发现,考试快结束的时候才发现出锅了. 改了过来是正解,但是出题人无良卡了线段树强制树状数组,T了一个子任务,卡常到飞起. T3暴力没什么问题. 卡常是一种习惯.要注意题目数据范围观察是否卡常. T1: 所有的决策都是一条一次函数. 分两类,斜率正或斜率非负. 如果第二类的直线里有在T=0时符合要求的,那么答案就是0,所以check(0)一下. 如果非负的直线都在

[考试反思]1003csp-s模拟测试58:沉淀

稳住阵脚. 还可以. 至少想拿到的分都拿到了,最后一题的确因为不会按秩合并和线段树分治而想不出来. 对拍了,暴力都拍了.挺稳的. 但是其实也有波折,险些被卡内存. 如果内存使用不连续或申请的内存全部使用的话,切记计算内存,一点都不能开大. T1: 直接根号筛,拿map也能过. 遍历map直接begin和end啊... 1 #include<cstdio> 2 int Cnt[202]; 3 struct hash_map{ 4 int cnt,fir[10000020],l[6666666],

[考试反思]0114省选模拟7:迷离

这次考得相对不错,但是没什么水准. 只不过记得T1这道原题而已.虽说我忘了怎么做,而且数据范围不一样...差不多是从头想的. 但是并没有AC,像个弱智一样,有两个细节写的完全不对还有80分运气也是真好. 其实挂了不止两个细节...以为是原题于是上来就写20分钟写完,然后过一会出一个锅... 然后看T2,感觉$O(nk^2)$也许差不多?常数很大...但也不会别的.挺好想但是不是很好写. 于是乎强烈谴责cbx没素质暴力水题考后还不改正解的无脸行径 于是就开始写,写了一个半小时. 看T3,绝对大神题

模拟测试68,69

68: 32 AlpaCa 41 03:08:20 31 02:46:16 30 02:46:28 102 03:08:20 69: 28 AlpaCa 20 02:51:15 60 03:05:32 0 01:39:45 80 03:05:32 彻底挂掉了呢. 不过也还好吧,至少之后的考试不会有那么大压力了吧(出第一机房是肯定的事了),那利用之后几场考试就调整好状态,下次再来嘛. 但是最近炸了那么多场还是要反思一下,其实不只是心态的问题,自身实力,考试技巧,时间分配等方面好像都有点问题. 69

微信在线信息模拟测试工具(基于Senparc.Weixin.MP)

目前为止似乎还没有看到过Web版的普通消息测试工具(除了官方针对高级接口的),现有的一些桌面版的几个测试工具也都是使用XML直接请求,非常不友好,我们来尝试做一个“面向对象”操作的测试工具. 测试工具在线DEMO:http://weixin.senparc.com/SimulateTool Senparc.Weixin.MP是一个开源的微信SDK项目,地址:https://github.com/JeffreySu/WeiXinMPSDK (其中https://github.com/Jeffrey

2016年上半年软考真题在线模拟测试,提前知晓你的成绩

2016年上半年软考于5月21日结束了,考试完想知道结果的急切心理,几乎每个经历过上学时代的人都能深刻体验到.如果你想知道你考的怎么样,如果你想要提前知道你的成绩,如果你想知道你哪个地方出错了,如果你想知道你哪个地方知识掌握的不够想要更深入的去理解,那就来希赛软考学院吧!希赛软考学院提供2016年上半年软考真题在线模拟测试,有标准的参考答案,有专业老师的解析视频,让你提前知晓你的成绩,让你再次巩固学习. 希赛授课专家介绍 张友生,计算机应用技术博士,软考培训教程系列丛书主编,考试指定教材<系统分

noip模拟测试11

T1:string 第一眼秒出思路,这不就是排序那道题的加强版吗? 然而歪?解复杂度虽然是对的,但常数过大,竟被卡到70 歪?解:(实际上std写的就是这个,但据说std被卡掉了 OAO) 因为字符集很小,所以我们可以把区间排序改为区间查询和覆盖 即:先查询区间内所有字符的个数,再从左端点开始按照大小关系依次将长度为字符个数的区间修改为该字符. 期望复杂度O ( 26*mlogn ),实际复杂度O ( 26*mlogn*(巨大的常数) ) 所以需要一(feng)定(kuang)的卡常 正?解: