- 博客(121)
- 收藏
- 关注
RSS订阅原创 pytorch中的广播机制
广播机制原理广播机制有什么用？广播机制实际上是在运算过程中，去处理两个形状不同向量的一种手段。什么情况下会出现广播机制？如果相加（数学运算）的两个数组的shape不同, 就会触发广播机制。广播机制的规则（条件）是怎样的？1> 两个张量都至少有一个维度# 反例x=torch.empty((0,))y=torch.empty(2,2)2> 按从右往左顺序看两个张量的每一个维度，x和y每个对应着的两个维度都需要能够匹配上匹配含义：a. 这两个维度的大小相
2023-04-09 23:05:00
314
原创 vs2019添加库文件
环境：windows10，64位，vs2019vs安装路径：D:\visual_studio需要的材料：头文件，编译好的库文件。
2022-09-15 22:20:07
2959
原创 Cmake小总结
CMake是一个跨平台的编译(Build)工具,可以用简单的语句来描述所有平台的编译过程。CMake能够输出各种各样的makefile或者project文件，能测试编译器所支持的C++特性,类似UNIX下的automake。CMake 不仅可以编译源代码、制作程序库、产生适配器（wrapper）、还可以用任意的顺序建构执行档。
2022-09-15 21:46:37
2005
原创 GMM基础
高斯混合模型就是用高斯概率密度函数(正态分布曲线)精确地量化事物，它是一个将事物分解为若干的基于高斯概率密度函数(正态分布曲线)形成的模型。假如我们我们现有的数据分布如红线所示，可以发现用一个高斯分布很难较好的描述这组数据的分布，所以我们能够用两个高斯分布加权平均得到一个新的分布，其中每一个数据点都有一定概率属于两个高斯分布中的一个。GMM只是若干个高斯分布的加权平均（上图中是两个高斯分布）P(x)=∑k=1KαkN(μk,Σk)∑k=1Kαk=1P(x)=\sum_{k=1}^K\alpha_kN(\
2022-08-19 20:32:52
816
原创 Calculus of Variations：变分计算
functional）指以函数构成的向量空间为定义域，实数为值域为的“函数”，即某一个依赖于其它一个或者几个函数确定其值的量，往往被称为“函数的函数”。函数代表了数到数的映射，而泛函代表了函数到数的映射，即给定一个函数，泛函能够得到一个数。当然也满足边界条件，也代表两点之间的一条路径，并且在确定的。不等于0，被积函数就是任意的，无法保证积分恒为0。下，其函数表达式只与ϵ有关。给定后是一个固定的函数。,我们首先来计算在任意确定的。当泛函取极值时，函数。为了要满足上面的条件。...
2022-07-31 16:56:24
501
原创 图像梯度(opencv-c++)
梯度的本意是一个向量（矢量），表示某一函数在该点处的方向导数沿着该方向取得最大值，即函数在该点处沿着该方向（此梯度的方向）变化最快，变化率最大（为该梯度的模）（拉普拉斯）算子是一种二阶导数算子，其具有旋转不变性，可以满足不同方向的图像边缘锐化（边缘检测）的要求。拉普拉斯算子（LaplaceOperator）是n维欧几里德空间中的一个二阶微分算子，定义为梯度（▽f）的散度（▽·f）这分别是图像在(x,y)点处x方向和y方向上的梯度，从上面的表达式可以看出来，图像的梯度相当于2个相邻像素之间的差值。...
2022-07-31 00:05:54
1821
原创 Codeforces Round #779 (Div. 2) D2. 388535 (Hard Version) （异或，01字典树）
D2. 388535 (Hard Version)题意：给出l，rl，rl，r，和一个序列aaa，aia_iai由l，rl，rl，r中的一个数异或xxx得到，问xxx可以取什么值。思路：l，rl，rl，r中的每个元素互不相同，我们只要确定一个xxx，min{x⨁a[i] ∣ i∈[1,r−l+1]}=lmin\{x \bigoplus a[i] ~|~i \in[1, r-l+1] \}=lmin{x⨁a[i] ∣ i∈[1,r−l+1]}=l，max{x⨁a[
2022-04-04 22:09:56
540
原创 pip 问题 记录一下
pip --default-timeout=100 install 库名称 -i http://pypi.douban.com/simple/ --trusted-host pypi.douban.com清华：https://pypi.tuna.tsinghua.edu.cn/simple阿里云：http://mirrors.aliyun.com/pypi/simple/中国科技大学 https://pypi.mirrors.ustc.edu.cn/simple/华中科技大学：http://pypi
2021-12-15 20:08:25
189
1
原创 Codeforces Round #752 (Div. 2)(A-E)
Codeforces Round #752 (Div. 2)A. Era题意：给出一个数组aaa，一种操作：每次能够在任意位置插入一个数，问最小的插入次数使∀i∈[1,n]\forall i \in[1,n]∀i∈[1,n]，都满足ai≤ia_i\le iai≤i思路：我们能够在最前面插入1，使每个数的下标都加1，如果i≤a[i]i\le a[i]i≤a[i]，只要让最小i−a[i]i-a[i]i−a[i]变成0的操作数就是答案。#include <bits/stdc++.h>
2021-11-01 21:17:30
293
1
原创 2018-2019 ACM-ICPC, Asia Nanjing Regional Contest（J Prime Game）
J Prime Game题意：给出一个序列，mul(l,r)=∏i=lraimul(l,r)=\prod\limits_{i=l}^ra_imul(l,r)=i=l∏rai，fac(l,r)fac(l,r)fac(l,r)表示mul(l,r)mul(l,r)mul(l,r)的中不同的质数的个数。问：∑i=1n∑j=infac(i,j)\sum\limits_{i=1}^n\sum\limits_{j=i}^nfac(i,j)i=1∑nj=i∑nfac(i,j)思路：统计在序列中出现过的质
2021-10-16 16:15:25
151
原创 Codeforces Round #746 (Div. 2)（A,B,C）
Codeforces_Round#746(Div.2)A.Gamer_Hemose题意：思路：代码：B.Hemose_Shopping题意：思路：代码：C.Bakry_and_Partitioning题意：思路：代码：Codeforces_Round#746(Div.2)A.Gamer_Hemose题意：一个 BossBossBoss 的血量为mmm，有nnn把武器，每次能够使用一把武器，每把武器不能连续使用，问最少要几次才能杀掉BossBossBoss思路：选两个伤害最大的武器，交替使用。.
2021-10-07 18:36:58
549
1
原创 Rinne Loves Sequence
Rinne Loves Sequence推式子∑i=1n∑j=i+1n[gcd(ai,aj)=1]∑i=1n∑j=1n[gcd(i,j)=1]×cnti×cntj−∑i=1n[gcd(i,i)=1]∑d=1nμ(d)∑i=1n∑j=1ncnti×cntj∑d=1nμ(d)∑i=1nd∑j=1ndcntid×cntjd\sum_{i=1}^n\sum_{j=i+1}^n[gcd(a_i, a_j)=1]\\\sum_{i=1}^n\sum_{j=1}^n[gcd(i,j)=1]\times cnt_
2021-09-21 10:32:44
140
原创 逆序对数列（简单dp）
P2513 [HAOI2009]逆序对数列题意：给出一个长度为nnn的排列，问有多少种排列方式的逆序对为k。思路：dp[i][j]dp[i][j]dp[i][j]表示长度为i的排列，逆序对为jjj的方案数。对于一个长度为iii的排列，将i+1i+1i+1插入到排列中有i+1i+1i+1种插入方法，最多能够多出iii个逆序对，所以可以得到转移方程dp[i][j]=∑k=0min(i−1,i(i−1)2)dp[i−1][j−k]dp[i][j]=\sum_{k=0}^{min(i-1, \fra
2021-09-13 19:00:31
327
原创 P1681 最大正方形II（简单dp）
P1681 最大正方形II题意给出一个n×mn\times mn×m的010101矩阵，问最大的01交错的正方形大小是多少。思路：dp[i][j]dp[i][j]dp[i][j]表示以(i,j)(i,j)(i,j)这个位置为右下角的最大的矩形的大小。考虑如何转移。当(i,j)(i,j)(i,j)和(i−1,j),(i,j−1)(i-1,j),(i,j-1)(i−1,j),(i,j−1)不相同时，(i,j)(i,j)(i,j)这个点才能向上左扩展。否则dp[i][j]=1dp[i][j]=1
2021-09-12 20:05:58
190
原创 EOJ Monthly 2021.9 Sponsored by TuSimple——A.Amazing.Discovery（分治or二次剩余）
EOJ Monthly 2021.9 Sponsored by TuSimpleA.Amazing.Discovery题意：思路：解法1：解法2：EOJ Monthly 2021.9 Sponsored by TuSimpleA.Amazing.Discovery题意：给出a，b，na，b，na，b，n，求S=(a+b)n+(a−b)n(mod p)S=(a+\sqrt b)^n+(a-\sqrt b)^n(mod~p)S=(a+b)n+(a−b)n(mod p)思.
2021-09-11 15:56:46
353
原创 二次剩余（简单理解）
二次剩余定义ppp是一个奇素数。x2≡n(mod p)x^2\equiv n(mod~p)x2≡n(mod p)判断n是否有解勒让德符号 ≡\equiv≡ 欧拉判别准则。勒让德符号(np)={1p∤n ∧ n是p的二次剩余−1p∤n ∧ n不是p的二次剩余0p∣n\left(\frac {n}{p}\right)=\begin{cases}1& p\nmid n~\wedge~n是p的二次剩余\\-1&p\nmi
2021-09-11 15:35:36
1057
原创 P2359 三素数数（dp)
P2359 三素数数题意：如果一个数的所有连续三位数字都是大于100100100的素数，则该数称为三素数数。比如113797113797113797是一个666位的三素数数。给出一个nnn，问有多少长度为n的三素数数，对1e9+91e9+91e9+9取模。思路：dp[i][j][k]dp[i][j][k]dp[i][j][k]表示长度为iii的，倒数第二位是jjj，最后一位是kkk的方案数。这样我们dp[i][j][k]dp[i][j][k]dp[i][j][k]怎么从dp[i−1][k][
2021-09-07 21:43:22
328
原创 P2401 不等数列（dp)
P2401 不等数列题意：给出一个长度为nnn的排列，每种排列方式在两个数之间都能产生一个大于号小于号，问有k个小于号的排列方式有多少个思路：现在是长度问n−1n-1n−1的一个排列，我们把n插入到这个排列中就成了一个长度问n的排列。如果我们把nnn放在序列的最前面，就会多一个大于号。如果我们把nnn放在序列的最后面，就会多一个小于号。如果我们放在小于号的中间，我们就能减掉一个小于号，然后再加上一个大于号和一个小于号，也就是加上一个大于号如果我么放在大于号的中间，我们就能减掉一个小于号
2021-09-05 16:53:58
191
1
原创 P1622 释放囚犯（区间dp）
P1622释放囚犯题意：思路：CodeCodeCodeP1622释放囚犯题意：p个牢房，q个要释放的人，两个监狱能够通信当且仅当两个监狱之间都有人，当释放某个监狱的人时，和他能够通信的监狱必须要吃肉，问要按什么样的顺序释放才能让所有人吃上的肉最少思路：区间dp，反过来思考，释放人，我们考虑先把要放的人全部清空，然后再抓进去（可以想到，最优情况的顺序下，正向和反向的答案相同）这样一来，我们每次抓进去一个人，就合并两个区间，且加上区间人数。dp[i][j]dp[i][j]dp[i][j]表示把第.
2021-09-04 21:50:17
159
原创 E. e-Government（AC自动机，fail树，树状数组）
E. e-Government题意：给出一些m个模式串和n个操作，操作分为3种：‘?’+字符串：查询字符串中 在模式串集合中的字符串 出现了多少次。‘+’+数字：将第几个模式串添加回模式串集合‘-’+数字：将第几个模式串从模式串集合中删除。思路：考虑暴力的做法：先对所有模式串建trie树，bfs找fail数组，用一个boo数组来表示某个字符串是否在模式串集合中。对于一个查询，字符串在trie树上跳fail，当跳到还在模式串集合中的节点时，答案++，输出答案。树状数组维护差分数组做法：
2021-09-03 15:25:14
215
1
原创 P3966 [TJOI2013]单词（AC自动机，fail树）
P3966 [TJOI2013]单词题意：给出n个字符串，问各个字符串在所有字符串中出现的次数思路：对于一个字符串来说，我们对它在trie树上的路径节点++，表示当前节点一共出现多少次。再建一棵fail树，对于fail树的一个子树来说，子树权值等于根节点所表示的字符串出现的次数。所以我们就可以dfs去统计每个子树的权值和，在输出每个字符串对应节点的权值，就是答案。#include<bits/stdc++.h>using namespace std;#define int lon
2021-09-03 14:21:52
97
原创 2021“MINIEYE杯”中国大学生算法设计超级联赛（9）——H - Integers Have Friends 2.0
[H - Integers Have Friends 2.0](https://vjudge.net/problem/HDU-7073)题意：思路：H - Integers Have Friends 2.0题意：给出一个数组aaa，找到一个最长的子序列bbb，存在m(m>1)m(m>1)m(m>1)使bi mod m=bj mod mb_i~mod~m=b_j~mod~mbi mod m=bj mod&nbs.
2021-08-31 17:33:57
165
原创 P2414 [NOI2011] 阿狸的打字机（AC自动机，fail树）
[P2414 [NOI2011] 阿狸的打字机](https://www.luogu.com.cn/problem/P2414)题意：思路：解释一些地方代码P2414 [NOI2011] 阿狸的打字机题意：在一串长字符串中，可以构建出多个字符串，n组询问，问第x个字符串在第y个字符串中出现的次数。思路：先把每个字符串插入到trie树上（注意插入的方法）再对每个询问按照y进行归类，相同的y则放到同一个数组中，方便得到答案。再求出每个节点的fail，建出一棵fail树，用dfs找出fail树的d.
2021-08-31 16:07:15
125
原创 2021牛客暑期多校训练营2——Product of GCDs
Product of GCDs问题：求给定的多重集合的 所有长度为k的子集 的gcdgcdgcd的 乘积。mx=max{multiple_set}ans=∏i=1mxikiki表示在集合中的长度为k且gcd为i的子集的个数。先求出i的倍数个数num，再求出(numk),有因为其中包含gcd为2i的集合，所以要减到这些集合产生的贡献。(容斥一下)\begin{aligned}mx&=max\{multiple\_set\}\\ans&=\prod_{i=1}^{mx}i^{k_i}\
2021-08-29 15:49:31
126
原创 P6217 简单数论题
题意：给出一个长度为 n 的序列 a，q次询问 ∏i=lrlcm(ai,x)\prod\limits_{i=l}^r \operatorname{lcm}(a_i,x)i=l∏rlcm(ai,x) 的值。答案对 109+710 ^ 9 + 7109+7取模。思路：∏i=lrlcm(ai,x)=∏i=lraixgcd(ai,x)=∏i=lraix∏i=lrgcd(ai,x)\begin{aligned}&\prod_{{i=l}}^rlcm(a_i,x)\\=&\prod.
2021-08-29 15:15:55
267
转载 隔板法（插板法）
转载自： 原文链接 理解隔板法【定义】【公式】【隔板应用】普通隔板法添元素隔板法添板插板法选板法分类插板逐步插板法理解隔板法【定义】隔板法就是在nnn个元素间的（n−1）（n-1）（n−1）个空中插入kkk个板，可以把nnn个元素分成k+1k+1k+1组的方法。应用隔板法必须满足3个条件： （1） 这nnn个元素必须互不相异；（2） 所分成的每一组至少分得1个元素；（3） 分成的组别彼此相异。【公式】把10个相同的小球放入3个不同的箱子，每个箱子至少一个，问有几种情况？Cn−1m−
2021-07-14 21:10:34
7071
原创 P3768简单的数学题(莫比乌斯，欧拉函数，杜教筛）
P3768简单的数学题解法一：推式子部分杜教筛部分：解法二：推式子部分杜教筛部分：P3768简单的数学题解法一：推式子部分∑i=1n∑j=1nijgcd(i,j)∑i=1ni∑j=1nj∑d∣i,d∣jd[gcd(i,j)=d]∑d=1nd3sumi=1ndi∑j=1ndj[gcd(i,j)=1]∑d=1nd3∑i=1ndi∑j=1ndj∑x∣gcd(i,j)μ(x)∑d=1nd3∑x=1ndx2μ(x)∑i=1ndxi∑j=1ndxj∑d=1nd3∑x=1ndx2μ(x)(⌊ndx⌋(1+⌊nd.
2021-07-13 23:51:52
130
1
原创 C. Need for Pink Slips
C. Need for Pink Slips问题：一个游戏，给出C,M,PC,M,PC,M,P三个事件，每个事件的概率是c，m，pc，m，pc，m，p，当出现事件PPP的时候，游戏结束。在每一轮游戏过程中c，m，pc，m，pc，m，p都会根据规则发生变化。问游戏进行多少轮的期望。规则：给定一个小数vvv。(c,m,p)=(0.2,0.1,0.7)(c,m,p)=(0.2,0.1,0.7)(c,m,p)=(0.2,0.1,0.7)， v=0.1v=0.1v=0.1， 在发生事件C后，因为c>v
2021-07-08 21:27:09
616
3
原创 E. Number Challenge(清晰地推式子)
E. Number Challenge问题：∑i=1a∑j=1b∑k=1cd(ijk)\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk)∑i=1a∑j=1b∑k=1cd(ijk)，d(n)d(n)d(n)表示n的因子个数。推式子：∑i=1a∑j=1b∑k=1cd(ijk)∑i=1a∑j=1b∑k=1c∑x∣i∑y∣j∑z∣k[gcd(x,y)=1][gcd(y,z)=1][gcd(x,z)=1]∑x=1a∑y=1b∑z=1c∑i=1a∑j=1b∑k=1c∑x∣
2021-05-19 23:31:53
163
原创 Happy 2004(代码+积性函数扩展)
Happy_2004题意：思路（详细解释在扩展）：CodeCodeCode：扩展（因子和，积性函数）nnn的全部因子之和是：积性函数σ\sigmaσHappy_2004题意：求2004k2004^k2004k的全部因子的和，对292929取模思路（详细解释在扩展）：积性函数，等比数列求和CodeCodeCode：#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N = 1e3+.
2021-05-11 21:58:48
365
原创 原根（模板+证明）
原根定义:原根判定定理原根的个数原根的存在定理题目：原根定义:设m∈N∗m\in\mathbb{N^*}m∈N∗,a∈Za\in\mathbb{Z}a∈Z.若(a,m)=1(a,m)=1(a,m)=1,且δm(a)=ϕ(m)\delta_m(a)=\phi(m)δm(a)=ϕ(m),则称aaa为模mmm的原根（注：a，ma，ma，m互质）.原根判定定理若ggg是模mmm的一个原根,则对于ϕ(m)\phi(m)ϕ(m)的任何大于111且不为自身的因数ppp,都有gϕ(m)p≡1(mod .
2021-05-11 00:13:40
1537
2
原创 F. Lunar New Year and a Recursive Sequence（矩阵快速幂+BSGS）
F. Lunar New Year and a Recursive Sequence题意：给出f1=f2=⋯=fk−1=1f_1=f_2=\cdots=f_{k-1}=1f1=f2=⋯=fk−1=1和b1,b2⋯bkb_1,b_2\cdots b_kb1,b2⋯bk，还有递推方程fi=fi−1b1fi−2b2⋯fi−kbkf_i=f_{i-1}^{b_1}f_{i-2}^{b_2}\cdots f_{i-k}^{b_k}\\fi=fi−1b1fi−2b2⋯fi−kbk问
2021-05-08 20:59:09
1363
2
原创 2020 ICPC Universidad Nacional de Colombia Programming Contest---A
A. Approach题意：给出444个点，A，B，C，DA，B，C，DA，B，C，D，两个人a,ba,ba,b，a从AAA到BBB，b从CCC到DDD，他们两个速度相同，问他们相距最近的距离大小（注：当某一个人到达终点时，他会停留在终点，而另一个人仍然会继续前进）思路：猜测：距离的他们两点的距离是一个关于时间的凹函数，然后就是三分去查找。注意细节：当一个人已经到达终点时，可以肯定当前距离函数仍然是一个凹函数或者一个单调递增的函数，所以我们要再用一次三分，来找到当前段的最小值，再和上一段进行比较。
2021-05-07 11:37:40
196
原创 2020 ICPC Universidad Nacional de Colombia Programming Contest---D
D.Dice题意：思路：Code：D.Dice题意：给定nnn个色子，每个色子有kkk个面，每个面都有一个数字，数字模mmm为000的面出现的概率为000，剩下的面出现的的概率相等，每个面在模mmm下，只有mmm中情况，我们可以算出一个色子的mmm个概率。思路：设aia_iai表示111个色子丢出的面模mmm是iii的概率。设bib_ibi表示111个色子丢出的面模mmm是iii的概率。则：{b0=a0∗a0+a1∗am−1+a2∗am−2...am−1∗a1b1=a0∗a1+a1∗a.
2021-05-06 21:29:54
145
原创 BSGS(扩展篇，思路+详解)
扩展篇问题：方法:题目：扩展篇问题：当p，ap，ap，a可能不互质时，求ax≡b(mod p)a^x\equiv b(mod~p)\\ax≡b(mod p)方法:当d1=gcd(a,p)≠1d_1=gcd(a,p)\neq1d1=gcd(a,p)=1，当d1∤bd_1\nmid bd1∤b时，无解，则原式变成：ad1ax−1≡bd1(mod pd1)\frac a{d_1}a^{x-1}\equiv \frac b{d_1}(mod~\frac p.
2021-05-05 22:36:06
154
原创 BSGS(进阶篇，解法+题目)
进阶篇问题：xa≡b(mod p)x^a\equiv b(mod~p)\\xa≡b(mod p)解法一：设ggg是p的一个原根，根据原根的性质存在ccc，满足gc≡xg^c\equiv xgc≡x成立，同理存在ttt，满足gt≡bg^t\equiv bgt≡b成立。∴(gc)a≡b(mod p)(ga)c≡b(mod p)ga已知，所以我们就能用基础篇直接求解c，也就是一个特解。\therefore (g^c)^a\equiv b(mod~p)\\(
2021-05-04 20:23:08
282
原创 BSGS(基础篇，题目+详解)
基础篇问题：思路：模板：题目：基础篇问题：给出a，b，pa，b，pa，b，p，其中gcd(a,p)=1gcd(a,p)=1gcd(a,p)=1，求xxx满足ax≡b(mod p)a^x\equiv b(mod~p)\\ax≡b(mod p)思路：设x=Ap−Bx=A\sqrt p-Bx=Ap−B其中A∈[1,p],B∈[0,p]A\in[1,\sqrt p],B\in[0,\sqrt p]A∈[1,p],B∈[0,p]，得到问题的变形aAp−B≡b(mod&nb.
2021-05-03 21:47:10
2268
3
原创 C. Chef Monocarp(思维，dp)
C. Chef Monocarp问题：给出nnn道菜，每道菜有一个预期出锅时间tit_iti，但每分钟只能出一道菜，设第iii道菜的出锅时间为TiT_iTi，则会得到∣ti−Ti∣|t_i-T_i|∣ti−Ti∣的不愉悦值，问怎么安排出菜时间，才能使不愉悦值最小，输出最小值。给出一个引理：ti<tj且Ti>Tj⇒∣ti−Ti∣+∣tj−Tj∣>∣ti−Tj∣+∣tj−Ti∣t_i<t_j且T_i>T_j\\\Rightarrow\mid t_i-T_i\mid
2021-04-28 20:32:04
336
原创 阶（数论，性质+证明）
阶定义：性质：阶定义：对于m>1m>1m>1且(a,m)=1(a,m)=1(a,m)=1，使an≡1(mod m)a^n\equiv 1(mod~m)an≡1(mod m)成立的最小的nnn，称为nnn模mmm的阶，记作δm(a)\delta_m(a)δm(a)。根据欧拉定理，可以证明至少存在一个nnn使an≡1(mod p)a^n\equiv 1(mod~p)an≡1(mod p)成立。性质：a,a2,...,aδm(a)a,a^.
2021-04-27 20:16:36
3412
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人