- 博客(156)
- 收藏
- 关注
RSS订阅原创 【Go自学版】02-goroutine
利用时间片分割进程，致使宏观上A,B,C同时执行（并发）CPU利用率包含了，进程/线程的数量越多，切换成本也会增大最大并行数：GOMAXPROCSwork stealing: 偷其他队列的Ghand off: 当前G1阻塞，创建/唤醒一个thread，移动当前本地队列至其上。
2023-12-08 16:17:06
226
原创 2021CCPC 桂林 F. Illuminations II （计算几何）
F. Illuminations IISolution内多边形每条线与外多边形做交点，利用前缀和计算即可。还是挺容易想到的，比赛的时候没时间写了，写计算几何再也不用cin了。Code/* /へ ／| /＼7 ∠＿/ / │ ／ ／ │ | ＿,＜ ／ /`ヽ │ ヽ / 〉 \ ` / /
2021-11-15 15:50:21
465
原创 2018 China Collegiate Programming Contest - Guilin Site L- Two Ants (计算几何)
L. Two AntsSOLUTION注意讨论CODE/* /へ ／| /＼7 ∠＿/ / │ ／ ／ │ | ＿,＜ ／ /`ヽ │ ヽ / 〉 \ ` / / ｲ● ､ ● ⊂⊃〈 / () へ | ＼〈 >
2021-10-30 14:12:12
345
原创 P1337 [JSOI2004]平衡点 / 吊打XXX 模拟退火
平衡点 / 吊打XXXsolution记录一下模拟退火的几个题目code/*SiberianSquirrel*//*CuteKiloFish*/#include <bits/stdc++.h>//#include<bits/extc++.h>#include<ext/rope>#include<ext/pb_ds/assoc_container.hpp>#include<ext/pb_ds/tree_policy.hpp>usi
2021-10-21 14:49:42
268
原创 SPOJ 10606 数位DP
SPOJ-10606solution坑点在于前导零对于统计 000 的个数的影响.dp[pos][num][appear]dp[pos][num][appear]dp[pos][num][appear]num−>100000101num->100000101num−>100000101 以此为例，表示 0,3,90,3,90,3,9 个数均为奇数.appear−>101100101appear->101100101appear−>101100101 以此为例，
2021-08-12 10:42:40
214
原创 Codeforces Beta Round #13 B. Letter A 计算几何
B. Letter Asolution循环枚举作为横线的边，计算另外两边的距离是否为0来判断两线段是否相交；利用向量计算另外两边的夹角，判断其是否非钝；分别令横线切边分两线段，长度为 x,yx,yx,y,(x<y),(x<y),(x<y) 化 xy≥14\frac{x}{y}\ge\frac{1}{4}yx≥41 为 4x−y≥0.4x-y\ge0.4x−y≥0.精度需 1e−81e^{-8}1e−8,1e−61e^{-6}1e−6 可能更合适 doubledoubled
2021-08-07 01:19:29
165
原创 P4124 [CQOI2016]手机号码 数位DP
P4124 [CQOI2016]手机号码solutiondp[pos][state][is4][is8]dp[pos][state][is4][is8]dp[pos][state][is4][is8]state=0,表示当前数字没有连续state=0,表示当前数字没有连续state=0,表示当前数字没有连续state=1,表示当前连续两个相同数字state=1,表示当前连续两个相同数字state=1,表示当前连续两个相同数字state=2,表示三个连续数字state=2,表示三个连续数字state=2
2021-08-06 20:18:29
167
原创 Codeforces Beta Round #2 C. Commentator problem 计算几何 + 模拟退火
C. Commentator problemsolution找点一个点使得，该点与三个圆的切线夹角都相等。找点问题，用模拟退火.让每个 α\alphaα (6个) 都相等即找到D点.code/*SiberianSquirrel*//*CuteKiloFish*/#include <bits/stdc++.h>//#include<bits/extc++.h>#include<ext/rope>#include<ext/pb_ds/assoc_c
2021-08-05 10:08:02
249
原创 计算几何小知识
位似图形对应几何性质完全相同∏r=1n−1sinrπn=n2n−1\prod\limits_{r=1}^{n-1}sin\frac{r\pi}{n}=\frac{n}{2^{n-1}}r=1∏n−1sinnrπ=2n−1n∏k=1n(1+2cos2kπn)=3\prod\limits_{k=1}^{n}(1+2cos\frac{2k\pi}{n})=3k=1∏n(1+2cosn2kπ)=3∏k=1n−1(1+2coskπn)={0n=3k,(−1)n−1n=3k+1,(−1)nn=3k+2.
2021-08-04 20:35:25
159
原创 Codeforces Round #259 (Div. 1) D. Little Pony and Elements of Harmony FWT
D. Little Pony and Elements of Harmonysolutionei[u]=∑vei−1[v]∗b[f(u,v)].e_i[u]=\sum_{v}e_{i-1}[v]*b[f(u,v)].ei[u]=v∑ei−1[v]∗b[f(u,v)].设 cnt[u⊕v]=b[f(u⊕v)]cnt[u⊕v]=b[f(u⊕v)]cnt[u⊕v]=b[f(u⊕v)]，有：ei[u]=∑vei−1[v]∗cnt[u⊕v].e_i[u]=\sum_{v}e_{i-1}[v]*cnt[
2021-08-04 16:23:26
167
原创 2021牛客暑期多校训练营5 F - Finding Points 计算几何 + 三分
Finding Pointssolution三分 xxx， 三分 yyy 即可，确保 x,yx,yx,y 均在圆内部，暴力遍历所有角度，得答案.code/*SiberianSquirrel*//*CuteKiloFish*/#include <bits/stdc++.h>using namespace std;#define gcd(a,b) __gcd(a,b)#define Inv(x) quick_pow(x, mod - 2)#define Polynomial ve
2021-07-31 21:24:39
278
原创 2021牛客暑期多校训练营5 B - Boxes 期望（线性）
solution期望题，两种情况：1.1.1. 花费 CCC ，得知n个箱子中的黑白球个数.2.2.2. 不花费 CCC，将所有箱子打开.第二种情况很简单，因此我们展开考虑第一种情况：假设 f(x)f(x)f(x) 表示开启 xxx 个箱子时完成要求，所需的费用.f(0)=0.f(0)=0.f(0)=0.先将箱子进行排序，优先开启费用小的箱子，对于黑球，可能有 0,1,2,..,n0,1,2,..,n0,1,2,..,n 个，因此我们开启第一个箱子的时候，花费 w1w_1w1，如果此时游戏.
2021-07-31 19:48:42
194
原创 2021“MINIEYE杯”中国大学生算法设计超级联赛（2） I love max and multiply HDU-6971 SOSDP
I love max and multiplysolution/*SiberianSquirrel*//*CuteKiloFish*/#include <bits/stdc++.h>using namespace std;#define gcd(a,b) __gcd(a,b)#define lcm(a,b) (1ll * a * b / gcd(a, b))#define Polynomial vector<int>#define Inv(x) quick_pow(
2021-07-28 11:25:25
158
原创 2021牛客暑期多校训练营4 Sample Game 生成函数
Sample Gamesoltuionf(x)=∑i=0∞Pixif(x)=\sum_{i=0}^{∞}P_ix^if(x)=∑i=0∞Pixi ,( PiP_iPi 表示构成长度 >i>i>i 的序列的可能性)由于序列为不下降序列，例：1112233451 1 1 2 2 3 3 4 5111223345，Pi=∏∑cntp=ipkP_i=\prod_{\sum cnt_p=i} p^kPi=∏∑cntp=ipkf(x)=∏j=1∞∑k=0ipjkxk=∏j=
2021-07-28 11:14:22
300
原创 递推 数列
递推solutionan=2an−1+n2a_n=2a_{n-1}+n^2an=2an−1+n2an+n2+4n+6=2[an−1+(n−1)2+4(n−1)+6]a_n+n^2+4n+6=2[a_{n-1}+(n-1)^2+4(n-1)+6]an+n2+4n+6=2[an−1+(n−1)2+4(n−1)+6]an+n2+4n+6=2n(a0+6)=6∗2na_n+n^2+4n+6=2^n(a_0+6)=6*2^nan+n2+4n+6=2n(a0+6)=6∗2nan=6∗2n−n2−
2021-07-20 19:16:31
175
原创 2021牛客暑期多校训练营2-Girlfriend 计算几何
Girlfriendsolution∣AP∣≥k∣BP∣,(k>1)|AP|≥k|BP|,(k>1)∣AP∣≥k∣BP∣,(k>1)，在二维平面上称为阿波罗尼斯圆，此题三维即阿波罗尼斯球.设A(x0,y0,z0),B(x1,y1,z1),P(x,y,z)A(x_0,y_0,z_0),B(x_1,y_1,z_1),P(x,y,z)A(x0,y0,z0),B(x1,y1,z1),P(x,y,z)，我们有（下面均取等号位置），(x−x0)2+(y−y0)2+(z−z0)2=
2021-07-19 17:00:13
239
原创 P4213 【模板】杜教筛（Sum） 杜教筛
P4213 【模板】杜教筛（Sum）solution杜教筛模板题杜教筛模板题杜教筛模板题code#include <iostream>#include <cstdio>#include <map>using namespace std;typedef long long ll;const double eps = 1e-7;const int MOD = 1e9 + 7;const int N = 5e6 + 10;int prime[N
2021-07-19 16:39:33
146
原创 E. Nikita and Order Statistics Codeforces Round #488 by NEAR (Div. 1) FFT
solution求任意区间内元素小于 kkk 的个数，并按顺序输出 k=0..nk=0..nk=0..n 的区间数.记录前缀和 ans[i]ans[i]ans[i] 表示前 iii 个数里小于 xxx 的元素个数.区间 l,rl,rl,r 内小于 xxx 的元素个数为 ans[r]−ans[l−1]ans[r] - ans[l - 1]ans[r]−ans[l−1]，即 xans[r]−ans[l−1]x^{ans[r]-ans[l-1]}xans[r]−ans[l−1] 的系数.由于 −ans[.
2021-07-19 10:43:32
175
原创 Hash Function 2021牛客暑期多校训练营1 数论 + NTT
Hash Functionsolution任意选择该集合中的两个数，a，b，由于a%mod<>b%mod任意选择该集合中的两个数，a，b，由于a\%mod<>b\%mod任意选择该集合中的两个数，a，b，由于a%mod<>b%mod推出(a−b)%mod<>0.推出(a-b)\%mod<>0.推出(a−b)%mod<>0.因此只需要得到所有(c=a−b)的值，将c的因子，即(c%d==0)的d全部标记即可。因此只需要得到所有(
2021-07-18 11:16:48
107
原创 无限手套 分治 + NTT + 生成函数
solution∏k=0∞aixi2+bixi+1\prod\limits_{k=0}^{∞}a_ix_i^2+b_ix_i+1k=0∏∞aixi2+bixi+1考虑生成函数：考虑生成函数：考虑生成函数：=>∑k=0∞aik2xk+bixixk+xk=>\sum\limits_{k=0}^{∞}a_ik^2x^k+b_ix_ix^k+x^k=>k=0∑∞aik2xk+bixixk+xk=∑k=0∞aik2xk+∑k=0∞bixixk+∑k=0∞xk,(三项分别当作A
2021-07-17 09:23:32
314
1
原创 Visible Lattice Points 莫比乌斯反演
Visible Lattice Pointssolution三维图形，分成四个部分（x,y,z分别=0，x,y,z同时不等于0）,最后加上三条正半轴上的点三维图形，分成四个部分（x,y,z分别=0，x,y,z同时不等于0）,最后加上三条正半轴上的点三维图形，分成四个部分（x,y,z分别=0，x,y,z同时不等于0）,最后加上三条正半轴上的点Answer=∑d=1nμ(d)[(3+nd)∗(nd)2]+3Answer=\sum\limits_{d=1}^{n}\mu(d)[(3+\frac{n}{d}
2021-07-16 14:51:46
68
原创 毕业生的序列游戏 期望DP
毕业生的序列游戏solutiondp[i][j]表示i个a，j个ab的概率，添加a的概率为p,添加b的概率为qdp[i][j]表示i个a，j个ab的概率，添加a的概率为p,添加b的概率为qdp[i][j]表示i个a，j个ab的概率，添加a的概率为p,添加b的概率为qif i+j<k(如果放b，ab个数+i，即加上先前a的个数):if~i+j<k(如果放b，ab个数+i，即加上先前a的个数):if i+j<k(如果放b，ab个数+i，即加上先前a的个数):dp[
2021-07-15 10:10:42
116
原创 Alice和Bob赌糖果 牛客 概率DP（赌徒破产模型）
Alice和Bob赌糖果solutionf0=0,fn+m=1f_0=0,f_{n+m}=1f0=0,fn+m=1fi=(1−p)fi−1+pfi+1f_i=(1-p)f_{i-1}+pf_{i+1}fi=(1−p)fi−1+pfi+1=>fi+1−fi=fi−fi−1=>f_{i+1}-f_i=f_i-f_{i-1}=>fi+1−fi=fi−fi−1令fi=gifi+1=>ki=p1−(1−p)ki−1令f_i=g_if_{i+1}=>k_i=\
2021-07-14 18:55:37
297
原创 Xieldy And His Password 牛客 DP
Xieldy And His Passwordsolutiondp[i][j]表示第i位为止有多少%3=j的方案数dp[i][j]表示第i位为止有多少\%3=j的方案数dp[i][j]表示第i位为止有多少%3=j的方案数记前一个状态为tempi−1,记前一个状态为temp_{i-1},记前一个状态为tempi−1,s[i]==1时，tempi=(tempi−1∗2+1)%3s[i]==1时，temp_i=(temp_{i-1}*2+1)\%3s[i]==1时，tempi=(tempi−1∗2
2021-07-14 15:18:22
100
原创 组合数总结
卡特兰数给定nnn个000和nnn个111，它们将按照某种排序成长度为2n2n2n的序列，求它们能排列成的所有序列中，满足任意前缀序列中000的个数都不少于111的序列有多少个？（即从X走到Y路线不穿越，y=x的方案数）|~~~~~~~~~~~~~~~~~~~.Y(n,n)||||||.X(0,0)|--------------------------------------------------------------->Answer=C2nn−C2nn−1=C2nnn+1A
2021-07-14 11:05:03
240
原创 E - Counting Sequences 2019年ICPC 上海网络赛 指数型生成函数
Counting Sequencessolution序列[a1,a2,a3...,an],(ai∈[1,m])，问序列中ai为偶数的个数非奇有几种序列[a_1,a_2,a_3...,a_n],(a_i∈[1,m])，问序列中a_i为偶数的个数非奇有几种序列[a1,a2,a3...,an],(ai∈[1,m])，问序列中ai为偶数的个数非奇有几种考虑生成函数：考虑生成函数：考虑生成函数：T=(1+x22!+x44!+x66!+...)m2∗(1+x1!+x22!+x33!+...)m+1
2021-07-14 10:06:36
170
原创 Interesting Series 计蒜客 FFT + CDQ分治 +生成函数
Interesting SeriessolutionF1=1F_1=1F1=1Fn=a∗Fn−1+1(n≥2)F_n=a*F_{n-1}+1(n≥2)Fn=a∗Fn−1+1(n≥2)=>Fn+1a−1=a∗(Fn−1+1a−1)=>F_n+\frac{1}{a-1}=a*(F_{n-1}+\frac{1}{a-1})=>Fn+a−11=a∗(Fn−1+a−11)Fn−1a−1=an−1(1+1a−1)=>Fn=an−1a−1F_n-\frac{1}{a-1}=a^
2021-07-13 11:21:35
145
原创 TrickGCD HDU-6053 莫比乌斯反演 容斥
TrickGCDsolutionF(n)表示[gcd==kn]的个数，k属于Z+,F(n)表示[gcd==kn]的个数，k属于Z_+,F(n)表示[gcd==kn]的个数，k属于Z+,f(n)表示[gcd==n]的个数.f(n)表示[gcd ==n]的个数.f(n)表示[gcd==n]的个数.F(n)=∏i=1nainF(n)=\prod\limits_{i=1}^{n}\frac{a_i}{n}F(n)=i=1∏nnaiF(n)=∑n∣df(d)=>f(n)=∑n∣dμ(dn)F
2021-07-10 10:50:34
132
原创 导数卷积 牛客 NTT
f(x)=∑i=0n−1ai∗xif(x)=\sum\limits_{i=0}^{n-1}a_i*x^if(x)=i=0∑n−1ai∗xif(x)=∑i=0n−1ai∗xif(x)=\sum\limits_{i=0}^{n-1}a_i*x^if(x)=i=0∑n−1ai∗xig(x)=∑i=0n−1f(i)(x)∗f(n−i+1)(x)g(x)=\sum\limits_{i=0}^{n-1}f^{(i)}(x)*f^{(n-i+1)}(x)g(x)=i=0∑n−1f(i)(x)∗f(n−i+1
2021-07-10 08:29:49
306
原创 乒乓球 牛客 NTT
乒乓球solution注意：以下 k=i−j+1表示i,j中间元素的个数，i<j注意：以下~k=i-j+1表示i,j中间元素的个数，i<j注意：以下 k=i−j+1表示i,j中间元素的个数，i<j考虑拿出顺序，对于i,j来说，它们之间有k个数(0<k<n−1),因此有2!∗k!种，剩下n−k−1个球插空放入(k+3)(k+3)...(n−1)n，且总方案为n!，易得：考虑拿出顺序，对于i,j来说，它们之间有k个数(0<k<n-1),因此有
2021-07-09 18:36:35
238
原创 RXD and functions HDU-6061 NTT
RXD and functionssolutiongm=f(x−∑i=1mai),令A=∑i=1mai，则有，g_m=f(x-\sum\limits_{i=1}^{m}a_i),令A=\sum\limits_{i=1}^{m}a_i，则有，gm=f(x−i=1∑mai),令A=i=1∑mai，则有，gm=f(x−A)g_m=f(x-A)gm=f(x−A)=∑i=0nci(x−a)i=\sum\limits_{i=0}^{n}c_i(x-a)^i=i=0∑nci(x−a)i=∑i=
2021-07-09 15:27:21
142
原创 离散余弦变换和离散傅里叶变换
DCT（离散余弦变换）：DCT（离散余弦变换）：DCT（离散余弦变换）：X[u,m]=∑u=0N−1∑m=0N−1x[u,m]∗c[u,m]X[u,m]=\sum\limits_{u=0}^{N-1}\sum\limits_{m=0}^{N-1}x[u,m]*c[u,m]X[u,m]=u=0∑N−1m=0∑N−1x[u,m]∗c[u,m]其中，其中，其中，c[u,m]={1N，u=02Ncos((2∗m+1)uπ2N),u!=0c[u,m]=\left\{\begin{matrix}\sqrt{\
2021-06-23 10:44:23
1638
原创 Intersection HDU-5120 计算几何+容斥定理
Intersectionsolution相当于四个圆：aR,ar,bR,br相当于四个圆：aR,ar,bR,br相当于四个圆：aR,ar,bR,brS总=S(aR)∩S(bR)−S(aR)∩S(br)−S(ar)∩S(bR)+S(ar)∩S(br)S_总=S(aR)∩S(bR)-S(aR)∩S(br)-S(ar)∩S(bR)+S(ar)∩S(br)S总=S(aR)∩S(bR)−S(aR)∩S(br)−S(ar)∩S(bR)+S(ar)∩S(br)code/*SiberianSquirrel*/
2021-05-24 19:08:45
569
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人