码迷,mamicode.com
首页 > 其他好文 > 详细

「考试」省选73

时间:2020-04-17 15:20:26      阅读:71      评论:0      收藏:0      [点我收藏+]

标签:nlog   isp   方案   nlogn   art   最大   情况   网络流   rac   

T1
看起来是个状压\(dp\)
\(dp[i][j][S]\)为前\(i\)行,放置了\(j\)个哲学家,第\(i\)行放的情况是\(S\)的方案数。
那么对于两个状态,能否由上一行转移到这一行会存在一些判断,这些判断暴力写一写就行了,也不是很恶心。
对于两个状态能够转移的话\(s_1\rightarrow s_2\),那我们就用状压\(dp\)暴力的在两行之间转移就可以了。
复杂度最劣是\(O(3n^22^6)\)的。
卡常可以过。
比如说对于一个\(i\)\(j\)枚举的下界是\(max(0,m-3(n-i+1))\)这样子。
然后看了一下神仙正解。
如果我们把\(dp(i,S)\)看成是\(\{dp[i][S][j]\}\)的生成函数。
那么有:\([x^j]dp(i,S)=dp[i][j][S]\)
同时可以用矩阵进行转移。
\(dp(i)\)为一个列向量,每一个位置存\(dp(i,S)\)就可以了。
然后转移矩阵\(A\)有:\(s_1\rightarrow s_2,a_{s_1,s_2}=x^{bit[s_2]}\)
这样就可以直接用矩阵快速幂了。
然而你直接用矩阵暴力乘多项式上去的复杂度就是\(O(512nlog^2n)\)的。
还不如暴力。
考虑直接用点值乘,这样省去了一个\(logn\)
\(A^{n}\)的复杂度就是:\(O(512nlogn)\)的了。

T2
是01分数规划。
我们考虑如何\(check\)一个答案。
相当于是有\(Amid-B>0\)则偏大,否则偏小。
可以二分。
\(check\)的时候算出来\(Amid-B\)的最大值就可以了,判断最大值是不是偏大。
很显然的一个最大权闭合子图。
两个点之间的贡献拆成两倍来做,这样\(a[i]\)也要乘二才能保证上下相除答案不变。
然后对于两个点的贡献分开割就可以了。
网络流得卡常。
我写了当前弧优化,然后被迫又把两条边合并。
这才卡过去。

T3
很奇怪的找到三种算法(两种可行)。
1.

\[ans=\sum\limits_{i=1}^{n}2^{f(i)} \]

\[g(n)=2^{f(n)} \]

\[(a,b)=1,f(ab)=f(a)+f(b),g(ab)=g(a)g(b),g(p^e)=2p^0 \]

这个东西可以发现是直接\(Min\_25\)筛就行。但是太难写了。
所以我后面两个都是想推出一个复杂度比较优秀的暴力出来。

\[g(d)=2^{f(d)}=\sum\limits_{g|d}[(g,\frac{d}{g})=1] \]

\[\begin{aligned} ans&=\sum\limits_{i=1}^{n}g(i)\&=\sum\limits_{i=1}^{n}\sum\limits_{d|i}[(d,\frac{i}{d})=1]\&=\sum\limits_{i=1}^{n}\sum\limits_{d|i}\sum\limits_{g|(d,\frac{i}{d})}\mu(g)\&=\sum\limits_{d=1}^{n}\sum\limits_{i=1}^{\frac{n}{d}}\sum\limits_{g|(d,i)}\mu(g)\&=\sum\limits_{g=1}^{n}\mu(g)\sum\limits_{d=1}^{\frac{n}{g}}\frac{n}{dg}\h(n)&=\sum\limits_{i=1}^{n}\frac{n}{i}\&=\sum\limits_{g=1}^{n}\mu(g)h(\frac{n}{g})\\end{aligned} \]

前面的部分整除分块+杜教筛就可以算出来,后面的\(h\)也是可以用整除分块的,复杂度在线性左右(不会证)。
这个算法中间的地方有一些性质没有挖掘出来。
其实本质上是和最终算法是相同的。
只不过最终算法把这个算法的枚举给合并同类项了。

\[g(d)=2^{f(d)}=\sum\limits_{g|d}[(g,\frac{d}{g})=1] \]

\[\begin{aligned} ans&=\sum\limits_{i=1}^{n}g(i)\&=\sum\limits_{i=1}^{n}\sum\limits_{d|i}[(d,\frac{i}{d})=1]\&=\sum\limits_{i=1}^{n}\sum\limits_{d|i}\sum\limits_{g|(d,\frac{i}{d})}\mu(g)\&=\sum\limits_{i=1}^{n}\sum\limits_{g^2|i}\mu(g)d(\frac{i}{g^2})\&=\sum\limits_{g=1}^{\sqrt{n}}\mu(g)\sum\limits_{i=1}^{\frac{n}{g^2}}d(i)\&=\sum\limits_{g=1}^{\sqrt{n}}\mu(g)\sum\limits_{i=1}^{\frac{n}{g^2}}\sum\limits_{d|i}I\&=\sum\limits_{g=1}^{\sqrt{n}}\mu(g)\sum\limits_{d=1}^{\frac{n}{g^2}}\frac{n}{g^2d}\&=\sum\limits_{g=1}^{\sqrt{n}}\mu(g)l(\frac{n}{g^2})\l(n)&=\sum\limits_{i=1}^{n}\frac{n}{i}\\end{aligned} \]

前面枚举\(g\),然后分块算\(l(\frac{n}{g^2})\)即可。
这个算法的话时间复杂度也很优秀。

\[O=\int_1^{\sqrt{n}}\sqrt{\frac{n}{x^2}}dx=\sqrt{n}\sum\limits_{x=1}^{\sqrt{n}}\frac{1}{x}dx=\sqrt{n}ln(\sqrt{n}) \]

这样就解决了。

「考试」省选73

标签:nlog   isp   方案   nlogn   art   最大   情况   网络流   rac   

原文地址:https://www.cnblogs.com/Lrefrain/p/12720071.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!