- 博客(565)
- 收藏
- 关注
RSS订阅转载 Markdown 设置字体大小颜色及背景色
一、更改字体、大小、颜色<font face="黑体">我是黑体字</font><font face="微软雅黑">我是微软雅黑</font><font face="STCAIYUN">我是华文彩云</font><font color=red>我是红色</font><font c...
2019-09-27 10:29:00
1132
转载 怎样把txt文档转换成csv文件？
其实csv就是逗号隔开的一行一行的数据，如果每行数据中都是用逗号分隔的，直接把文件后缀txt改成csv就行了。用python搞定：import numpy as npimport pandas as pd txt = np.loadtxt('file.txt')txtDF = pd.DataFrame(txt)txtDF.to_csv('file.cs...
2019-09-24 22:02:00
11908
转载 【模板】A*B Problem（FFT快速傅里叶）
题目：给出两个n位10进制整数x和y，你需要计算x*y。（$n \leq 60000$）分析：两个正整数的相乘可以视为两个多项式的相乘，例如 $15 \times 16 = 240$，可写成 $(5+x)*(6+x) = 30 + 11x + x^2$，$x=10$这样得到多项式 $A(x)$ 和 $B(x)$，并且能用FFT求出 $C(x)=A(x)B(x)$，...
2019-09-23 22:50:00
240
转载 FFT入门
可以先了解下用分治法求多项式的值 链接。求 $A(x)$ 在 $w_j,j=0,1,2...,n-1$ 处的值：设多项式$$A(x)=a_0+a_1*x+a_2*{x^2}+a_3*{x^3}+a_4*{x^4}+a_5*{x^5}+ \dots+a_{n-2}*x^{n-2}+a_{n-1}*x^{n-1}$$按下标的奇偶性分类，设$$A_1(x)=a_0+a_2*...
2019-09-23 22:36:00
161
转载 分治法求多项式的值
问题：设多项式 $A(x) = a_0 + a_1x + ...+a_{n-1}x^{n-1}$，求多项式在某点的值。分析：将多项式按奇偶分类：设 $A_{even}(x)$ 为偶数项系数构造的多项式，$A_{odd}(x)$ 为奇数项系数的多项式$\begin{aligned}A(x) &= a_0 + a_1x + ...+a_{n-1}x^{n-1} \\&...
2019-09-23 15:54:00
563
转载 [TJOI2015]概率论——期望&&母函数
题意求一个含有 $n$ 个结点的有序二叉树的叶子节点的期望个数。（$n \leq 10^9$）分析一堆推导.....得 $ans = \frac{n(n+1)}{2(2n-1)}$#include<bits/stdc++.h>using namespace std;double n;int main(){ scanf(...
2019-09-22 15:22:00
253
转载 hdu1171&&P2000——母函数
hdu1171题意：有 $n$ 种设施，每种有价值 $v_i$ 和数量 $m_i$，求一种方案使得分成价值尽可能相近的两组。（$n \leq 50, v_i \leq 50, m_i \leq 100$）分析：可以用背包做，这里讲母函数的做法。直接用样例说明一下：310 1 20 230 1其母函数为$$(1+x^{10})(1+x^{20}+x^{40})...
2019-09-22 14:29:00
125
转载 母函数入门【模板】
正整数拆分hdu1028解：对于正整数 $n$ 的拆分，其母函数为$$f(x) = (1+x+x^2+...)(1+x^2+x^4+...)(1+x^3+x^6+x^9+...)...$$答案就是多项式展开后 $x^n$ 项的系数。Code://其实就是模拟，从前往后一一合并#include<bits/stdc++.h>using na...
2019-09-21 23:31:00
76
转载 BZOJ2159 Crash的文明世界——树上DP&&第二类Stirling数
题意给定一个有 $n$ 个结点的树，设 $S(i)$ 为第 $i$ 个结点的“指标值”，定义为 $S(i)=\sum_{i=1}^{n}dist(i,j)^k$，$dist(i, j)$ 为结点 $i$ 到结点 $j$ 的最小距离。请输出每个结点的指标值。($n \leq 5000, k \leq 150$)分析一个常用的转化$$n^k=\sum_{i=0}^{k}S(k...
2019-09-21 19:49:00
107
转载 hdu3625
hdu3625题意：酒店发生一起谋杀案。作为镇上最好的侦探，您应该立即检查酒店的所有N个房间。但是，房间的所有门都是锁着的，钥匙刚锁在房间里，真是个陷阱！您知道每个房间里只有一把钥匙，并且所有可能的分配可能性均等。例如，如果N = 3，则有6种可能的分布，每种分布的概率为1/6。为方便起见，我们将房间编号从1到N，房间1的键编号为键1，房间2的键编号为2，依此类推。要检查所有房间...
2019-09-21 15:31:00
134
转载 Stirling数入门
第一类Stirling数定义$$\begin{aligned}(x)_n & =x(x-1)...(x-n+1)\\&= s(n, 0) + s(n,1)x +..+s(n,n)x^n\\\end{aligned}$$例如，$n=3$ 时，$(x)3 = x(x-1)(x-2)$$(x)3 = x^0 + 2x -3x^2 + x^3$于是 $s(...
2019-09-21 14:32:00
143
转载 自然数幂和——第一类Stirling数和第二类Stirling数
第一类Stirling数首先设$$S_k(n)=\sum_{i=0}^ni^k$$根据第一类斯特林数的定义（P是排列数，C是组合数,s是Stirling）$$C_n^k={P_n^k\over k!}={\sum_{i=0}^k(-1)^{i+k}s(k,i)n^i\over k!}$$变形得$$n^k ={\sum_{i=0}^{k-1}(-1)^{i+k}s...
2019-09-21 09:43:00
390
转载 51Node1228序列求和 ——自然数幂和模板&&伯努利数
伯努利数法伯努利数原本就是处理等幂和的问题，可以推出$$ \sum_{i=1}^{n}i^k={1\over{k+1}}\sum_{i=1}^{k+1}C_{k+1}^i*B_{k+1-i}*(n+1)^i $$因为$$\sum_{k=0}^nC_{n+1}^kB_k=0(B_0=1)$$所以$$B_n={- {1\over{n+1}}}(C_{n+1}^0B_...
2019-09-20 21:50:00
200
转载 CF622F——自然数幂和模板&&拉格朗日插值
题意求 $ \displaystyle \sum_{i=1}^ni^k \ mod (1e9+7), n \leq 10^9, k \leq 10^6$.CF622F分析易知答案是一个 $k+1$ 次多项式，我们找 $k+2$ 个值代进去，然后拉格朗日插值。$n+1$ 组点值对 $(x_i, y_i)$，得到 $n$ 次多项式 $f$ 的拉格朗日插值公式为：$$...
2019-09-20 11:58:00
284
转载 Bell数入门
贝尔数贝尔数是以埃里克·坦普尔·贝尔命名，是组合数学中的一组整数数列，开首是（OEIS的A000110数列）：$$B_0 = 1, B_1 = 1, B_2 = 2, B_3 = 5, B_4 = 15, B_5 = 52, B_6 = 203, ...$$$B_n$ 的含义是基数为 $n$ 的集合划分成非空集合的划分数。例如， $B_3=5$ 是因为3个元素的集合有5种...
2019-09-19 14:49:00
695
转载 hdu2643&&hdu2512——斯特林数&&贝尔数
hdu2643题意：$n$ 个人的排名情况数（$n \leq 100$）分析：考虑 $n$ 个有区别的球放到 $m$ 个有区别的盒子里、无空盒的方案数为 $m!\cdot S(n, m)$。这题中 $m$ 可取 $1 \sim n$（可能排名相同），累加即可。#include<bits/stdc++.h>using namespace std;...
2019-09-19 11:23:00
98
转载 hdu4767 Bell——求第n项贝尔数
题意设第 $n$ 个Bell数为 $B_n$，求 $B_n \ mod \ 95041567$.（$1 \leq n \leq 2^{31}$）分析贝尔数的概念和性质，维基百科上有，这里用到两点。若 $p$ 是任意素数，有 $B_{p+n} = B_n + B_{n+1}(mod \ p)$每个贝尔数都是相应第二类斯特林数之和，即 $B_n = \sum_...
2019-09-18 23:25:00
471
转载 Uva11762 Race to 1——有向无环图&&记忆化搜索
题意给出一个整数 $N$，每次可以在不超过 $N$ 的素数中等概率随机选择一个 $P$，如果 $P$ 是 $N$ 的约数，则把 $N$ 变成 $N/P$，否则 $N$ 不变。问平均情况下需要多少次随机选择，才能把 $N$ 变成1呢？分析本题可以画出一个状态转移图，例如 $n=6$ 时，$n$ 的每个约数都对应一个状态，每个状态转移都有一定概率，从每个状态出发转移的...
2019-09-18 17:18:00
119
转载 P3232 [HNOI2013]游走——无向连通图&&高斯消元
题意一个无向连通图，顶点从1编号到N，边从1编号到M。 小Z在该图上进行随机游走，初始时小Z在1号顶点，每一步小Z以相等的概率随机选 择当前顶点的某条边，沿着这条边走到下一个顶点，获得等于这条边的编号的分数。当小Z 到达N号顶点时游走结束，总分为所有获得的分数之和。 现在，请你对这M条边进行编号，使得小Z获得的总分的期望值最小。（2<=N<=500）分析直接算边的...
2019-09-18 16:03:00
154
转载 Random Walk——高斯消元法
题目有一个 $N \times M$ 大小的格子，从(0, 0)出发，每一步朝着上下左右4个格子中可以移动的格子等概率移动。另外有些格子有石头，因此无法移至这些格子。求第一次到达 $(N-1, M-1)$ 格子的期望步数。（$2 \leq N,M\leq 10$）分析设 $E(x, y)$ 表示从 (x, y) 出发到终点的期望步数。我们先考虑从 $(x, y)$ 向上...
2019-09-18 08:37:00
215
转载 B君的历史——复数乘法&&爆搜
题意设 $r = \frac{-1+\sqrt7 i}{2}$，对任意整数 $x, y$ 都可以找到一个有限的整数集合 $S$，使得$$x + y\sqrt7 i = \sum_{k \in S }r^k$$换句话说，就是将 $x + y \sqrt7 i$ 转换成 $r$ 进制，集合中不能包含相同的元素，可以证明 $S$ 是唯一的。输入 $x$ 和 $y$，从小到大输出...
2019-09-17 16:18:00
64
转载 复数快速幂【模板】
实数范围内的乘法都是满足结合律的，像 $1+\sqrt2,1+2i, 1+2\sqrt2 i$ 这些没有直接可用的乘法，但是相乘后形式保持不变，所以可以模拟乘法。#include<bits/stdc++.h>using namespace std;typedef long long ll;const ll mod = 1e9 + 7;struc...
2019-09-17 16:02:00
296
转载 UVa11542Squre——异或方程组&&高斯消元法
题意给出 $n$ 个整数，从中选出1个或多个，使得选出的整数乘积是完全平方数。一共有多少种选法？（$1 \leq n \leq 100$，$1 \leq a_i \leq 10^{15}$ 且不含大于500的素因子）分析“不含大于500的素因子”提示我们考虑每个数的素因数分解，用01向量表示一个数，再用 $n$ 个01变量 $x_i$ 来表示我们的选择，其中 $x_i=1$ ...
2019-09-17 09:32:00
254
转载 UVa10828 Back to Kernighan-Ritchie——概率转移&&高斯消元法
题意给出一个程序控制流图，从每个结点出发到每个后继接结点的概率均相等。当执行完一个没有后继的结点后，整个程序终止。程序总是从编号为1的结点开始。你的任务是对于若干个查询结点，求出每个结点的期望执行次数。结点个数 $n < 10$.分析如果是有向无环图，可以直接解出递推关系，再采用记忆化搜索求解。当这题可能有环，只能列出方程，用高斯消元解方程组。设结点 $i$ 的...
2019-09-16 17:07:00
97
转载 高斯消元法【模板】
高斯消元法，消成行阶梯型矩阵。下面两种消元法的时间复杂度都是 $O(n^3)$.#include<bits/stdc++.h>using namespace std;const int maxn = 100+10;typedef double Matrix[maxn][maxn];//要求系数矩阵可逆//这里的A是增广矩阵，即A[i][...
2019-09-16 11:40:00
115
转载 转成p进制算法C语言
今天打比赛的时候竟然下一没有想起来，实际上是非常简单的。举例说明：$64 = 2 \times 3^3 + 1 \times 3^2 + 3^0$根据秦九韶算法每次提出3，即 $3(2 \times 3^2 + 1 \times 3^1) + 1 = 64$,可见得到最低位系数只需要对3取余，将n除以3，重复上述步骤。#include<stdio.h...
2019-09-16 09:27:00
185
转载 2019牛客多校第九场AThe power of Fibonacci——扩展BM
题意求斐波那契数列m次方的前n项和，模数为 $1e9$.分析线性递推乘线性递推仍是线性递推，所以上BM。由于模数非质数，上扩展版的BM。递推多少项呢？本地输入发现最大为与前57项有关（而且好像有周期性？）；当然这个算法飞快，搞1000都没问题。#include <cstdio>#include <cstdlib>#inclu...
2019-09-15 22:22:00
195
转载 BM递推杜教版【扩展】
也就是模数不是质数的时候，//下面的板子能求质数和非质数，只需要传不同的参数。#include <cstdio>#include <cstdlib>#include <cassert>#include <cstring>#include <bitset>#include <cmath&...
2019-09-15 22:13:00
111
转载 2019牛客多校第二场BEddy Walker 2——BM递推
题意从数字 $0$ 除法，每次向前走 $i$ 步，$i$ 是 $1 \sim K$ 中等概率随机的一个数，也就是说概率都是 $\frac{1}{K}$。求落在过数字 $N$ 额概率，$N=-1$ 表示无穷远。分析设落在过 $i$ 的概率为 $p_i$，则 $p_i = \frac{1}{K}p_{i-1} + \frac{1}{K}p_{i-2}...+\frac{1}{K}...
2019-09-15 20:54:00
173
转载 hdu6172&&hdu6185&&P5487——BM算法
hdu6172模板的简单应用先根据题中的表达式求出前几项，再上BM，注意一下n的大小关系。#include <bits/stdc++.h>using namespace std;#define rep(i,a,n) for (long long i=a;i<n;i++)#define per(i,a,n) for (long long...
2019-09-15 17:09:00
181
转载 铺砖头问题（完美）——爆搜&&插头DP
题意给定一个 $n \times m$ 的格子，每个格子被染成了黑色或白色。现在要用 $1 \times 2$ 的砖块覆盖这些格子，要求块与块之间互不重叠，且覆盖了所有白色的格子，但不覆盖任意黑色格子。求一共有多少种覆盖方法。结果对 $M$ 取余。（$1 \leq n\leq 15, 1 \leq m\leq 15$）分析爆搜，从最左上方开始放置，从左至右，从上到下，每找到一...
2019-09-15 10:36:00
189
转载 BM算法【实数模板】
BM递推杜教版是在整数取模的情况下的，这个可以求解实数系数，但是可能有精度误差。若一个问题的结论是通过推线性递推式来解，考虑到实际的情况，可以用BM算法的模板，先输入项数再依次输入项，项越多越准确（按道理，前k项的递推，只需要2*k 个初始项就能确定）#include <bits/stdc++.h>using namespace std;#define...
2019-09-15 00:25:00
127
转载 BM递推杜教版
#include <bits/stdc++.h>using namespace std;#define rep(i,a,n) for (long long i=a;i<n;i++)#define per(i,a,n) for (long long i=n-1;i>=a;i--)#define pb push_back#define mp...
2019-09-15 00:11:00
104
转载 DockerToolbox安装docker - Windows 10家庭版
由于本机使用的是win10家庭版操作系统，无法直接Docker for Windows安装，因此只好使用Docker Toolbox。在此记录一下过程，以供参考。下载因为toolbox安装包的官网下载需要翻墙,所以从阿里云提供的镜像地址下载.地址:http://mirrors.aliyun.com/docker-toolbox/windows/或者去github：http...
2019-09-14 20:08:00
453
转载 npm中的学习课程
我也不知道叫什么名字好进入https://nodeschool.io/zh-cn/，你能够看到许多课程。这些课程大多面向零基础的，非常有趣（类似于许多关卡的小游戏）。快速开始首先，确定你的电脑已经安装好Node.js，因为我们之后都将会使用到它。然后使用命令npm（Node.js 自带的工具）去安装下面的教学课程。一旦安装完成后，这些教学课程的名称都将是一个个电...
2019-09-14 19:44:00
126
转载 利用GitHub+Node.js+Hexo搭建个人博客
本篇是自己在搭建Hexo博客平台时的一个过程记录。（2019。9.13实测有效）GitHub账号注册因为此文所搭建的个人博客是基于GitHub平台服务的，所以首先是注册GitHub，当然已有账号的跳过此步；创建Repository点击开始后，跳转至Create a new repository，这边Repository name命名规则就是username.githu...
2019-09-13 21:36:00
159
转载 更丰富的符号工具包 Font Awesome
我时常想要在此类文档中通过一些图形符号来表达更丰富的含义或是对段落进行标注，例如使用 Emoji。然而 Emoji 在这方面仍然有存在一些不足，如：颜色与文字风格不统一，在不同系统的平台上显示不统一，种类有限，很多时候找不到合适的表情。这时，可以在文档写作时引入另一个更丰富的符号工具包Font Awesome。简介Font Awesome...
2019-09-13 17:57:00
290
转载 Markdwon入门2
插入表情这里是指广义的表情包，包括表情、物体、动物等。:+1::smile: :s :scream: :kissing_heart: :yum: :cry: :blush::frog: :cow: :pig::computer: :pill: :house::one: :two: :three: :four:更多的表情用法：ht...
2019-09-13 16:29:00
98
转载 Codechef：Fibonacci Number/FN——求通项+二次剩余+bsgs
题意定义 $F_n$ 为$$F_n = \left\{\begin{matrix}0, n=0\\ 1, n=1 \\F_{n-1} + F_{n-2}, n > 1\end{matrix}\right.$$现给你一个素数 $p$ 和一个非负整数 $C$，你需要最小的非负整数 $n$，使得 $F_n \equiv C (mod \ p)$.分析因为题目保证 $p...
2019-09-12 20:37:00
130
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人