• 等级
  • 517333 访问
  • 975 原创
  • 3 转发
  • 1414 排名
  • 66 评论
  • 63 获赞

[中等题] Project Euler 608 Divisor Sums

这个题怎么Difficultyrating80%80\%啊,送经验的感觉啊D(m,n)=======∑d|m∑k=1nσ0(kd)∑d|m∑k=1n∑a|k∑b|d[(a,b)=1]∑a=1n⌊na⌋∑d|m∑b|d[(a,b)=1]∑a=1n⌊na⌋∑d|m∑b|d∑i|a,i|bμ(i)∑i|mμ(i)×(σ0∗1)(mi)∑id≤n⌊nid⌋∑i|mμ(i)×(σ0∗1)(mi

2018-01-20 18:16:47

[简单题] Project Euler 601 Divisibility streaks

(k+1)|(n+k)(k+1)|(n+k)就是(k+1)|(n−1)(k+1)|(n-1)这个函数就是最大的kk使得1,2,⋯,k1,2,\cdots,k都整除n-1吧随便容斥下咯#include#include#includeusingnamespacestd;typedeflonglongll;inlinellGcd(lla,ll

2018-01-20 17:57:53

[简单题] Project Euler 603 Substring sums of prime concatenations

直接考虑每一位的贡献,应该是一个ai×i×(1+10+⋯+10n−i)a_i\timesi\times(1+10+\cdots+10^{n-i})的形式,这就是个等比数列求和然后因为是循环串,还是个等比数列求和,就好了#include#include#includeusingnamespacestd;typedeflonglongll;constintP

2018-01-20 17:54:56

[杜教筛] Codechef January Challenge 2018 #SQRGOOD Simplify the Square Root

二分转化为μ2\mu^2的前缀和。然后转化为O(n13)O(n^{1\over3})的运算,但是需要预处理μ\mu的前缀和,大力杜教筛求和。然后感谢阿爷教我把二分改成了迭代,小范围内一个一个挪,用rho求μ(n)\mu(n),然后就能卡进去了。复杂度似乎是萎的吧。#include#include#include#include#include#include#defi

2018-01-19 14:16:38

[整体二分] Codechef January Challenge 2018 #MONSTER Killing Monsters

整体二分,然后问题变成,子集加,单点查询,然后像CTSC吉夫特可以用经典的二进制分高位低位的搞搞。调个参,大概是高5位低12位。不知道在线怎么做。#include#include#include#include#definepbpush_backusingnamespacestd;typedeflonglongll;inlinecharnc(){

2018-01-19 14:11:46

[后缀数组 后缀树] Codechef January Challenge 2018 #KILLKTH Killjee and k-th letter

建出后缀树,记录每个子串的出现次数,然后二分下答案在哪个子串中就好了退役选手不会写后缀自动机#include#include#include#include#definepbpush_backusingnamespacestd;typedeflonglongll;inlinecharnc(){staticcharbuf[100000],*p1=b

2018-01-19 14:07:23

高三日志

原标题高三志开坑2017.09.0223:47高三,还是记录一下吧。DATEEVENTHINT07.22退役07.25-07.27CodeMFinal,Beijing点赞08.04-08.18CodechefAugustLongChallenge2017高三可能都不会这么认真的打比...

2017-09-03 00:08:24

AFO & 博客已搬 manchery.co

抱歉,退役了,Blog下的留言我可能不会回答了吧。UPD.博客已搬[雯舞-Manchery'sSpace](manchery.co)——2018.05.25

2017-07-24 10:19:51

[数学 FFT] Codechef July Challenge 2017 #APRPS Irrational Root

跟51Nod1356代数数的次数是一样的不过这里都是质数也就是就是2n2^n关键是输方案这个不一定有二次剩余感谢sxt一个一个数加进答案转化成已知F(x)F(x),求F(x+a√)F(x+\sqrta)和F(x−a√)F(x-\sqrta)的系数这个推一下就是一个FFT拷了myy的板子自己的太慢了QAQ#include<cstdio>#include<cs

2017-07-24 10:18:34

[杂题] Codechef July Challenge 2017 #MULDIG Multiplication Program

大意是把二元运算化成一元运算感谢sxt的帮助剩下的看代码注释吧//我觉得蛮清楚了唯一不爽的是没想到临时变量这么多mmp#include<cstdio>#include<cstdlib>#include<algorithm>#include<vector>#definepbpush_backusingnamespacestd;intB;structabcd{i

2017-07-24 10:12:29

[杂题] Codeforces 830C Round #424 Div1 C. Bamboo Partition

随便搞搞吧#include<cstdio>#include<cstdlib>#include<algorithm>usingnamespacestd;typedeflonglongll;#defineread(x)scanf("%d",&(x))constintmaxn=200005;constintP=1e9+7;constintN=505;intn,a[N]

2017-07-24 10:05:45

[FWT] UOJ #310. 【UNR #2】黎明前的巧克力

这是若干个2xai+12x^{a_i}+1的东西的卷积然后这个FWT一下发现每一项只有−1-1或33那么卷积的FWT每一项就是若干个−1-1和33的乘积这个不好求直接加在一起FWT,那么我们得到了每一项−1-1和33的和因为只有这两个取值,可以直接解方程有多少个−1-1和33然后快速幂再乘起来#include<cstdio>#include<

2017-07-24 09:59:50

[DP] UOJ #311. 【UNR #2】积劳成疾

fi,jf_{i,j}表示长为ii的区间最大值是jj的答案转移就枚举最左边的最大值在区间的位置前缀和优化下就好了好像也可以fi,j,kf_{i,j,k}表示前ii个,末尾KK个中最大值在jj,最大值是kk,有点复杂?#include<cstdio>#include<cstdlib>#include<algorithm>usingnamespacestd;typ

2017-07-23 14:10:32

[色多项式] UOJ #308. 【UNR #2】UOJ拯救计划 & SRM 717 div1 AcyclicOrientation

一个图的kk染色数是关于kk的nn次(?)多项式称为色多项式那么这里模6我们只要知道模2和模3的值然后分类讨论下就好了一张图的0染色数是0,1染色数等于[m=0][m=0],2染色数与二分图的联通块个数有关#include<cstdio>#include<cstdlib>#include<algorithm>#include<cstring>#define

2017-07-23 14:07:33

[线代小记] 树形图求和

未经允许搬了过来代码略丑#include<cstdio>#include<cstdlib>#include<algorithm>usingnamespacestd;typedeflonglongll;inlinecharnc(){staticcharbuf[100000],*p1=buf,*p2=buf;returnp1==p2&&(p2=(p1=buf)+f

2017-07-23 13:55:01

[高斯消元 矩阵的秩] 51Nod 1356 代数数的次数

请教了比利首先2n2^n次项是没有问题的,比如2√+3√\sqrt2+\sqrt3,可以构造f(x)=(x+2√+3√)(x+2√−3√)(x−2√+3√)(x−2√−3√)f(x)=(x+\sqrt2+\sqrt3)(x+\sqrt2-\sqrt3)(x-\sqrt2+\sqrt3)(x-\sqrt2-\sqrt3)这个不是最优解最优解应该是构造这样一个矩阵每一行代表一个质因子

2017-07-23 13:50:12

[数学 二分图匹配] SRM 456 div1 FunctionalEquation

本来想自己再推一下的,但是退役了也就弃坑了未经允许的搬了搬题人的题解//BEGINCUTHERE#include<conio.h>#include<sstream>//ENDCUTHERE#include<cstdio>#include<cstdlib>#include<cstring>#include<cmath>#include<iostream>#in

2017-07-23 13:41:16

[杂题] Codechef SnackDown 2017 Onsite Final #MINIMAX Minimax

最小值的最大值小于等于最大值的最小值那么只要有一行一列最小值等于最大值,那么矩阵就合法如果我们枚举这一行和这一列,那么最小值最大值一定就是他们的交点那么再枚举交点最后的值是多少这样复杂度是O(n4)O(n^4)换一下先枚举值再枚举位置可以去一个代价最小值O(n2logn)O(n^2\logn)#include<cstdio>#include<cstdlib>#in

2017-07-23 13:32:53

[数论] LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具

这种非简单路一般转化成任意一条路加上若干环这里大概是任意一条路加上若干环长的gcd任意一条路可以弄出任意一颗生成树?但是这里实际上只需要一个带权并查集详见官方题解#include<cstdio>#include<cstdlib>#include<algorithm>usingnamespacestd;typedeflonglongll;typedefpair<

2017-07-23 13:26:30

[反演] LOJ #509. 「LibreOJ NOI Round #1」动态几何问题

μ2(n)=∑d2|nμ(d)\mu^2(n)=\sum_{d^2|n}\mu(d)然后就是xjb推反正退役了我也就弃坑了95分代码复杂度分析及优化详见官方题解#include<cstdio>#include<cstdlib>#include<cmath>#include<algorithm>usingnamespacestd;typedeflonglongl

2017-07-23 13:22:57

里阿奴摩西

关注
奖章
  • 持之以恒