du_lun的博客

#define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include using namespace std; int N; // 1<= N <= 20000const int maxn = 20000;vector tree[maxn + 5]; // tree[i

参考于该博客 int Mark[MAXSIZE]; int prime[MAXSIZE]; //判断是否是一个素数 Mark 标记数组 index 素数个数 int Prime(){ int index = 0; memset(Mark,0,sizeof(Mark)); for(int i = 2; i &lt; MAXSIZE; ...

传送门#include&lt;bits/stdc++.h&gt;using namespace std;#define fir first#define sec secondconst int N=2e5+10;typedef pair&lt;int, int&gt; pii;pii e[N];set&lt;int&gt; G[N];bool removed[N];int ...

传送门题意：给一个只有‘1’~‘9’的串，求不同回文串代表数值的和思路：马拉车找回文串，hash标记不同，hash方法求数值#include&lt;bits/stdc++.h&gt;using namespace std;const int N=2e6+10;const int mo=1e9+7;typedef long long ll;typedef unsigned l...

传送门思路是参照着这个博客写的https://blog.csdn.net/a1035719430/article/details/80488083dp的一个状态是点集，这个是很好想的，但是加上路径上宝藏的数量就不好想了，因为所有路径都挖完是一颗树，我们可以用一维表示状态中最大的深度，即dp[i][s]表示深度为i，包含s集合点的时候最小值转移就是dp[i][t]=min(dp[i][t...

https://www.luogu.org/problemnew/show/P3773这个题就用到了上篇博客https://blog.csdn.net/du_lun/article/details/82414086，遍历子集的小技巧，题意是找有多少个长度&gt;=2的下降子序列，满足mod2的lucas，并且值大于0由lucas可以知道C(x,y)=C(a1,b1)*C(a2,b2)…...

#include&lt;bits/stdc++.h&gt;using namespace std;const int N=1e5+10;struct Tire{ int nxt[N][2], cnt[N], ed[N]; int L, root; void init(){ L=root=0; nxt[root][0]=nxt[root...

题意：有n个人，每个人猜一个[1,6]长度为L的序列，然后无限制的掷骰子，先掷到的那个序列，那个人就是winner，问每个人赢得概率。思路：骰子是掷无限次的，但是有许多状态是重复的，每一种掷法最终必定是以某一个人的序列结束的，也就是说，谁都不赢的序列是无穷的，概率趋近于0，所以每个人概率相加是1.我们可以用ac自动机节点编号当做某种情况的状态，nxt数组就是这个状态转向的下个状态。...

传送门思路：首先需要了解Matrix_Tree定理1.设 G是无向图的邻接表，D是无向图各个点的度数2.令M=D-G，则M的任意n-1阶余子式的行列式值的 绝对值， 就是无向图生成树的个数应注意的是邻接表的值不只是0,1值，也就是说，当图出现重边的话，邻接表该位置累加。这道题思路，对于最小生成树，克鲁斯卡尔算法是将边从小到大排序，对于同一大小的边，在最小生成树中用到的数量是...

 题意：有一个坐标轴，从s点开始，为了到达e点，走的时候是周而复始循环的，比如 0 1 2 3 3 2 0.。。，从一个点开始，每走k步都有pk的概率停下，问正好停在目标点的期望是多少；思路：设dp[i]是在i点到e点的期望，dp[i]=p1*dp[i+1]+p2*dp[i+2]....+pk*dp[i+k](因为期望都是倒着推)gauss 求解就可以了因为有不能到的点，如果不能到的点...

title: 二维线段树模板date: 2018-10-31 15:21:44tags:- 二维线段树categories: “算法”这个模板是根据HDU-4819编写的。功能是用二维线段树进行二维点修改和区间查询最小者和最大值。二维线段树的思想跟一维一样，只不过一维是二分二叉树，二维是四分四叉树。具体实现看代码。#include&lt;bits/stdc++.h&gt;us...

转载博主链接前言：    在ACM比赛中，大家可能会遇到一些读入相当大的题目，此时前人告诉我们：遇到这些题目，我们最好不要用cin/cout，因为cin/cout在读入输出的过程中相对与scanf/printf来说比较慢。    但是还是有那么些个善（丧）良（心）的出题人还是喜欢卡读入的常数，可能我们的scanf/printf，甚至普通的read（）函数也无法支撑。此时要如何办呢？于是...

传送门huffman模板，如果每次取k个能取完，num=k，取不完先将n%（k-1）取完，再每次取k#include&lt;bits/stdc++.h&gt;using namespace std;const int N=1e5+10;typedef long long ll;int arr[N], n, limit;/****************************/...

#include&lt;iostream&gt;#include&lt;cstdio&gt;#include&lt;cstring&gt;using namespace std;int nxt[110];char s[110], s1[100];void getnxt(){ int j=nxt[0]=-1, i=0; int len=strlen(s); whi...

//以i为终点 嵌套矩形问题 int dfs(i){ if(d[i]&gt;0) return d[i]; d[i]=1; for(int j=1; j&lt;=n; j++) { if(G[i][j]) d[i] =max(d[i], dfs(j)+1); } return d[i];}

#include&lt;cstdio&gt;using namespace std;typedef long long ll;ll extend_gcd(ll a, ll b, ll &amp;x, ll &amp;y){ if(!a&amp;&amp;!b) return -1; if(b==0){ x=1; y=0; return a; } ll d=extend_gcd(...

传送门题意是求，n是1e15，枚举1的个数，数位dp#include&lt;stdio.h&gt;#include&lt;iostream&gt;#include&lt;string.h&gt;#include&lt;algorithm&gt;using namespace std;typedef long long ll;const int N=1e5+10;const ...

传送门偏序关系二维树状数组这个题用dp[i][j]表示前i个数，以h[i]结尾，用了j次修改，最长的非递减序列。但是这样直接dp的话如下for(int i=1; i&lt;=n; i++){ for(int j=0; j&lt;=K; j++){ for(int a=0; a&lt;i; a++){ fo...

陌上花开，可缓缓归矣偏序关系就是按照某个条件的排序关系，这个题是三维偏序关系，然后给这些花划分等级，也就是说,A(x, y, z) 比A'(x',y',z')等级高的话 那么就满足x&gt;=x', y&gt;=y', z&gt;=z'。我们只要统计一下比自己等级低的有几个，按照这个数量做一个等级归纳排序，输出每个等级的个数就解决了。解决这类问题可以用cdq+树状数组，首先我们对x排序，然...

传送门思路：用KMP和AC思路是一样的。ac的思路就是跑文本串，记录下每个节点的最小值，这样会省时间，因为fail指针的存在，我们不需要管模式串，kmp的话就是每个模式串都要匹配，取最小值；#include&lt;bits/stdc++.h&gt;using namespace std;const int N=1e5+10;int val[N];int nxt[500010][2...

传送门挑选m个数 ，和为k，数的取值为0~n-1，相当于将每个数加1，使挑的m个数和是k+m，可以看成k+m个1，用m-1个隔板，隔成m分，所以总数为C(m+k-1, m-1)，减去不合法的，就是大于n的情况。先看只有一个位置大于n，有C(m,1)*C(m+k-n-1, m-1)种情况，但是有重复，比如1，3，2，1   。  3，1，2 ，1这两种情况，n是2，k是8的话，3、3、2、1是一...

 传送门 #include&lt;bits/stdc++.h&gt;using namespace std;typedef long long ll;const int N=1e5+10;const int mo=1e9+7;ll n, P, F[N], C, D, inv[N];struct Matrix{ ll m[3][3];};void init(){...

HDU1695莫比乌斯反演两个公式思路：该题的题意是给你两个范围，1~n，1~m，求x属于1~n，y属于1~m， 且GCD（x，y）==k，这样的有多少对。等同于求1~n/k, 1~m/k, gcd(x,y)==1,的对数。令f（t）为gcd（x， y）等于t的个数，F（n）等于gcd（x， y）=t的倍数的个数。所以可以用第二个公式，就是说gcd（x,y）等于t的倍数的所...

传送门这道题真的是坑，竟然用ull，ll都会爆。首先这是一个仙人掌图。题意让求任意两点最大流，再进行异或。先说最大流，对于环上任意两点最大流，就是两条路径最小的边和，这就等效于求出环中最小的边，使环上任意一条边都加上这条边，然后去掉最小边。这样会求出一棵树。我们将边按从大到小排序，每次向空图里加一条边（用并查集维护），那么加的这条边就是两个连通分量的最大流。还是老方法。对于仙人掌图，（...

该题有很多种做法，乱搞可以过，最好的方法应该是扫描线，本题是扫描线模板题官网数据第一个错了，所以83分就是满分#include<bits/stdc++.h>using namespace std;const int N=1e4+10;#define mp(x,y) make_pair(x,y) #define fir first #define sec sec...