DongResponse-CSDN博客

原创半平面交

#include <cstdio>#include <cmath>#include <algorithm>using namespace std;const double eps = 1e-8;int Sign(double x) { if (fabs(x) < eps)return 0; return x < 0 ? -1 : 1;}struct Point { double x, y;} ans[505];.

2020-09-07 21:10:04 138

原创树状数组求最值+单点修改

#include <cstdio>#include <algorithm>using namespace std;const int max_n = 200005;int n, a[max_n], tree[max_n];inline int low_bit(int x) { return x & -x; }inline void Change(int p, int x) { a[p] = x; for (int i = p; i &l.

2020-09-01 21:27:19 219

原创可持久化fhq

#include <cstdio>#include <cstdlib>#include <algorithm>using namespace std;const int max_n = 5e5 + 5;struct Fhq { int ch[2], val, key, sze;} fhq[max_n * 50];int tot = 0, root[max_n];inline int New_Code(int val) { int.

2020-08-25 14:35:49 269

原创极角排序+凸包

#include <cstdio>#include <algorithm>#include <cmath>using namespace std;int n;const int MAXN = 1e5 + 5;struct Point { double x, y;} p[MAXN], Stack[MAXN];int Top = 0;bool cmp1(Point u, Point v) { if (u.y == v.y)retu.

2020-08-17 13:05:54 182

原创海伦秦九韶

#include <cstdio>#include <cmath>using namespace std;const double eps = 1e-8;struct Point { double x, y;};struct Vector { double x, y;};double Distance(Point A, Point B) { double dx = A.x - B.x; double dy = A.y - B.

2020-08-15 19:45:25 256

原创模拟退火

#include <cstdio>#include <cmath>#include <algorithm>using namespace std;const double delta = 0.996;const double eps = 1e-15;double X, Y;struct Point { double x, y;} a[1005], ans;int n;double ans_val = 0;double Min_Va.

2020-08-14 20:54:28 103

原创两圆相交面积

#include <cstdio>#include <cmath>using namespace std;const double PI = acos(-1);struct Point { double x, y;};struct Round { Point o; double r;};double Cos(double a, double b, double c) { return (a * a + b * b - c *.

2020-08-14 16:50:14 143

原创逆矩阵

#include <cstdio>#include <algorithm>using namespace std;#define int long longconst int mod = 1e9 + 7;int n, a[405][405], b[405][405];int Power(int q, int c) { if (c == 0)return 1; if (q == 0)return 0; int ans = Power(q .

2020-08-12 20:35:55 93

原创点分治

#include <cstdio>#include <algorithm>using namespace std;const int max_n = 1e4 + 5;int n, k;struct Edge { int to, next, val;} edge[max_n * 2];int head[max_n];int ans = 0;bool vis[max_n];int mx_son[max_n], sze[max_n];int F.

2020-08-07 19:19:23 75

原创三点确定一个圆

struct Point { double x, y;} a[max_n];struct Round { double r; Point o; bool On_Round(Point p) { return (p.x - o.x) * (p.x - o.x) + (p.y - o.y) * (p.y - o.y) <= r * r; }};Round Make_Round(Point p) { return {0, p};.

2020-08-03 20:09:54 394

原创后缀数组

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;const int max_n = 1e5 + 5;int n, m;char s[max_n];int des[max_n], d = 0;int Rank[max_n], Sa[max_n], Height[max_n], tax[max_n], tp[max_n];void Rank_S.

2020-07-17 20:56:56 89

原创广义后缀自动机

#include <cstdio>#include <cstring>#include <algorithm>#include <queue>using namespace std;const int max_n = 1e6 + 5;char s[max_n];struct Trie { int ch[26]{}, fth, pos; Trie() { for (int &i:ch)i = 0;.

2020-07-09 12:01:30 148

原创 Manacher

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;const int max_n = 1e5 + 5;char s[max_n], a[max_n * 2];int n = 0, t, len[max_n * 2], pre[max_n], suf[max_n];signed main() {// freopen("in", "r", std.

2020-06-02 20:33:00 106

原创回文自动机

#include <cstdio>#include <cstring>using namespace std;const int max_n = 500005;struct Palindromic_Tree { int next[max_n][26]{}, fail[max_n]{}; //next指针，next指针和字典树类似，指向的串为当前串两端加上同一个字符构成 //fail指针，失配后跳转到fail指针指向的节点 int cnt.

2020-05-19 21:38:00 134

原创大整型

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define int long longconst int base = 1e18;struct BigInt { int num[1000]{}, len; BigInt() { len = 0; } void Assign(int x) { len = .

2020-05-15 15:27:54 111

原创高斯消元

#include <cstdio>using namespace std;int n;double a[105][105], b[105];signed main() { scanf("%d", &n);//NOLINT for (int i = 1; i <= n; i++) { for (int j = 1; j <= n; j++)scanf("%lf", &a[i][j]);//NOLINT sc.

2020-05-15 11:28:13 157

原创高次同余方程

int Power(int q, int n, int p) { if (n == 0)return 1 % p; if (q == 0)return 0; int ans = Power(q * q % p, n / 2, p); if (n & 1)ans *= q;//NOLINT return ans % p;}int BSGS(in...

2020-04-26 18:38:22 240

原创最小斯坦纳树

#include <cstdio>#include <algorithm>#include <queue>using namespace std;const int INF = 0x3f3f3f3f;int n, m, k, start[105], a[15], dp[105][1 << 10];//NOLINT;struct E...

2020-04-23 21:29:32 151

原创线段树合并

#include <cstdio>#include <vector>#include <algorithm>#include <queue>using namespace std;const int MAXN = 1e5 + 5;int n, m;struct Edge { int to, next;} edge[MAX...

2020-04-21 20:49:46 79

原创后缀自动机

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;const int MAXN = 1e6 + 5;char s[MAXN];int n = 0;int tot = 1, last = 1;struct Suffix_Automaton {...

2020-04-17 11:31:28 107

原创 Splay

int n, root = 0;struct Tree { int lef, rgh, num, sze, fth;} tree[MAXN];inline void Rotate(int x) { int y = tree[x].fth; int z = tree[y].fth; if (root == y)root = x; else if (t...

2020-04-06 10:54:39 83

原创杜教筛

#include <cstdio>#include <unordered_map>using namespace std;const int N = 5000000;int prime[N + 5], min_prime[N + 5], cnt = 0;long long pre_miu[N + 5], pre_phi[N + 5];unordered_ma...

2020-04-05 20:55:26 106

原创 Mobius 2

求 1 <= i <= n, 1 <= j <= m 范围的 d(i * j) 的和，d(x) 为 x 的因子个数#include <cstdio>#include <algorithm>using namespace std;#define int long longconst int MAXN = 50005;int prim...

2020-04-05 16:44:52 87

原创 Mobius

求 1 <= x <= n, 1 <= y <= m 中 gcd(x, y) 为质数的个数#include <cstdio>#include <algorithm>using namespace std;const int MAXN = 1e7 + 5;int is_prime[MAXN], miu[MAXN], prime[MAXN...

2020-04-05 15:25:53 95

原创拓展Lucas

#include <cstdio>using namespace std;#define int long longint P[50], D[50], cnt = 0, B[50];int Power(int q, int n, int mod) { if (n == 0)return 1; if (q == 0)return 0; int res...

2020-04-04 13:54:45 118

原创 Lucas

int Power(int q, int n, int p) { if (n == 0)return 1 % p; return Power(q * q % p, n / 2, p) % p * (n & 1 ? q : 1) % p;}int Factorial(int n, int p) { if (n >= p)return 0; if ...

2020-04-03 19:02:44 117

原创线性同余方程

int Extend_Gcd(int a, int b, int &x, int &y) { if (b == 0) { x = 1, y = 0; return a; } int d = Extend_Gcd(b, a % b, x, y); int z = x; x = y; y = z - y ...

2020-04-03 18:27:29 149

原创 fhq

struct Treap { int lef, rgh, pos;} fhq[100005 * 205];int root = 0, tot = 0;queue<int> que;int tmp[100005], cnt = 0;int Build(int l, int r) { if (l > r)return 0; int ans = q...

2020-03-26 14:01:52 429

原创左偏树

#include <cstdio>#include <algorithm>using namespace std;#define int long longconst int MAXN = 100005;int n, m;int fth[MAXN], Val[MAXN], Pow[MAXN];struct Edge { int to, next;...

2020-03-25 21:17:07 67

原创 Treap模板

#include <iostream>#include <algorithm>const int INF = 0x7fffffff;const int MAXN = 100005;using namespace std;int root = 0;struct { int lef, rgh; int val, key; int cnt...

2020-01-12 15:48:02 81

原创等比数列求和模板

#define mod ....//随便填int qpow(int n, int m) { int ans = 1; while (m) { if (m & 1)ans *= n, ans %= mod; m /= 2; n *= n; n %= mod; } return ans;}int ers(int q, int n) { if (n == 0)retur...

2019-09-23 21:00:10 153

dong_pope的博客