- 博客(73)
- 问答 (2)
- 收藏
- 关注
RSS订阅
原创 计算机图形学几何工具算法详解（第2章）：(hua)笔(shui)记
还是去那里看好了第2章 矩阵和线性系统2.3.2 转置性质：(AB)^T=B^T A^T2.3.4 矩阵乘法性质：没有交换律：A(BC)!=(AB)C奇妙的转置：2.4.1 线性方程2.4.3 一般线性系统增广矩阵。2.4.4 减行、消元、轶等等减行、基本行运算、线性无关2.5.1 对角矩阵
2017-10-01 18:52:26
1190
1
原创 每日亿题之20220904
正好sz*(sz-1)/2=0+1+2+…sz-1，所以可以对k从大到小的因数一次赋予权值0,1,2,…,sz-1，放在主席树上统计即可得到lcm==k的三元组个数之和。给定n,k=L的因数个数sz_k，则sz_k。(sz_k-1)/2为(i,j)pair的个数。因为lcm为k的倍数，所以仅lcm。
2022-09-05 00:04:15
226
原创 每日十(?)题之20220903
经过超过一年半的颓废生活，某csdn用户不可避免地从一个拉跨的OIer编程了彻底的fw。已经到来的大三生活让他感到一些压力(…)，暂且为了今年的icpc不要那么划水🤨，某人决定重拾OI刷题，开始这个每日x题计划。今天的每日10题，是从cf随便选的10题。由于博主已经太久没有碰过OI，于是没有自信地先选了一道cf1700分题，以及一个1900分的。第1题。
2022-09-04 04:00:00
219
原创 记录4 最近贺题时做的题 cf611e,cf873f,cf264d,cf1320d,cf1438c,hdu6833,hdu3602,cf1187g,cf1051e,cf925d,cf840c,etc
cf611e简单题，三个火枪手分别有a，b，c 的攻击，打怪，每个怪有生命值，每回合三个火枪手各选一个怪打，怪只有一回合里被各个火枪手打的累计伤害大于等于生命值才会死。求最少多少回合打怪#include<bits/stdc++.h>#define FOR(i,s,t) for(int i=s;i<=t;++i)#define REP(i,t,s) for(int i=t;i>=s;--i)#define RESET(a) memset(a,0,sizeof a)using
2020-12-16 23:58:29
342
原创 最近做题遇到的坑·三 cf429c,cf377d,cf288d,cf294e,cf958b2,cf765e,cf812e,gym101981d,gym101981g
cf429c给一个序列a1..ana_1..a_na1..an,能否找到一棵n个节点的，且每个非叶子节点有至少2个叶子的树，使每个节点子树大小一一对应于序列a中的一个数n<=24暴力都不会写。。直接搞搞就好了#include<bits/stdc++.h>using namespace std;int n,a[25],tot;int st[12],top;void dfs(int d) { if(d>tot) cout<<"YES",exit(0)
2020-12-01 23:45:36
336
原创 最近做题遇到的坑·第二弹 cf976d,agc022_c,cf1000f,cf1039d,cf1042f,cf850c,cf727f,cf761f,cf501e,cf735e
cf976d一道yzy两年前就秒了的题我今天还不会而且竟然看不懂yzy写的题解构造题，构造方法不少，yzy的看上去最nb，不过我哪一种都不会，写了个random_shuffle竟然过了agc022_c没什么营养的 思维题，没有什么难度cf1000f5×1055 \times 10^55×105级别，每次找一个区间只出现了一次的数重新学习主席树原来主席树每颗树里的叶子存的是pre[]，前一次出现的位置原来主席树modify(int &u,int v \quad u和v可以是 一样
2020-11-13 23:36:57
256
原创 我用treap写了5k的“大数据结构题”被动态开点线段树爆艹
https://codeforces.com/contest/803/submission/97759561过程中产生了若干bug先是mn0没置为inf（小问题）treap中的split不管运算符重载为什么，都要是***<***，不能出现大于对一个块两头都 真包含 所求块 的询问没处理大的线段树如果已经覆盖了，就不能用小的平衡树里查询的结果数组开小了整体上bug还是比想象中少，想必是treap比较好写（想必是因为treap是贺板子的）cf竟然没有submission按memery排的
2020-11-06 16:32:42
140
原创 最近做题遇到的坑 cf865d,cf838e,cf317d,cf1394b,poj3322,cf223c,hdu2435,cf47d
https://vjudge.net/problem/CodeForces-865Dn<=3e5n<=3e5n<=3e5,a[i]<=1e6a[i]<=1e6a[i]<=1e6,每天可以花a[i]a[i]a[i]买一个单位，或卖一个单位赚a[i]a[i]a[i]，最后不留下单位，求最大赚的值（最后不留下单位 这条件是没用的）（每天计为卖了以后还可以买一次，相当于没有卖）（不买不卖相当于当天买当天卖）（可看作每天都买或卖）（买了还可以卖1次（退回），又一次（
2020-10-30 19:26:36
325
原创 UNKNOWN ERROR 诡异错误
#includeusing namespace std;struct Long { int *s; Long() { s=new int[5]; memset(s,0,20); } ~Long() { delete[] s; } void operator=(const Long &A) { delete[] s; s=new i
2018-05-08 08:52:05
7260
1
转载 HDU 4658 Integer Partition（生成函数）
https://blog.csdn.net/V5ZSQ/article/details/79943615
2018-04-16 11:30:17
237
原创 bzoj3028 食物
https://blog.csdn.net/popoqqq/article/details/42805975#includeint main() { int x=0; char c; while ((c=getchar())>'-') x=(x*10+c-'0')%10007; printf("%d\n",x*(x+1)%10007*(x+2)%10007*1668%10007)
2018-04-15 19:27:59
217
原创 洛谷 P3224 [HNOI2012]永无乡
#include#define N 100005using namespace std;int n,m,q;int a[N],pre[N];int rt[N],lc[N],rc[N],sz[N];int find(int x){ return x==pre[x] ? x : find(pre[x]);}void up(int u){ sz[u]=sz[lc[u]]+sz[r
2018-02-28 20:17:22
292
1
原创 非旋treap
#include#define N 100005using namespace std;int rt,val[N],lch[N],rch[N],pro[N],sz[N],cnt;int rnd(){ static int sd=4234; return sd*=3423;}int node(int v){ val[++cnt]=v; pro[cnt]=rnd(); sz[
2018-02-25 10:27:36
229
原创 无向连通图的割点和桥
割点#include#define N 100005using namespace std;int n,cnt;int hd[N],nxt[N<<1],to[N<<1],tot;int dfn[N],low[N],now;int ans[N];void add(int a,int b){ to[++tot]=b; nxt[tot]=hd[a]; hd[a]=tot;
2018-02-25 09:03:57
243
原创 分块*2
洛谷 P3203 [HNOI2010]BOUNCE 弹飞绵羊#include#include#define N 200255#define M 450int n,q,a[N];int stp[N],nxt[N];void build() { memset(nxt,66,sizeof nxt); scanf("%d",&n); for(int i=
2018-01-06 13:08:53
229
原创 洛谷 P3806 【模板】点分治1
#include#include#include#define K 10000005#define N 10005using std::max;int n,m;int hd[N],nxt[N<<1],to[N<<1],w[N<<1];int vis[N];int rt,sz[N],siz[N];int son[N],dfn[N],cnt;bool mp[K];void ad
2018-01-06 13:06:48
554
1
原创 树链剖分：洛谷 P2590 [ZJOI2008]树的统计
一共三部分：树剖，线段树，LCA难打主要是因为代码长#include#include#include#define mid (l+r>>1)#define N 30005#define M 19#define INF 1128481604int n,q;int mx[N<<2],sum[N<<2];int head[N],next[N<<1],to[N<<1];int
2017-12-30 17:59:50
238
原创 树上区间和1
#include#define lowbit(x) (x&-x)#define N 100005int n,q,w[N];int head[N],next[N<<1],to[N<<1];int fa[N],son[N],dep[N],size[N],top[N],num[N],pos[N];void add(int a,int b) { static int cnt=0; to[+
2017-12-28 18:58:39
246
原创 洛谷P2633 Count on a tree(RE)
#include#include#define N 100005using namespace std;int n,m,x,y,z,last,w[N],s[N],len;int head[N],next[N<<1],to[N<<1];int pre[N][20],dep[N];int rt[N],ch[N*20][2],size[N*20],tot;void add(int a,i
2017-12-16 14:29:53
382
原创 洛谷P2617 Dynamic Ranking 树状数组套主席树
#include#include#include#define lowbit(x) (x&-x)#define N 10005using namespace std;int n,m,a[N],b[N*2],s[N][3],len;char c;int rt[N],val[N*400],lch[N*400],rch[N*400],cnt;int L[20],R[20],llen,r
2017-12-15 19:17:35
285
原创 郁闷的出纳员
splay真的很慢#include#include#define N 100005#define INF 100000007int m,lowst,o,tmp,leave;char c;int root,fa[N],ch[N][2],val[N],size[N],cnt;int son(int t) { return t==ch[fa[t]][1];}void pushu
2017-12-08 19:07:04
270
原创 文艺平衡树
#include#define N 100005int n,m;int root,fa[N],ch[N][2],size[N],rev[N];void swap(int &a,int &b) { a^=b, b^=a, a^=b;}int son(int t) { return t==ch[fa[t]][1];}void pushup(int t) { si
2017-12-07 18:22:19
241
原创 树状数组 区间修改 区间查询
#include#include#define lowbit(x) (x&-x)#define LL long long#define N 100005int n,q,x,y,z;LL a0[N],a1[N],a2[N];char c;void read(int &t) { t=0; int f=1; for(;!isdigit(c=getchar());) if(
2017-12-01 22:48:59
196
原创 untitled(1)
#include#define LL long long#define N 233333#define mod 11111111inline int pow(int a,LL x) { int ret=1; while(x) { if(x&1) ret=(1LL*ret*a)%mod; x>>=1; a=(1LL*a*a)%mod; } return ret;}inli
2017-12-01 18:10:48
310
原创 神奇的东西
有一集合s:s={x∈N|x=(2^n1)*(3^n2),n1∈N,n2∈N}有以下程序：k=1, tmp=0(REPEAT FOR n TIMES){EITHER k←k+tmpOR tmp←k}则：MAX k = nTH ELEMENT IN s（所有程序中k的最大值为s中的第n个元素）
2017-11-30 23:08:33
194
转载 假的2048
#include<stdio.h>#include<Windows.h>#include<time.h>#include<conio.h>#include<algorithm>using namespace std;namespace __2048{ const int d[][2]={{0,1},{0,-1},{1,0},{-1,0}}; char map[4][4]
2017-11-11 11:52:45
171
原创 我的第二棵平衡树: 普通Treap
这又是洛谷平衡树的模板题...[SBT的写法]Treap的代码很短,只是在普通二茬查找树的插入操作和删除操作中加入维护堆性质的rotate,形式也和堆很像,插入时将新节点不断向上抬,删除时将节点向下降直至弹出.随机种子的选取可以用一些玄学的方法(来自网上某位大佬)inline int random() { static int seed=48271; return
2017-11-09 13:30:22
254
原创 我的第一棵平衡树-SB树
数组实现的SB树(而且竟然没有结构体)功能就和普通的二茬查找树一样:insert: 插入元素erase: 删除元素pre: 查询一个元素的前驱suc: 查询一个元素的后继select: 返回第几小的元素的值rank: 返回一个值的位置除了maintain全部非递归实现,并且有assert(似乎并没有什么用)maintain写得不是很高效,并且缺储存平衡树中元
2017-11-07 16:32:45
2245
1
原创 草稿
#include#include#define N 100005int root,val[N],lch[N],rch[N],w[N],size;void rrot(int &t) {int x=lch[t]; lch[t]=rch[x]; rch[x]=t;w[x]=w[t]; w[t]=w[lch[t]]+w[rch[t]]+1;t=x;}void lrot(
2017-11-06 21:00:16
242
原创 [洛谷]最大流 模板
P3376 【模板】网络最大流开始我打的是假的网络流,深搜宽搜的临时变量存的都是边明显变比点要多得多这是EK#include#include#define INF 1000000009#define N 10010int n,m,s,t,x,y,z,pre[N*20],way[N*20];int head[N],next[N*20]
2017-11-02 18:54:00
358
转载 P3382 【模板】三分法
#include#define eps 1e-7int n;double l,r,a[14];double getval(double x) { double ret=0; for(int i=n;i>=0;i--) { ret*=x; ret+=a[i]; } return ret;}double solve() { while
2017-11-02 18:44:05
329
原创 RPFS(RP First Search) 模拟退火算法 未完
POJ 2420求费马点代码来自网络#include#includeconst double inf=1e9;const int N=100;const double threshold=1e-7,delta=0.98,inital_temperature=100;const double d[4][2]={{0,1},{1,0},{0,-1},{-1,0}};int n;
2017-10-30 19:48:05
1451
1
原创 LCA的三种求法 - HDU 2586
Tarjan一遍 大法师 处理所有询问, 每次递归处理子树, 把子树与当前节点的并查集合并. 然后如果一个询问的两个节点都被访问, 直接输出答案#include#include#includeusing namespace std;int T,n,m,a,b,c,p[40010],vis[40010],dep[40010],ans[210];vector >g[40010];v
2017-10-30 17:37:29
254
原创 模板:强连通分量和缩点
洛谷 P2002 消息扩散缩点后输出入度为0的点(tarjan)#include#include#include#define N 100010#define M 500010using namespace std;int n,m,x,y,Q[M],dfn[N],low[N],inq[N],p[N],in[N],tail,now,size,ans;vectorg[N];
2017-10-30 14:46:13
332
原创 扩欧习题
•https://vjudge.net/problem/LightOJ-1306•http://acm.hdu.edu.cn/showproblem.php?pid=1576•http://acm.hdu.edu.cn/showproblem.php?pid=5651•http://www.lightoj.com/volume_showproblem.php?problem=130
2017-10-29 21:16:08
400
空空如也
Ubuntu虚拟机Firefox“server not found"
2020-09-20
Ubuntu虚拟机Firefox“server not found"
2020-09-20
TA创建的收藏夹 TA关注的收藏夹
TA关注的人