bitwater的空间

//并查集算法 void init(int vn){ for(int i = 0;i <= vn;++i) Father[i] = i; } int find(int x){ return (Father[x] == x) ? x : Father[x] = find(Father[x]); } inline v

经典的NP完全问题,只有暴力解,时间复杂度O(n2^n)对于无向图来说所谓最大团, 其实就是找一个最大完全子图,最大就是包含的点最多.而最大独立集 == 补图的最大团这里使用深度优先搜索实现,对于每一个结点,考虑要与不要两种状态,则问题构成一个子集树,本质上与01背包一样,只不过多了联通性的判断#include const int SIZE = 55;int Gra

/* 图使用向量数组保存,注意初始化edge * 每个点对应的 最大分量的序号 保存在sccNub中 * 可以增加一个vector 存储每个 最大分量中的点 * 所有下标都是 1 ~ n */const int SIZE_N = 10050;vector Edge[SIZE_N];int pre[SIZE_N] , lowlink[SIZE_N] , sccN

参考来源:http://blog.sina.com.cn/s/blog_700906660100v7vb.html这里讨论无向联通图全局最小割方法一  因为两点的最小割其实就是最大流,所以我们可以用dinic算法,枚举源点与汇点, 枚举时间复杂度O(n^2) , Dinic算法 时间复杂度O(n^2m),总的时间复杂度O(N^4M),显然当n很大的时候,需要的时间特

HDU 1402A * B Problem PlusTime Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 18696    Accepted Submission(s): 4246Problem Descrip

参考自http://blog.csdn.net/cclsoft/article/details/5467743循环字符串的最小表示法的问题可以这样描述：对于一个字符串S，求S的循环的同构字符串S’中字典序最小的一个。由于语言能力有限，还是用实际例子来解释比较容易：         设S=bcad，且S’是S的循环同构的串。S’可以是bcad或者cadb,adbc,dbca。

Description满足a^x≡1(mod n)的最小正整数x称为a模n的阶。现给出两个正整数，求x。Input第一行输入k，表示有k组数据之后k行每行两个数a,n(2Output对于每组输入，用一行输出x的值，若不存在输出-1Sample Input22 32 4Sa

#define SIZE 2100struct point_t{ llt x,y;}P[SIZE];//叉积，OA×OBinline llt cross(point_t const&O,point_t const&A,point_t const&B){ llt xoa = A.x - O.x; llt yoa = A.y - O.y; llt xob = B.x - O.

import java.math.*;import java.util.Scanner;import java.util.*;public class Main { public static void main(String[] args) { Scanner cin = new Scanner(System.in); BigInteger []Catalan =

凸包的严格凸多边形Time Limit: 1000ms, Special Time Limit:2500ms, Memory Limit:32768KBTotal submit users: 39, Accepted users: 35Problem 11326 : No special judgementProblem description

http://acm.hrbust.edu.cn/index.php?m=ProblemSet&a=showProblem&problem_id=2315http://acm.hunnu.edu.cn/online/?action=problem&type=show&id=10498以上是验证模板的题目/* * 模板有风险,使用需谨慎 , * 在获取某年的天数时,可以打

PairsDescriptionGiven N integers,count the number of pairs of integers whose sum is less than K. And we have Mqueries. InputThe first behavioris an integer T (1 The next M lines,ea

中国剩余定理又称孙子定理, 主要是为了解线同余方程组x ≡ a1 ( mod  m1 )x ≡ a2 ( modm2 )x ≡ a3 ( modm3 )        .........x ≡ an ( modmn )成立条件是m1,m2, ... ,mn 两两互质, 则对任意的整数：a1,a2, ... ,an，方程组S 有解 并且可以

#include #include #include #include #include #include #include //stl#include #include #include #include #include #include #include #include #include #include #include #include #inc

#include #include using namespace std;int const SIZE = 30100;int const BLOCK_SIZE = 200;//分块大小为接近根号n的整数，这样容易调试struct _t{ int s,e; int idx;};bool operator < (_t const&lhs,_t const&rh

扩展欧几里德算法//非递归的扩展欧几里德算法 //返回a、b的gcd，同时x、y满足ax+by=gcd int_t exEuclid(int_t a,int_t b,int_t&x,int_t&y){ int_t x0 = 1, y0 = 0; int_t x1 = 0, y1 = 1; x = 0; y = 1; int

//基本版#include#include#include#include#includeusing namespace std;const int SIZE = 50;int a[SIZE][SIZE];//增广矩阵int x[SIZE];//解集bool free_x[SIZE];//标记是否是不确定的变元/*void Debug(int equ,int var)

//接口：Gauss();解存在x数组中。equ，var分别表示方程个数和变量的个数，（因为是变量 的个数不算常数的个数）要赋值。a就是行列式，其中常数是在方程右侧的符号。 int const maxn = 2000; const double eps=1e-10; double a[maxn][maxn],x[maxn]; int equ,var; int Gauss()

虽然看不懂原理,但是实现还是不难,留一发模板....//POJ 2417//baby_step giant_step// a^x = b (mod n) n为素数，a,b < n// 求解上式 0 <= x < n的解#include #include #include #define MOD 76543using namespace std;int hs[MOD], hea

#include #include typedef long long llt ;typedef long long int_t;using namespace std;/* * 先写好三个模板, 矩阵快速幂 , 扩展欧几里德 , 中国剩余定理 *///非递归的扩展欧几里德算法//返回a、b的gcd，同时x、y满足ax+by=gcdint_t exEuclid(int_t a

#include using namespace std;typedef int weight_t;const int SIZE_E = //500;const int SIZE_V = //205;int Start,End;struct edge_t{ int node; weight_t c; edge_t *next; edge_t *redg

1.rand()功能:随机数发生器用法:int rand(void)所在头文件: stdlib.hrand()的内部实现是用线性同余法做的，它不是真的随机数，因其周期特别长，故在一定的范围里可看成是随机的。rand()返回一随机数值的范围在0至RAND_MAX 间。RAND_MAX的范围最少是在32767之间(int)。用unsigned int 双字节是655

部分内容参考 http://blog.csdn.net/angelazy/article/details/384892931.点到线段的距离 矢量算法矢量算法过程清晰，如果具有一定的空间几何基础，则是解决此类问题时应优先考虑的方法。当需要计算的数据量很大时，这种方式优势明显。由于矢量具有方向性，故一些方向的判断直接根据其正负号就可以得知，使得其中的一些问

#include #include #include #include #include using namespace std;//S :: p~z//S :: NS;//S :: CSS | DSS | ESS | ISSconst int SIZE = 1111;char s[SIZE];int idx;bool Parse(){ char ch = s[

Area of CirclesTime Limit: 1000ms, Special Time Limit:2500ms, Memory Limit:32768KB Problem description  There are two circles on the plane. Now you must to calculate the area w

普通筛法const int SIZE = 65000;bool isComp[SIZE] = {false};int P[SIZE] = {2};int Pcnt = 1; // 从1开始void sieve(){ for (int i = 2;i < SIZE;++i){ if ( isComp[i] )continue; P[Pcnt+

可以用hdu 1159 测试代码#define N 500/* *下标要从0开始 0~len-1 */int dp[N+1][N+1];const int maxn = 500;int LCSL(char const str1[],char const str2[],int len1,int len2){ int len=max(len1,len2); for(

倍增算法#include #include /* * da函数中的m要小于SIZE； * m可以开很大，前提是SIZE 也要开很大 */int const SIZE = 1000;//分隔符，多串连接时需要用到，第0个为结束符，肯定用到char const DELIMETER[] = {'#'};int const DELIMETER_CNT = 1;//字母表的字母个

代码参考http://blog.csdn.net/zz198808/article/details/7485214//c 的printf 提供了16进制的格式化输出printf("%x",x);#define N 100#define k 16void tran(int num){ int arr[N],i; for (i=0;i <N;i++

#include using namespace std;typedef int weight_t;const int SIZE_E = //500;const int SIZE_V = //205;int Start,End;struct edge_t{ int node; weight_t c; edge_t *next; edge_t *redg

本文部分转自http://blog.csdn.net/acm_cxq/article/details/52038059  二分图又称作二部图，是图论中的一种特殊模型。 设G=(V,E)是一个无向图，如果顶点V可分割为两个互不相交的子集(A,B)，并且图中的每条边（i，j）所关联的两个顶点i和j分别属于这两个不同的顶点集(i in A,j in B)，则称图G为一个二分图。

/*注意初始化init(); *点从0开始计算; 读入数据的时候要-1;输出的时候要加1 * * */int Ecnt = 0;int n,m;int const SIZE_E = //1000010;int const SIZE_V = //10010;struct edge_t{ int v; edge_t *next;}Edge[SIZE_E];

const int SIZE_V //= 600;const int SIZE_E //= 600;typedef int weight_t;int n,m;int start //= 1;int end //= m;//边的结构struct edge_t{ int node; weight_t c;//c为容量 edge_t* next; edg

/*Kruskal 算法 *注意：点是从1到n编号 * */typedef int weight_t;const int SIZE_V = 101;const int SIZE_E = SIZE_V / 2 * SIZE_V;int Father[SIZE_V];int Ecnt = 0;//int n; // 此处写点的全局变量//并查集算法void init(int

参考自 http://www.cnblogs.com/ikids/articles/4660011.htmlhttp://blog.csdn.net/ggggiqnypgjg/article/details/6645824/算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串

#include #define MEM(a,x) memset(a , x , sizeof(a))using namespace std;typedef int weight_t;const int SIZE_E = //1000000;const int SIZE_V = //1050;int n;int ECnt = 0;struct edge_t{ int

#include #include #define M 310#define inf 0x3f3f3f3fint n,nx,ny;int lx[M],ly[M]; //lx,ly为顶标，nx,ny分别为x点集y点集的个数int _visx[M],_visy[M],w[M][M],_slack[M],_link[M];int DFS(int x){ _visx[x]

利用三点坐标可以直接求面积向量法向量AB = [x2 - x1, y2 - y1]向量AC = [x3 - x1, y3 - y1]三角形ABC的面积 = 0.5|AB||AC|sin(角BAC)【|AB|为向量AB的模，也就是向量AB的长度】= 0.5*|向量AB与向量AC的叉积|= 0.5*|(x2 - x1)(y3 - y1) - (y2 - y1)(x3 - x

/* * 多边形的交，多边形的边一定是要按逆时针方向给出 * 还要判断是凸包还是凹包，调用相应的函数 * 面积并，只要和面积减去交即可 */#include using namespace std;const int maxn = //300;const double eps = 1e-8;int dcmp(double x){ if(x > eps) return

#include /*区间从0 ～ n标号 *注意Max的二维长度 */int A[100050];int Max[100050][20];void ST(int n){ int i,j; //初始化，长度为1的区间 for(i=0;i<n;i++) Max[i][0] = *(A+i); //类似倍增法 for(j=1;(1<<j