qq_35160381的博客

原创 nlp pattern安装失败试试这个

pip install --only-binary :all: mysqlclient

原创 一行代码写鸡尾酒问题

鸡尾酒问题[W,s,v] = svd((repmat(sum(x.*x,1 ), size(x, 1),1).*x)*x’);<以后再做研究

原创 map+dijkstra+邻接表

hdu2112最短路#include<iostream>#include<cstdio>#include<algorithm>#include<vector>#include<map>using namespace std;struct Edge{ int to; int v; Edge(int t,int vv)...

原创 KL散度笔记

传送

原创 牛顿法与Hessian矩阵

好博客讲解链接

原创 deeplearningbook-note（伪逆）

对于非方矩阵而言，其逆矩阵没有定义。假设在下面的问题中，我们希望通过矩阵A的左逆B来求解线性方程： Ax=yAx=yAx = y等式两边左乘左逆B后，我们得到： x=Byx=Byx = By是否存在一个唯一的映射，将A映射到B，取决于问题的形式。 如果矩阵A的行数大于列数，那么上述方程可能没有解；如果矩阵A的行数小于列数，那么上述矩阵可能有很多解。 Moore-Penros...

原创 数学之美笔记一——文字和语言vs数字和信息

待续

原创 hdu4876(dfs取数字组合+剪枝）

传送门题意:给一个数组，任意取k个数，给定l，求能填满[l,r]的最大r./** @Author: biggf* @Date: 2018-07-05 10:42:50* @Last Modified by: biggf* @Last Modified time: 2018-07-05 13:26:38*/#include<iostream>#includ...

原创 hdu5745(字符串dp+bitset优化）

传送门题意：交换相邻字符，用二进制表示匹配。/** @Author: biggf* @Date: 2018-07-04 19:38:55* @Last Modified by: biggf* @Last Modified time: 2018-07-04 20:11:26*/#include<iostream>#include<cstring&g...

原创 后缀数组hdu5558

传送门 题意：求解所有后缀的中小标比当前后缀小的最长前缀匹配，要求lcp相等时输出下标最小的。 题解略有错误，求大佬给测试用例/** @Author: biggf* @Date: 2018-07-03 12:09:33* @Last Modified by: biggf* @Last Modified time: 2018-07-03 19:24:41*/#in...

原创 机器学习入门之dlib安装----看dlib官文档

因为一次随意的复制粘贴，找bug找了一贴。cd examplesmkdir buildcd buildcmake ..cmake --build . --config Release cmake .. 用指令 cmake -G “Visual Studio 15 2017 Win64” -T host=x64 .. 更换 如果是vs2015则改成cm...

原创 Codeforces Round #486 (Div. 3) C

传送题意：找出两个序列，各删除一个数后序列和相等，问删除哪两个序列的哪两个数 题解：预处理出删除每个数后的序列和，set存一下。#include<bits/stdc++.h>using namespace std;int k,a,sum[200005];vector<int> b[200005];set<pair<int ,int> ...

原创 Codeforces Round #486 (Div. 3) D

传送题意：给你2∗1052∗1052*10^5个数，要你找出最大的一个子集，每两个数差的绝对值是2的指数幂。题解：思维，在思维，排序，以第一个数a为基准，考虑一下与a的差为2的指数幂的数小于32个，再考虑一下这32个数中两两都要为2的指数幂时的最大个数时，似乎有点复杂，再想想，|2x−2y||2x−2y||2^x-2^y|要等于2的幂，似乎只有|x-y|=2的时候能满足，其实确实就是如此。...

原创 Avito Code Challenge 2018 C

传送 题意：给定一棵树，将树拆分成一些没有重叠的路径，要求任何两条路径之间必须要有一个交点。 题解：稍微分析，除了只有一条路径的情况，有多条路径时，要满足两两香叫，必然所有的路径相加与一点，否则会成环，那就不是树了，排除这两种情况后就是没法划分的树了。要找出这样的路径：端点度数为1，或者度数等于度数为1的顶点个数。#include<bits/stdc++.h>using ...

原创 Avito Code Challenge 2018 E

传送 题意：给定一些区间增加值的集合，要求任选若干子集之后的最大值，求出所有可能值。 题解：仔细看题目数据可以发现，它要求得值[1,n]，即可以以值枚举区间，对每个值枚举出其可能覆盖的区间，然后其它新加入的值加入这样的区间时进行组合。#include<bits/stdc++.h>using namespace std;int r[10005];struct N{ ...

原创 Avito Code Challenge 2018 D

传送 题意：给你一串数字，叫你给它划分为k段，对每段求和，求每段相&的最大值。 题解：dp[i][j]作为一个标记，1表示前j个数划分成i段满足当前要求（即当前位可为1）。#include<bits/stdc++.h>#define ll long longusing namespace std;bool dp[55][55];int n,k;int ma...

原创 2018江苏省赛G题（hdu6282）

hdu6268 题意：给你aa,bb,abab这样的三个字符串，和两个字符串s，t，让你在t中任意删除或者插入任意多个aa,bb,abab，最后得到t。 题解：思维题，简单手画画就可以得出一个结论，aa,bb,abab,可以得到任意由，a和b组成的且a和b的个数为都为偶数个的字符串。这样s可以变成t就要求c分隔后的每一段a和b同为奇或偶。上代码#include<iostream&...

原创 2018江苏省省赛c题 （hdu6278)

2018江苏省省赛 c题 （hdu6278)#include<bits/stdc++.h>#define lson l,m#define rson m+1,rusing namespace std;const int maxn=100005;struct Node { int...

原创 2018江苏省赛 F （hdu 6281）

注意精度：转换后用long long#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath&

原创 Educational Codeforces Round 44 (Rated for Div. 2) F

传送 题意：在一个字符串里面选取任意两个起始位置不同但长度相同的字符串，判断从前往后字符是否可以一对一映射。题解： 上面两个表格是分别对adcae和cebca的分析，通过hash值可以发现他们的一一映射关系，所以他们是同构的。利用这样的做法hash即可。#include<bits/stdc++.h>#define ll long longusing name...

原创 Educational Codeforces Round 44 (Rated for Div. 2)E

传送 题目大意：给定一个序列，要你把这个序列任意地分成若干个集合，如果有满足以下两个条件的集合就输出YES，否则NO，1）任意一个集合不少于k个元素，2）任意一个集合中任意两个数的差绝对值不超过d题解：先排个序再dp，dp[i]=1表示第i个数（不包括第i个数）之前的所有数满足上述条件，反之不满足。 转移方程：dp[i]=1,(若dp[j]=1且，x<=j<=i-k+1,x为序...

原创 Educational Codeforces Round 44 (Rated for Div. 2) D

传送 题目大意：给你一个数n，和H，将n分解成多个数的和，将这些数排成一个序列。但是有以下要求第一个数一定小于H，相邻两个数的差的绝对值小于1，求这样的序列最小的长度。 要使序列长度最小，很明显要使maxh尽量大，简单说明就是h越大会减少其它的拆分。因此总结出，长度最小的话有如图： 两种情况。 序列先严格增后严格减或严格减。最后如果有剩余，尽量填充到maxh高度。最后，二分ma...

原创 1-N中素数的个数（N非常大）

贴个博客以后学习

原创 数据结构与第k大

HDU2665 题意：给你一个数列长为n，（1<=n<=100000）,q次询问，(1<=q<=100000) 题解：高级数据结构 hdu2665 + 尝试一 数列分块（求小于某个数的个数），再加二分得到答案。 结果（TLE）,时间复杂度（O(q∗n−−√log2(n−−√)∗log2(num)O(q∗nlog2(n)∗log2(num)O(q*\sqrt n ...

原创 数论-素数（miller-rabin，埃式筛，欧拉筛（线性筛求欧拉函数的系数）+大数因子分解）

首先贴一个最基本的O(n−−√)n)\sqrt n)素数判定方法直接进行因子检测即可#include<bits/stdc++.h>using namespace std;int n;bool prime(int n){ if(n<2)return false; for(int i=2;i*i<=n;i++) {

原创 tarjan学习以及模板总结

tarjan无向图割点：删掉某点后的联通块增加 桥：删掉某边后的联通快增加 点双联通图：删除某点后依然联通的图（或者没有割点的图） 边双联通图：删除某边后依然联通的图（或者没有桥的图） 双联通分量：点或者边的极大双联通分量 缩点：将一个联通分量等价成一个点有向图强连通：在有向图中如果两点间至少存在一条路径，则称两个顶点强连通 强连通图：任意两点都强联通的图 强...

原创 网络流入门

贴一个网络流详解gap优化之最大流（能解决几乎所有最大流问题）简述一下：首先反向给节点标号，然后，正向跑，过程中用gap判断是否有断层，即达到优化的效果。const int maxn = 1200;const int maxe = 500000;const int inf = 1<<30;struct Edge{ int v,cap,next;}edge...

原创 manacher(最大回文串算法+手写模板）

poj3976 题意：求每个串的最大回文串长度。 这是一份通过manache算法理解过后手动模拟的算法，纯手写，便于理解#include<iostream>#include<cstdio>#include<cstring>using namespace std;char ch[2000005];int hui[2000005];char ca...

原创 cf983B dp

cf983B blemset/problem/983/B)代码：#include<bits/stdc++.h>using namespace std;int dp[5005][5005];int a[5005];int main(){ int n; scanf("%d",&n); for(int i=1;i&amp

原创 Guass(高斯消元+模板）

复杂度（O(n3)O(n3)O(n^3) 过程：先将增广矩阵转上三角矩阵，然后解x[i];模板：#include<iostream>#include<cstdio>#include<cstring>#include&am

原创 bitset学习笔记

基本操作对于一个叫做foo的bitset： foo.size() 返回大小（位数） foo.count() 返回1的个数 foo.any() 返回是否有1 foo.none() 返回是否没有1 foo.set() 全都变成1 foo.set(p) 将第p + 1位变成1 foo.set(p, x) 将第p + 1位变成x foo.reset() 全都变成0 foo.reset...

原创 ural1002(字典树+dp)

题意：给你一些数字对应一些字母（一对多），然后给一个数字串和一些字母串，问数字串最少可由哪几个字母串组成。 题解：将字母串对应的数字串插入字典树。然后设dp[i]为数字串第i个数字到最后一个数字最少能用几个字母串构成，dp[i]=min(dp[j]),i+1<=j#include<iostream>#include<cstring>#include<c...

原创 数列分块求区间众数（学习笔记）

数列分块1. 区间加法，区间求和简单分析一下：将数列划分为n−−√n\sqrt{n}个块，对于区间修改，整块的进行atag标记一下，非整块的最多只有2n−−√2n2\sqrt n个，暴力加一下 复杂度：O(n−−√)O(n)O(\sqrt n)#include<iostream>#include<cstdio>#include<cmath&gt...

原创 splay操作详解

贴一个好博客传送门

原创 后缀数组学习与应用

学习步骤首先明确几个概念：sa[i]：即，排名为i的后缀的起点下标 rank[i]：即，第i个后缀的排名 height[i]：即，sa[i]和sa[i-1]的最长公共前缀 h[i]：即，height[rank[i]]，即第i个后缀与前一名的最长公共前缀一些结论结论一： 定义：LCP(i,j)=lcp(suffix(sa[i],suffix(sa[j])LCP(i,j)=...

原创 FFT+NTT自学笔记（大数乘积模板）

学习步骤1.回顾一下复数的指数表示形式：常规表示： z=a+ib；z=a+ib；z=a+ib； 指数表示： z=reiθ；z=reiθ；z=re^{iθ}； 其中，r为z的模，θ为辐角主值。 z=r(cosθ+isinθ)=reiθz=r(cosθ+isinθ)=reiθz=r(cosθ+isinθ)=re^{iθ} ，欧拉公式2.快速傅里叶变换要解决的问题：计算多项式乘...

原创 kmp简洁模板（水题入门）

poj2406 题意：求循环节的个数#include<iostream>#include<cstdio>#include<cstring>using namespace std;const int maxn=10000005;char s[maxn];int next[maxn];void getnext(){ int len=...

原创 容斥原理+快速幂求逆元+lucas

原创 最长公共子序列、最长上升子序列、最长公共上升子序列

1.最长公共子序列：核心代码：for(int i=1;i<n;i++) { for(int j=1;j<m;j++) { if(a[i]==b[j]) { dp[i][j]=dp[i-1][j-1]+1; } e...

转载 hdu5296

hdu5296题意：给定一颗边权树，动态地在点集中增加和删除点。求当前所在点集中使得所有点都连通的最小边权和。 解法： 结论：树上点u到以x和y为端点的链的最短距离可以在O(n)预处理和O(logn)时间复杂度内求出。证明：设z=lca(x,y)，dis[u]为树中点u到跟节点的路径权值和。则dis< u,(x,y) >=dis[u]-dis[lca(u,x...