- 博客(239)
- 收藏
- 关注
原创 手搓全连接层神经网络识别手写数字
需要对28x28的灰度数字图像进行数字识别。输入层 28x28个输入节点,分别代表像素点上的灰度大小输出层10个输出,分别代表0-1的概率使用sigmoid函数由于求导过于繁琐,因此这里仅添加了一层隐藏层。代码如下#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>using namespace std;char file[]="train-images-id
2021-11-28 23:39:28
775
1
原创 聪明钱因子(R语言实现)
聪明钱因子最近回到了csdn,打算记录一下自己在量化策略方面的学习历程。正好最近也快R语言的期末考试了,这几天就用R实现了一个简单的聪明钱因子,但不得不说R语言的读入数据是真的慢。。。根据研报上的信息,将α=∣pctchg∣volume\alpha=\frac{|pct_{chg}|}{\sqrt{volume}}α=volume∣pctchg∣作为一个中间变量,其中pctchgpct_{chg}pctchg为涨跌幅,volumevolumevolume为交易量,对每分钟的交易数据处理α\al
2021-05-03 21:19:34
1187
5
原创 #长期填坑# 一个奇怪的静态树上联通块处理技巧
第二次做CC的那道边上gcd的题目想出来一个奇怪的技巧 可以把理论复杂度从O(nw‾‾√logn+q2w‾‾√logn)O(n\sqrt{w} logn+q^2\sqrt{w}logn)变成O((n+q)w‾‾√lognw)O((n+q)\sqrt{w}lognw) 在线询问 支持可持久化 空间的话。。貌似有点大大。。和时间复杂度同阶 然而这个是理论复杂度。。。并查集的loglog和线段树的
2017-03-16 08:52:54
1292
原创 BZOJ4762: 最小集合
感谢度教。。。转化条件1 : 使得Or集为全集使得Or集为全集 其他不变一个合法集合S必不存在一个大小为|S|−1的子集T满足条件1一个合法集合 S 必不存在 一个大小为 |S|-1的子集T满足条件1 即一个集合合法必不存在一对长度和小于|S|的前缀与后缀使得Or集为全集即一个集合合法 必不存在一对长度和小于|S|的前缀与后缀使得Or集为全集然后考虑暴力 g[i][j][k]表示已经做了前
2017-03-13 12:10:16
1260
原创 BZOJ4768: 2555加强版之wxh loves substring
很显然的后缀平衡树 一开始以为要可持久化 发现根本不用。。 treap的常数要死人啊?寄刀片寄刀片 #include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<map>#include<cmath>using namespace std;#define ld long longchar
2017-03-10 21:55:27
1141
原创 Anarchy的解题报告
题目大意: 假设高斯定理在mm维空间成立 已知mm维空间所有整点电荷 aja_j 给出mm维空间下x,yx,y两点距离公式 以及x点在y点引发的电势公式 降序输出求前100个点的电势T≤51≤m≤18n≤3∗105T≤5 1≤m≤18 n≤3*10^5 时限7s7s考试的时候看这道题的时候 感觉这道题根本不可做啊??? 部分分很友好啊?除了一个5分的点 其他我都不会。。看
2017-03-10 15:49:42
669
原创 2017.3.1
首先这是一篇检讨 这场比赛我因为思考不充分便开始打T2 bug如雪球一般越滚越大 最后GG这场比赛我没有仔细阅读题面 而是看一眼题目 自信心爆膨 便开始刚T2 导致其他题目暴力没打 五个小时的时间里 我的T2遇到了很多问题 而我并没有从根本上去质疑我的算法 而是不断地打着补丁 并作死。。。这场我觉得是对我最有启发的一场比赛 1、仔细阅读所有题目之后再开始 2、思考时需要保证
2017-03-01 17:15:17
595
原创 2017.7.22模拟赛被虐记录
感觉自己要完 成功被虐了两倍分。。 其实故事是这样的 看了第二题 发现是个SB树分 然后就没有然后了 突然Manchery说自己的代码会被卡常 然后我拿了大样例 没开优化跑了0.8s 时限2s 感觉很稳啊 看了第一题 完全没思路 看了看大样例 我日 SB题。。。 但是要特判几个hack点。。。 第三题 弃疗 暴力都不会。。。 然后。。。本机跑0.8s的我在评测机上2s TL
2017-02-22 18:25:57
554
原创 补2016.2.21
貌似昨天的又忘了写了? 做了一上午的题 下午我在干嘛?好像是在推一个生成函数的题目。。然后发现二项式定理和暴力乘法答案不一样 弃疗。。。。感觉自己已经是一条咸鱼了晚上TC 成功贡献一记Hack 有小哥叉了我的随机化的代码是闹哪样?
2017-02-22 18:08:43
553
原创 补2016.2.20
来 我来补一下昨天 昨天讲了一些组合数学的内容 感觉好像没什么营养的样子? 好像这些都是说给高一的人听的吧?多项式什么的感觉好像没有劼老师讲的愉悦啊? 最后的混吃等死。。。
2017-02-21 11:03:16
545
原创 2017.2.18模拟赛总结反思
第一次模拟赛被绍一大爷虐惨了。。。。看T1 woc这不是SB虚树么?打上ST表之后发现自己丝薄了。。 听着别人键盘噼里啪啦 蜜汁心慌一小时后我才发现T1可以O(n32)O(n^{\frac{3}{2}})的做 然后感觉复杂度太高 看着旁边Manchery已经码起来了 感觉很害怕 问了一下复杂度 居然和我一样 于是我也不虚了然后调出大样例手造数据发现自己代码跑的还蛮快的 之后开始刚第三
2017-02-18 18:44:44
527
原创 BZOJ4695: 最假女选手
我好像被卡空间了啊? 惨啊.. 代码好丑啊QAQ 考虑更新 类似VEB Tree的记录当前最大值 不过这样子还是不容易更新最小值的 于是记录一个次大值 和最大值个数便于更新 对于最大值用类似方法维护#include<cstdio>#include<iostream>#include<cstring>#include<cmath>using namespace std;char
2017-02-17 11:35:51
924
原创 Hackerrank How many substrings
好难啊QAQ 今天又去想了一下 然后向敦爷求证了一下 发现思路还是对的代码太难实现了。。 LCT维护last并在Segment Tree上更新。。。 感觉自己被日死了 最后看着板子打了一发。。 先离线 考虑每个串在最后出现的开始位置有贡献 然后就是一个区间加 区间减的操作 在Access的时候更新线段树 时间复杂度O(nlognlogn)O(nlognlogn)你看800+sub
2017-02-16 20:40:51
706
原创 HDU5766 Filling
Burnside+插头DP 旋转有四个群 不动 旋转 π2\frac{π}{2} ππ 3π2\frac{3π}{2} 显然要染色在某置换下不动则每一块结构相同 考虑分别计算 显然 旋转π2\frac{ π}{2} 与旋转 3π2\frac{ 3π}{2} 答案相同我们只需要考虑旋转π2\frac{ π}{2}的情况 此时将n*n的格子分成[n+12][\frac{n+1}{
2017-02-15 15:40:53
343
原创 Stirling数学习笔记
劼爷上的课现在才去整理…第一类Stirling数 s(n,m)表示n个元素组成m个圆排列s(n,m)表示n个元素组成m个圆排列 由以上定义我们可以得出递推公式:由以上定义我们可以得出递推公式: s(n,m)=∑n−1i=0s(n−i,m−1)∗Ci−1n−1∗(i−1)! s(n,m)=\sum_{i=0}^{n-1}s(n-i,m-1)*C_{n-1}^{i-1}*(i-1)! 以及
2016-12-18 15:39:16
1921
原创 Hackerrank XOR Subsequences
求连续子串的Xor和出现次数最多的Xor和 Xoryi=xAi=(Xorx−1i=1Ai)Xor(Xoryi=1Ai)Xor_{i=x}^{y}Ai=(Xor_{i=1}^{x-1}Ai) Xor (Xor_{i=1}^{y}Ai) 最后发现其实就是个卷积形式 上FWT即可 中间处理过程可能会爆int 所以加上取模会免掉不必要的麻烦#include<cstdio>#include<iost
2016-12-15 18:19:41
762
原创 FFT优化进制转换 十进制转二进制
其实进制转换都能这么做.. ∑ni=0xi∗Basei\sum_{i=0}^{n}x_i*Base^i ∑n2−1i=0xi∗Basei+Basen2∗∑n2i=0xi+n2∗Basei\sum_{i=0}^{\frac{n}{2}-1}x_i*Base^i+Base^{\frac{n}{2}}*\sum_{i=0}^{\frac{n}{2}}x_{i+\frac{n}{2}}*Base^i
2016-12-14 19:49:52
1834
原创 Hackerrank Hard Homework
题目链接: https://www.hackerrank.com/contests/w26/challenges/hard-homework 题目大意: 给定nn 求正整数x,y,zx,y,z满足x+y+z=nx+y+z=n 最大化sin(x)+sin(y)+sin(z)sin(x)+sin(y)+sin(z)一开始往死里搞没拆出来.. 搞了一个O(nlogn)O(nlogn)的做法结果死
2016-12-12 19:25:16
488
原创 ZOJ2674 Strange Limit
欧拉函数….#include<cstdio>#include<cmath>#include<iostream>using namespace std;int p;int Fai(int x){ int t=sqrt(x),res=1; for(int i=2;i<=t&&i<=x;i++) if(x%i==0) { x/=i;
2016-12-11 20:35:45
779
原创 Ural1540 Battle for the Ring
区间DP+SG函数看到代码淦过去的时候我的内心是崩溃的… 明明O(n^4)啊。。#include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<cmath>#include<algorithm>using namespace std;char c;inline void read(int&a)
2016-12-11 19:22:25
342
原创 Hackerrank Random Number Generator
题目链接: https://www.hackerrank.com/challenges/random-number-generator-1 题目大意: Pi>=0P_i>=0 且1=Σni=1Pi1=Σ_{i=1}^{n}P_i 最大化 Σni=1Pi∗(1−Pi)∗iΣ_{i=1}^{n}P_i*(1-P_i)*i一开始的感觉很对 但是发现拉格朗日出现负数就不会处理 后来问了uo
2016-12-09 17:10:55
566
原创 BZOJ3270: 博物馆
类似之前usaco的那道题。。 可能简单一点 少几个变量。。#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<iostream>using namespace std;char c;inline void read(int&a){a=0;do c=getchar();while(c<'0
2016-12-08 20:32:04
542
原创 Hackerrank Connect the country
题目链接: https://www.hackerrank.com/challenges/connect-the-country由于n的整数拆分不超过6000 我们可以考虑用每个整数拆分表示一个状态 预处理每个状态之间的转移 然后暴力跑复杂度不是很会证明 反正跑得过去#include<cstdio>#include<map>#include<cmath>#include<cstring
2016-12-08 19:33:20
432
原创 BZOJ1778: [Usaco2010 Hol]Dotp 驱逐猪猡
设概率转移矩阵为An∗nA_{n*n} 在爆炸概率矩阵PP 其中P1,iP_{1,i}为 在点 ii爆炸概率 设行向量T1∗n=T_{1*n}={1,0,0..1,0,0..}(1∗n) (1*n) 设B=p/qB=p/q 即在每一时刻炸弹爆炸概率 有P1∗n=B∗T1∗n∗Σ∞j=1(1−B)j−1∗Aj−1n∗nP_{1*n}=B * T_{1*n}*Σ_{j=1}^{∞}(1-B)^{
2016-12-07 16:56:52
643
原创 AFO
打出这三个字母或许需要莫大的勇气 因为从此OI也只能成为回忆了 爆炸不能找理由 只是因为自己弱看到因为自己两题的两个小错误 到手的200分没了 幻想最终破灭 不甘心 可也只剩下不甘心了 自己不成熟的心态使我沦落到现在的成绩 我无话可说请假回来后和父母交换了一下意见 还是去冲一下中科大吧 虽然很不情愿 但也没办法 只能怪自己弱 考挂不如狗 博客置顶帖中的意气扬扬现在在我眼中是那
2016-11-21 19:00:05
787
1
原创 HDU5852Intersection is not allowed!
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5852https://en.wikipedia.org/wiki/Lindström–Gessel–Viennot_lemma 行列式的sgn刚好是容斥系数 糖教说这个可以应用到区间图所以直接建行列式然后做成对角线就可以了#include<cstdio>#include<iostream>#i
2016-11-18 15:18:19
447
原创 HDU5608 function
设 G(i)=i2−3i+2G(i)=i^2−3i+2 有莫比乌斯反演可得出F(i)F(i)关于G(i)G(i)的关系 带入得 ∑ni=1F(i)=∑nd=1u(d)∑n/di=1G(i)∑_{i=1}^nF(i)=∑_{d=1}^nu(d)∑_{i=1}^{n/d}G(i) 后面可以O(1) 前面就是经典的杜教筛了#include<cstdio>#include<iostream>
2016-11-17 12:52:40
759
原创 51Nod 1220 约数之和&&PE439
∂(i)∂(i)表示i的约数和 ∑ni=1∑nj=1∂(ij)∑_{i=1}^{n}∑_{j=1}^{n}∂(ij) =∑ni=1∑nj=1∑w|i∑v/jwv∗∑d|w,d|(i/v)u(d)=∑_{i=1}^{n}∑_{j=1}^{n}∑_{w|i}∑_{v/j}wv*∑_{d|w,d|(i/v)}u(d) 我们尝试把d提前 会发现有 =∑nd=1d∗u(d)(∑(n/d)i=1∂(i
2016-11-17 07:46:22
911
原创 SPOJ DIVCNT2
题目链接: http://www.spoj.com/problems/DIVCNT2/根据rzz的课件 可以分三段做度教搞出来一种分一次做的方法 看起来很神的样子#include<cstdio>#include<iostream>#include<vector>#include<cstring>#include<cstdlib>#include<algorithm>#include<
2016-11-16 16:19:20
863
原创 Hackerrank Fibonacci Numbers Tree
题目地址 https://www.hackerrank.com/challenges/fibonacci-numbers-tree 大意: 子树按深度加一段连续的Fib数列 询问两点之间的权值和注意到对于任意数列FF Fi=Fi−1+Fi−2F_{i}=F_{i-1}+F_{i-2} 当a>0,b>0a>0,b>0我们有 Fa+b=Fa−1∗Fibb+Fa∗Fibb+1
2016-11-16 09:49:54
447
原创 Hackerrank SG 函数小结
今天去做了一下Hackerrank上的博弈题 感觉某些Hard有点偏容易?Simple Game: https://www.hackerrank.com/challenges/simple-game 思路:分k=2 k=3 k>=4讨论 计算出SG函数之后直接暴力DPMove the Coins: https://www.hackerrank.com/challenges/move-t
2016-11-14 20:08:27
535
原创 支配树(Dominator tree)学习笔记 及HDU4694Important Sisters
前置技能: LCA (Tarjan)今年冬令营上讲的东西现在才学 惭愧惭愧。。 首先做出图G的dfs树T 其中S为起点 定义T中节点 u 的半必经点sdom[u] 为 u的祖先中能经过若干条非树边到达u的(除了sdom[u]直接引出的边 其余边必须为非树边) 深度最小的节点 则对于任意一条边 (v->u) 这里我们默认sdom idom 记录的是dfn
2016-11-08 18:41:12
2420
原创 后缀平衡树学习笔记
给定一个空串SS 操作1 代表在SS前加入一个字母使之成为新SS 操作2 代表在询问在当前SS中有多少连续子串等于给定串TT假设我们已经有了串SS的后缀平衡树 插入一个字母cc 我们用SiS_i 代表原串SS从第ii个字符开始的后缀 则后缀 cScS与 任意一个后缀 SiS_i的大小关系可以用 c与SiS_i 的第一个字母的大小关系 以及 后缀S1S_1与后缀Si+1S_i+1
2016-11-06 19:48:53
1423
原创 HDU5823color II
怎么现在FWT啥的都烂大街了呀 将集合S分为两个没有连边的集合A BF[i][S]表示集合S能否用i种颜色染色 则有 F[i][S]=[∑A∈SF[i−1][A]∗F[1][B]>0]F[i][S]=[∑_{A∈S}F[i-1][A]*F[1][B] >0] 多美的一个卷积形式啊#include<cstdio>#include<cstring>#include<cstdlib>#in
2016-11-02 19:17:38
490
原创 51NOD1525 重组公司
这SB题我居然会WA?#include<cstdio>#include<iostream>#include<cstring>#include<algorithm>#include<cstdlib>#include<cstdio>using namespace std;int n,q;char c;inline void read(int&a){a=0;do c=getchar()
2016-10-29 12:19:57
500
原创 HDU5751Eades
先处理连续一段中最大值出现位置然后用FFT计算出各个最大值出现次数时的线段树 累加答案#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<complex>#include<algorithm>#include<iostream>#include<queue>#include<vector
2016-10-28 11:47:23
481
原创 HDU5808 Price List Strike Back
分治 类似ZJOI2016的最短路分治思想 妈的有毒 极大值有毒#include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<queue>#include<algorithm>#include<stack>#include<vector>using namespace std;char
2016-10-28 08:14:52
411
原创 HDU 5910Advanced Traffic System
由于边数太多且大多相同我们考虑用点去扩展 类似想法就是把原来的边看做点 原来的点看做一群边 把原来这个点附近的边代表的点相连 把加边操作挂上线段树上 然后类似Dij做一遍 并查集处理做过的点#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<algorithm>#include<queue
2016-10-27 13:43:55
440
原创 HDU5909 Tree Cutting
裸的树上FWT FWT取模很方便#include<cstdio>#include<iostream>#include<cstdlib>#include<cmath>#include<cstring>#include<vector>#include<queue>#include<stack>#include<algorithm>using namespace std;char c;
2016-10-26 15:31:36
426
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人