- 博客(28)
- 收藏
- 关注
RSS订阅
原创 研学-元胞自动机模拟疫情传播
蒙特卡罗方法模拟疫情传播论文提要2020年伊始，新冠肺炎席卷全球。长时间潜伏期、极高的传染率，使得全人类面对这次疫情显得有些措手不及。为了更好地认识、防控传染病，本文采用蒙特卡罗方法中的一种——元胞自动机，借鉴著名数学家约翰·康威所发明的生命游戏规则，自创规则进行疫情模拟。最终发现严格的防控、警惕的意识与强大的医疗水平均是抗击疫情必不可少的内容。政府越早采取防控措施防控效果越好；并且医疗水平的支持更是防控的关键。主题词蒙特卡罗方法；元胞自动机；疫情传播正文1 研究背景1.1 传染病介绍在人
2020-08-11 10:28:38
7836
16
原创 多项式公式处理
目的用来处理多项式公式：对于一个仅含有自变量x的多项式，将其按次数展开。大致思路一、拆括号二、读取每个无括号的表达式三、记录结果四、以后再遇到这个括号直接读取结果具体内容在实现过程中，以下顺序是按照我们一开始的思考顺序，又简入深polystruct poly{ int dag; int pre_op;//多项式前面的符号 1234代表 +*/^ pair<int,int>c[N];//c[a]分数表达c^a的系数 int pre,aft;//指针 前后的多项式 i
2021-01-30 13:28:46
917
原创 公式处理-数字四则运算不含括号
首先发现，计算机的精度有限，导致1/3+1/3+1/3！=11/3+1/3+1/3！=11/3+1/3+1/3！=1所以就把1看成1/10.5看成1/2 每个运算都是分数的加减法字符串转数字形式将字符串转化为数字形式第一行：共几个数字第二行开始：每行描述一个数字（包括该数字前面的运算符）x y zx表示符号 123代表+*/ 然后x的正负代表数字正负然后用y/z描述该数字大小#include<bits/stdc++.h>//不能乘负数 using namespace std
2021-01-25 18:54:08
246
1
原创 森林火灾模拟
森林火灾模拟代码import numpy as npimport matplotlib.pyplot as pltimport randomimport copyclass Forest: def __init__(self): self.forest_size = (100,100) self.forest = np.zeros(self.forest_size,int) self.forest_temp = copy.deepcopy(
2020-11-23 19:06:14
1489
原创 CSP2020游记
CSP/J初三开始学信竞，第一次普及因经验不足错失省一。今年高二的我为了防止考提高组的上午没事情干，报名了个普及来玩玩。T1非常简单的位运算拆分二进制，判断奇偶性然后拆分即可，五分钟打完过大样例。#include<bits/stdc++.h>using namespace std;int read(){ char s; int x=0,f=1; s=getchar(); while(s<'0'||s>'9'){ if(s=='-')f=-1; s=get
2020-11-14 14:58:11
217
原创 小程序基本语句01
app.json易错点1.忘记加，2.“ “用英文字符”:“代替 ”=“作用-特殊字段：pages字段：1.新建页面+确定页面出现顺序"pages":[ "pages/mine02/mine02", "pages/index/index", "pages/logs/logs", "pages/mine/mine" ]window字段：2.设置初始页面属性置顶区域"window":{ "backgroundTextStyle":"d
2020-07-22 14:43:58
387
原创 P2486 [SDOI2011]染色
题目题目链接题目大意：一棵树，每个点权是一个颜色。支持两种操作：1.修改任意一条链，使得这条链上的颜色均为c；2.询问任意一条链上的颜色段数。例如：11221 颜色段数为3分析大体思路由于是链上操作，考虑树链剖分+线段树的做法这两个操作都是区间修改，1操作直接普通的区间修改+dag标记即可实现线段树合并细节处理下面重点来讲讲2：线段树的查询中，假设查询的区间为[l,r][l,r][l,r]定义两个全局变量lstl、latrlstl、latrlstl、latr表示col[l]、col
2020-07-16 16:35:32
154
原创 P3730曼哈顿交易-莫队算法
题目题目链接题目大意：一个区间内有n个数字，有m个询问，每次询问一个区间。回答这个区间的第k大的热度值。 热度值定义：数字i的热度值为：这个区间内i的出现次数。分析正解前置知识权值分块；莫队算法首先暴力做法：对于每一个区间，暴力搜索出每个数字的出现次数，把出现次数push到一个数组里面，寻找到出现次数的第k大的数字。优化一：权值分块优化寻找第k大数字，可以采用权值分块的方法。一共n个数字，因此热度值的取值范围是[0,n][0,n][0,n]，所以把这个区间进行分块，维护每个块中
2020-07-13 15:16:12
208
原创 CF489E：hiking-01分数规划
题目题目链接题目大意一个数轴上有n个点，每个点的坐标均为正整数，每个点有一个美丽值。给定常数p，从原点出发，只能向正方向移动，每次必须停在一个点上，假设经过的距离为d，那么这次hike的疲劳值为abs(d−p)\sqrt{abs(d-p)}abs(d−p)。现在要从原点前往最左端的点，求总疲劳值/总美丽值的最小值。总美丽值为停留的点的美丽值之和。题目分析分数规划注意到求两个和式的比值，并且以一个点为终点就要对疲劳值、美丽值都产生贡献。所以可以考虑01分数规划1.二分一个答案，然后令新点
2020-06-15 19:00:51
258
原创 01分数规划学习笔记
何为01分数规划01分数规划解决的是这么一类问题：从给定的数列ai、bi{a_i}、{b_i}ai、bi，求∑i=1ai∗pi∑i=1bi∗pi\frac{\sum_{i=1}{a_i}*p_i}{\sum_{i=1}{b_i}*p_i}∑i=1bi∗pi∑i=1ai∗pi的最值问题其中p为一个01数列。也就是从一个二元组列表中选一些二元组，求x项之和/y项和的最值。常常是求最大中的最小、最小中的最大。这类问题我们通常用二分+分数规划来做详解直接求解不容易，我们可以考虑如何
2020-06-14 19:35:21
191
原创 P1730 最小密度路径-01分数规划
题目大意题目链接题目大意：一个有向无环图（DAG），边权均为正整数。定义一条路径的密度为 路径上的边权和/边数。先有q组询问，每次给出xy，求从x到y的若干条路径中的最小路径密度是多少。题目分析令c为边权，d为路径条数那么所求即为：∑cd\frac{\sum{c}}{d}d∑c这是一个和式分式，对于这种题目，我们通常用一种方法叫做 01分数规划我们二分答案，然后判断。判断过程：把每条边边权-二分答案然后判断最短路的正负性。如果最短路是负的那么答案不合法，修改r否则修改l代码下附
2020-06-14 19:32:54
234
原创 HDU5201-The Monkey King
题目题目大意简化题意：1~n个不同盒子，放入m个相同小球，盒子可空，但是1号盒子的球数必须严格最多，问合法方案数。分析简化问题如果没有1号盒子严格最多的限制，就变成了一个插板的问题。也就是x1+x2……+xn=mx_1+x_2……+x_n=mx1+x2……+xn=m这类不定方程求非负整数解的问题，插板法解决。容斥如果枚举猴王（1号盒子）的球数为i那么问题就变成不定方程并且0<=x<=i0<=x<=i0<=x<=i，加了一个上限。上限通常会用容斥转
2020-05-31 21:03:45
208
原创 Kings Colors-二项式反演
题目题目链接题目大意：给定一个n个节点的树，给它染色并且使得相邻节点异色。问恰好用k种颜色的染色方案数题目分析恰好k种不是很好求，因为我们很难保证每种颜色都用到，于是我们先考虑求最多k种颜色。那么就让每个点和它的父亲节点异色就可以了。也就是k∗(k−1)n−1k*(k-1)^{n-1}k∗(k−1)n−1那么我们令f(i)f(i)f(i)表示最多用i种颜色的染色方案数 g(i)g(i)g(i)表示恰好用i种的染色方案数那么假设我们要最多用n种，那么不妨从中随意选其中i种，用且仅用
2020-05-30 15:53:51
202
原创 洛谷P4859已经没有什么好害怕的了-二项式反演
题目题目链接题目大意：有两组长度为n的数ab，组间数字互不相同。两组数字之间两两配对，不重不漏，给定一个数字k，求满足a中数字>b 比 b中数字>a 的组数多k的分组方案数题目分析1.问题转化因为组间数字两两不同，所以其实就是告诉你了有多少组a比b大。令新k=n+k2\frac{n+k}{2}2n+k，那么如果n+k是一个奇数答案为0.现在就变成了询问a>ba>ba>b的数目为k的方案数2.动态规划首先对ab数组排序直接求正好为k的不太好求，因为我们考虑从
2020-05-30 14:05:50
285
原创 二项式反演-学习笔记
定义如果两个函数g，f满足：g(n)=∑i=1nCnif(i)g(n)=\sum_{i=1}^{n}C_{n}^{i}f(i)g(n)=∑i=1nCnif(i)那么f(n)=∑i=1n(−1)n−iCnig(i)f(n)=\sum_{i=1}^{n}{(-1)^{n-i}C_{n}^{i}g(i)}f(n)=∑i=1n(−1)n−iCnig(i)证明以上推下为例：1.将g函数表达式代入：f(n)=∑i=1n(−1)n−i∗Cni∗[∑j=1iCij∗f(i)]f(n)=\sum_{i=1
2020-05-30 13:22:36
309
原创 莫比乌斯反演-能量采集题解
题目题目链接重点：在下面的解答中默认n<=mn<=mn<=m前置知识莫比乌斯函数与欧拉函数莫比乌斯反演熟练地和式运算整除分块迪利克雷卷积分析步骤一：转化答案很容易想起这道题：仪仗队从（0，0）（0，0）（0，0）到x，y）x，y）x，y）之间经过的点数=gcd（x，y）−1=gcd（x，y）-1=gcd（x，y）−1为什么呢？因为xy同时乘除某一个数后...
2020-05-05 19:20:16
162
原创 数表题解
题目题目链接初步转化令f(n)f(n)f(n)表示n的所有因数和如果d∣i并且d∣jd|i 并且 d|jd∣i并且d∣j，那么d∣gcd(i,j)d|gcd(i,j)d∣gcd(i,j)那么对于每一个询问，答案为：∑i=1n∑j=1m[f(gcd(i,j))<=a]∗f(gcd(i,j))\sum_{i=1}^{n}\sum_{j=1}^{m}{[f(gcd(i,j))<=a...
2020-05-05 19:18:41
201
原创 莫比乌斯反演-数字表格题解
题目题目链接分析初步变形这道题其实是求：∏i=1n∏j=1mfibgcd(i,j)\prod_{i=1}^{n}\prod_{j=1}^{m}{fib_{gcd(i,j)}}∏i=1n∏j=1mfibgcd(i,j)和别的莫比乌斯反演题目不同之处在于：它是累乘不是累加但是转念一想，底数累乘实际上就是指数的累加，所以我们可以这样做：=∏d=1nfibd∑d∣in∑d∣jm[gc...
2020-05-05 19:16:22
199
原创 [HNOI2005]狡猾的商人：一题两解
题目链接题目大意：一段区间，给定m个信息，每个信息表示一段子区间中的区间和，问你这个区间是否合法（也就是说是否存在这种区间）第一种解法：差分约束：由于是一段区间和，那我们不妨把它看成是前缀和形式的两个相减。然后将它用差分约束的方式转换成一张图，那么这个区间是否合法就变成了这张图上是否有最短路，就变成了用SPFA判负环的问题。1.建图：对于一条信息：at−as==x(t>=...
2020-01-29 10:13:27
131
原创 apple tree POJ2486 树形dp
题目链接http://poj.org/problem?id=2486题目大意：给你一颗树，树上每个点有一个点权，让你从1号点走k步，问可以得到的最大点权和是多少？思路：这道题既然是个树结构，我们很容易想到是树形dp设计状态：由于我们走到一个点有可能是要回到这个点走到另一颗子树上去，也有可能不回来。所以要设两个状态：1.dp[0][x][s]dp[0][x][s]dp[0][x][s]表示...
2019-08-30 20:24:27
105
原创 八数码问题-启发式搜索
八数码想必大家小时候都玩过，今天我们就一起来解决这道题目，如何输出八数码的答案。今天先介绍启发式搜索，也就是A*，后续我会继续推出双向广搜与迭代加深启发式搜索（IDA*）。题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1043首先，...
2019-08-13 21:09:16
1232
原创 Examining the Rooms 组合数学
题目：Examining the Rooms 题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3625题目大意：给出N个房间，每个房间的钥匙随机放在某个房间内，概率相同。有K次炸门的机会，求能进入所有房间的可能性为多大。但是1号房间不能打开我们可以先打开k个房间，也就是获得k把钥匙，我们要用这些钥匙打开所有房间。因为要求概率，所以可以用古典概型。...
2019-08-04 21:43:58
370
原创 组合数学学习笔记
错位排序n个物品，每个物品都不在对应的位置上，即对于任意的i∈[1,n],均有第i个物品不在i位置上的方案数。令dp[i]表示前i个物品的错排方案数，递推解决问题。对于第n个物品，我们假设dp[1]------dp[n-1]都已经处理好了，那么这个物品会有两种情况：1.前n-1个物品已经错排好了，那么这个物品只用和前面任意一个物品交换位置即可，即dp[n]=(n-1)*dp[n-1]2....
2019-08-04 21:26:35
269
1
原创 New Year Tree 线段树+dfs序+状态压缩
知识储备：线段树、图的存储与遍历、状态压缩题意翻译你有一棵以1为根的有根树，有n个点，每个节点初始有一个颜色c[i]。有两种操作：1 v c 将以v为根的子树中所有点颜色更改为c2 v 查询以v为根的子树中的节点有多少种不同的颜色题目链接：https://www.luogu.org/problem/CF620E首先我们讲一下dfs序在这道题里面的应用定义：一棵树从根节点开始进行深度...
2019-08-03 21:55:32
291
1
原创 树上问题：基础版T2：Ilya And The Tree
树上问题种类非常多，今天的主要讲述的是以搜索为框架，加一些树上优化的题目。先定义一些通用的符号：size[i]:以i号结点为根的子树的结点个数。dist[i]:i号结点到根节点的距离（无边权情况下默认为一）p[i][j]：只有最后一题出现了，用于树上倍增，表示第i个点向上的第2^j点。dep[i]一个点i的深度（也就是无边权情况下下的dist)上期文章 T1：Kay and Snowflak...
2019-08-02 21:42:30
126
原创 树上问题：基础版T1：Kay and Snowflake
树上问题种类非常多，今天的主要讲述的是以搜索为框架，加一些树上优化的题目。先定义一些通用的符号：size[i]:以i号结点为根的子树的结点个数。dist[i]:i号结点到根节点的距离（无边权情况下默认为一）p[i][j]：只有最后一题出现了，用于树上倍增，表示第i个点向上的第2^j点。dep[i]一个点i的深度（也就是无边权情况下下的dist)我们来看第一道题::题目链接：https://...
2019-08-02 21:25:59
197
原创 HISTOGRA - Largest Rectangle in a Histogram
题目：如图所示，在一条水平线上有n个宽为1的矩形，求包含于这些矩形的最大子矩形面积（图中的阴影部分的面积即所求答案）。输入格式：有多组测试数据，每组数据占一行。输入零时读入结束。每行开头为一个数字n（1<=n<=100000）,接下来在同一行给出n个数字h1h2…hn（0<=hi<=1000000000）表示每个矩形的高度。输出格式：对于每组数据，输出最大子矩阵...
2019-07-31 14:20:46
240
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人