2014年04月_utoppia

原创模板，凸包

int TuBao(P p[],int n,P q[]){ int m = 0; sort(p,p+n); for(int i=0;i<n;++i) { while (m>1 && dcmp(cross(q[m-2],q[m-1],p[i]))<=0) m--; q[m++] = p[i]; } int k = m;

2014-04-30 10:24:45 631

原创 HDU 1817 , POJ 1286

polya // hdu 1817 poj 1286 // by utoppia // polya #include typedef long long ll; ll pow(ll a,int n) { ll ret = 1ll; while(n--) ret *= a; return ret; } int gcd(int a,int b) { retur

2014-04-22 16:16:57 874

原创模板，扩展，离散对数，数论

扩展版的离散对数基本思路还是基于 Shank的大小步算法(Shank's Baby-St

2014-04-22 12:01:13 728

原创模板，FFT 快速傅里叶变化

大神的模板。。

2014-04-21 20:17:32 4566

原创 CSU 1416 Practical Number 很有意思的数学题

基本上就是一个结论题，然后wiki一下就出了 // csu 1416 // by utoppia // http://planetmath.org/PracticalNumber // http://en.wikipedia.org/wiki/Practical_number // http://oeis.org/A005153 #include #include #include #inc

2014-04-21 12:12:24 967

原创一些有意思的网站

找规律的 oeis.org

2014-04-21 10:55:11 627

原创模板，无向图强连通缩点，tarjan

int first[nMax],nxt[nMax],to[nMax],e; void addadge(int u,int v) { to[e]=v;nxt[e]=first[u];first[u]=e;e++; to[e]=u;nxt[e]=first[v];first[v]=e;e++; } int low[nMax],dfn[nMax],dep,belong[nMax],vis

2014-04-17 16:48:45 2559 1

原创模板，pollard_rho大数因数分解

大数分解，没有模板你就OUT啦。。。

2014-04-16 21:31:27 1577

原创模板，树状数组

树状数组的神奇在于省空间和时间；空间复杂度O（

2014-04-16 21:12:18 533

原创模板，线段树，矩形轮廓周长

原理和面积并是一样的，注意有了两边的长要计算；

2014-04-16 11:17:01 569

原创模板，线段树，面积并

具体的思路有空补上，先贴代码 #include #include #include using namespace std; int const nMax = 100000; struct T{ int y[2]; double f[2]; double x; int left; void make(double a,double b,double c,

2014-04-16 10:41:59 549

utoppia is another word .