<?xml version="1.0" encoding="utf-8" ?><rss version="2.0"><channel><title><![CDATA[雯舞]]></title><description><![CDATA[Love three things int he world -- Manchery]]></description><link>https://blog.youkuaiyun.com/u014609452</link><language>zh-cn</language><generator>https://blog.youkuaiyun.com/</generator><copyright><![CDATA[Copyright &copy; u014609452]]></copyright><item><title><![CDATA[[中等题] Project Euler 608 Divisor Sums]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79115841</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79115841</guid><author>u014609452</author><pubDate>Sat, 20 Jan 2018 18:16:47 +0800</pubDate><description><![CDATA[这个题怎么Difficulty rating 80%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]]></description><category></category></item><item><title><![CDATA[[简单题] Project Euler 601 Divisibility streaks]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79115753</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79115753</guid><author>u014609452</author><pubDate>Sat, 20 Jan 2018 17:57:53 +0800</pubDate><description><![CDATA[(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
#include
using namespace std;
typedef long long ll;

inline ll Gcd(ll a,ll]]></description><category></category></item><item><title><![CDATA[[简单题] Project Euler 603 Substring sums of prime concatenations]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79115741</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79115741</guid><author>u014609452</author><pubDate>Sat, 20 Jan 2018 17:54:56 +0800</pubDate><description><![CDATA[直接考虑每一位的贡献，应该是一个 ai×i×(1+10+⋯+10n−i)a_i\times i\times (1+10+\cdots+10^{n-i}) 的形式，这就是个等比数列求和 
 然后因为是循环串，还是个等比数列求和，就好了

#include
#include
#include
using namespace std;
typedef long long ll;

const int P]]></description><category></category></item><item><title><![CDATA[[杜教筛] Codechef January Challenge 2018 #SQRGOOD Simplify the Square Root]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79106226</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79106226</guid><author>u014609452</author><pubDate>Fri, 19 Jan 2018 14:16:38 +0800</pubDate><description><![CDATA[二分转化为μ2\mu^2的前缀和。 
然后转化为O(n13)O(n^{1\over 3})的运算，但是需要预处理μ\mu的前缀和，大力杜教筛求和。 
然后感谢阿爷教我把二分改成了迭代，小范围内一个一个挪，用rho求μ(n)\mu(n)，然后就能卡进去了。 
复杂度似乎是萎的吧。

#include
#include
#include
#include
#include
#include
#defi]]></description><category></category></item><item><title><![CDATA[[整体二分] Codechef January Challenge 2018 #MONSTER Killing Monsters]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79106152</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79106152</guid><author>u014609452</author><pubDate>Fri, 19 Jan 2018 14:11:46 +0800</pubDate><description><![CDATA[整体二分，然后问题变成，子集加，单点查询，然后像CTSC吉夫特 可以用经典的二进制分高位低位的搞搞。调个参，大概是高5位低12位。 
不知道在线怎么做。



#include
#include
#include
#include
#define pb push_back
using namespace std;
typedef long long ll;

inline char nc(){]]></description><category></category></item><item><title><![CDATA[[后缀数组 后缀树] Codechef January Challenge 2018 #KILLKTH Killjee and k-th letter]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/79106100</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/79106100</guid><author>u014609452</author><pubDate>Fri, 19 Jan 2018 14:07:23 +0800</pubDate><description><![CDATA[建出后缀树，记录每个子串的出现次数，然后二分下答案在哪个子串中就好了

退役选手不会写后缀自动机

#include
#include
#include
#include
#define pb push_back
using namespace std;
typedef long long ll;

inline char nc(){
  static char buf[100000],*p1=b]]></description><category></category></item><item><title><![CDATA[高三日志]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/77807566</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/77807566</guid><author>u014609452</author><pubDate>Sun, 03 Sep 2017 00:08:24 +0800</pubDate><description><![CDATA[原标题 高三志

开坑 2017.09.02 23:47 
高三，还是记录一下吧。




  DATE
  EVENT
  HINT



  07.22
  退役
  


  07.25-07.27
  CodeM Final，Beijing
  点赞


  08.04-08.18
  Codechef August Long Challenge 2017
  高三可能都不会这么认真的打比...]]></description><category></category></item><item><title><![CDATA[AFO & 博客已搬 manchery.co]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75981677</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75981677</guid><author>u014609452</author><pubDate>Mon, 24 Jul 2017 10:19:51 +0800</pubDate><description><![CDATA[抱歉，退役了，Blog下的留言我可能不会回答了吧。

UPD. 博客已搬 [雯舞 - Manchery's Space](manchery.co) ——2018.05.25]]></description><category></category></item><item><title><![CDATA[[数学 FFT] Codechef July Challenge 2017 #APRPS Irrational Root]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75980207</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75980207</guid><author>u014609452</author><pubDate>Mon, 24 Jul 2017 10:18:34 +0800</pubDate><description><![CDATA[跟 51Nod 1356 代数数的次数 是一样的 
不过这里都是质数 也就是就是 2n2^n  
关键是输方案 这个不一定有二次剩余感谢sxt 
一个一个数加进答案 转化成 已知F(x)F(x)，求F(x+a√)F(x+\sqrt a)和F(x−a√)F(x-\sqrt a)的系数 
这个推一下就是一个FFT拷了myy的板子 自己的太慢了QAQ#include<cstdio>
#include<cs]]></description><category></category></item><item><title><![CDATA[[杂题] Codechef July Challenge 2017 #MULDIG Multiplication Program]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75976809</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75976809</guid><author>u014609452</author><pubDate>Mon, 24 Jul 2017 10:12:29 +0800</pubDate><description><![CDATA[大意是把二元运算化成一元运算 
感谢sxt的帮助剩下的看代码注释吧 //我觉得蛮清楚了 
唯一不爽的是 没想到临时变量这么多 mmp#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#define pb push_back
using namespace std;int B;struct abcd{
  i]]></description><category></category></item><item><title><![CDATA[[杂题] Codeforces 830C Round #424 Div1 C. Bamboo Partition]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75974550</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75974550</guid><author>u014609452</author><pubDate>Mon, 24 Jul 2017 10:05:45 +0800</pubDate><description><![CDATA[随便搞搞吧#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;#define read(x) scanf("%d",&(x))const int maxn=200005;
const int P=1e9+7;const int N=505;int n,a[N]]]></description><category></category></item><item><title><![CDATA[[FWT] UOJ #310. 【UNR #2】黎明前的巧克力]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75972013</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75972013</guid><author>u014609452</author><pubDate>Mon, 24 Jul 2017 09:59:50 +0800</pubDate><description><![CDATA[这是若干个 2xai+12x^{a_i}+1 的东西的卷积 
然后这个FWT一下发现每一项只有 −1-1 或 33 
那么卷积的FWT每一项就是若干个 −1-1 和 33 的乘积 
这个不好求 
直接加在一起FWT，那么我们得到了每一项 −1-1 和 33 的和 
因为只有这两个取值，可以直接解方程有多少个 −1-1 和 33 
然后快速幂再乘起来#include<cstdio>
#include<]]></description><category></category></item><item><title><![CDATA[[DP] UOJ #311. 【UNR #2】积劳成疾]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75908822</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75908822</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 14:10:32 +0800</pubDate><description><![CDATA[fi,jf_{i,j}表示长为 ii 的区间 最大值是jj 的答案 
转移就枚举最左边的最大值在区间的位置 前缀和优化下就好了好像也可以fi,j,kf_{i,j,k}表示前 ii 个，末尾 KK 个中最大值在 jj ，最大值是 kk，有点复杂？#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typ]]></description><category></category></item><item><title><![CDATA[[色多项式] UOJ #308. 【UNR #2】UOJ拯救计划 & SRM 717 div1 AcyclicOrientation]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75908287</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75908287</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 14:07:33 +0800</pubDate><description><![CDATA[一个图的 kk 染色数是关于 kk 的 nn 次(?) 多项式 
称为色多项式  
那么这里模6 我们只要知道模2和模3的值  
然后分类讨论下就好了 
一张图的0染色数是0，1染色数等于[m=0][m=0]，2染色数与二分图的联通块个数有关#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define]]></description><category></category></item><item><title><![CDATA[[线代小记] 树形图求和]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907411</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907411</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:55:01 +0800</pubDate><description><![CDATA[未经允许搬了过来代码略丑#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+f]]></description><category></category></item><item><title><![CDATA[[高斯消元 矩阵的秩] 51Nod 1356 代数数的次数]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907391</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907391</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:50:12 +0800</pubDate><description><![CDATA[请教了比利 
首先 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) 
这个不是最优解 
最优解应该是构造这样一个矩阵 每一行代表一个质因子]]></description><category></category></item><item><title><![CDATA[[数学 二分图匹配] SRM 456 div1 FunctionalEquation]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907295</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907295</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:41:16 +0800</pubDate><description><![CDATA[本来想自己再推一下的，但是退役了也就弃坑了 
未经允许的搬了搬题人的题解// BEGIN CUT HERE  
#include<conio.h>
#include<sstream>
// END CUT HERE  
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#in]]></description><category></category></item><item><title><![CDATA[[杂题] Codechef SnackDown 2017 Onsite Final #MINIMAX Minimax]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907167</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907167</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:32:53 +0800</pubDate><description><![CDATA[最小值的最大值小于等于最大值的最小值 
那么只要有一行一列最小值等于最大值，那么矩阵就合法 
如果我们枚举这一行和这一列，那么最小值最大值一定就是他们的交点 
那么再枚举交点最后的值是多少 
这样复杂度是O(n4)O(n^4) 
换一下 先枚举值 再枚举位置 可以去一个代价最小值 
O(n2logn)O(n^2\log n)#include<cstdio>
#include<cstdlib>
#in]]></description><category></category></item><item><title><![CDATA[[数论] LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907105</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907105</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:26:30 +0800</pubDate><description><![CDATA[这种非简单路 一般转化成任意一条路加上若干环 
这里大概是任意一条路加上若干环长的gcd 
任意一条路 可以弄出任意一颗生成树？ 但是这里实际上只需要一个带权并查集 
详见官方题解#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef pair<]]></description><category></category></item><item><title><![CDATA[[反演] LOJ #509. 「LibreOJ NOI Round #1」动态几何问题]]></title><link>https://blog.youkuaiyun.com/u014609452/article/details/75907075</link><guid>https://blog.youkuaiyun.com/u014609452/article/details/75907075</guid><author>u014609452</author><pubDate>Sun, 23 Jul 2017 13:22:57 +0800</pubDate><description><![CDATA[μ2(n)=∑d2|nμ(d)\mu^2(n)=\sum_{d^2|n} \mu(d) 
然后就是xjb推  
反正退役了 我也就弃坑了 95分代码 
复杂度分析及优化详见官方题解#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long l]]></description><category></category></item></channel></rss>