孤鸿子_-CSDN博客

原创 HDU 4609 3-idiots(多项式计数 + FFT卷积优化)

题目链接3-idiots分析这是第一道FFT,本想拿模板试试手的谁知道直接来了一个错误的模板……debug好久,具体分析就不说了，比较像构建生成函数，然后计数，中途用了卷积优化kungbin大佬的blog写的很好我先用内建complex交了一发，1300ms 大概比手写慢500ms(233333)AC codeProblem : 4609 ( 3-idiots ) Judge Status

2017-08-21 22:11:58 428

原创 HDU 6143 Killer Names(容斥原理)

题目链接Killer Names分析记 f(x)f(x) 为长度为n的字符串恰好用　xx 个字符填充有多少种取法，用容斥原理f(x)=xn+∑i=1n−1(−1)i(xi)(x−i)nf(x) = x^n + \sum_{i=1}^{n-1}(-1)^i\binom{x}{i}(x-i)^n 暴力计算 n2lognn^2logn,遗憾的是我居然超时了，所以我先暴力计算了　xnx^n 打了

2017-08-18 17:46:51 435

原创 HDU 6134 Battlestation Operational (mobius +前缀和)

题目链接Battlestation Operational计算f(n)=∑i=1n∑j=1i⌈ij⌉[(i,j)=1]f(n)=\sum_{i=1}^n\sum_{j=1}^i\lceil \frac{i}{j} \rceil[(i,j)=1]分析赛后看q神给了个简洁的题解才发现，自己走错方向了，计算前缀和要简单的多啊…记 g(i)=∑j=1i⌈ij⌉h(i)=∑j=1i⌈ij⌉[(i,j)

2017-08-17 23:41:57 432

原创 deeplearning Note : Practical aspects of Deep Learning

作者: Dylan_frank(滔滔)这是吴恩达 coursera Deep Learning Specialization 的第二门课程《Improving Deep Neural Networks: Hyperparameter tuning, Regularization and Optimization》的第一周笔记，主要讲在实现神经网络过程中所遇到的问题，和处理方法，具体来说是3点初

2017-08-17 01:18:44 3017

原创 HDU 6118 度度熊的交易计划 (最小费用最大流模板题)

题目链接度度熊的交易计划分析这是一个最小费用最大流的模板题，建边也很直接，唯一需要注意的是当跑出来的费用为正(dist[t]>0)的时候就不要跑了，当时就是这个地方弄错了，还以为我的模板错啦…..AC code/*Problem : 6118 ( 度度熊的交易计划 ) Judge Status : Accepted*RunId : 21751738 Language : G++

2017-08-16 13:56:04 334

原创 HDU Just do it (找规律暴力)

题目链接Just do it分析简单打一下表，我们会发现，在第　mm 次，第　xx 项对第　yy 项的贡献是　(y−x+m−1m−1)\binom{y-x+m-1}{m-1} 由于是异或所以只与这个组合数是否为奇数有关，由lucus定理　这个数是奇数 ⇔(y−x+m−1 and m−1)=1\Leftrightarrow (y-x+m-1\ and\ m-1)=1 为了方便后面的计算，我们简单

2017-08-16 13:07:42 352

原创 HDU 6121 Build a tree(找规律+implement)

题目链接HDU 6121 Build a tree分析规律是很好找的，就是实现比赛的时候感觉很难，没有去敲. 我用的是题解的方法,从根部向下递归，因为只有k个儿子，而每个儿子只有个是不满的，并且在这棵不满的树左边的儿子一定都是 LL 层的满树，而右边的节点数都是 L−1L-1 层的满树这个可以自行画图找规律.而对于满树来说计算他对答案的贡献是很简单的，这个规律很强，也可以看我的代码.那么现在的问题

2017-08-15 22:59:05 294

原创 UVA 10652 - Board Wrapping(凸包模板题)

题目链接10652 - Board Wrapping分析凸包模板题紫书272AC code#include<bits/stdc++.h>#define pb push_back#define mp make_pair#define PI acos(-1)#define fi first#define se second#define INF 0x3f3f3f3f#define INF6

2017-08-14 16:10:25 306

原创 codeforces 839 E. Mother of Dragons(最大团)

题目链接E. Mother of Dragons分析这个题的算法实现其实不难，难的是得出其中的结论，题解可以在这里看到 http://codeforces.com/blog/entry/53815AC code#include<bits/stdc++.h>#define pb push_back#define mp make_pair#define PI acos(-1)#define f

2017-08-13 20:11:30 498

原创 HDU 6109 百度之星2017　1002　数据分割(并查集+set)

题目链接百度之星2017　1002分析个比赛的时候想到了并查集+set的做法，只是评测实在太慢了，也没仔细想set处理合并的复杂度，赛后发现原来直接暴力合并set就行了　(不能分析复杂度是多少)AC code//Problem : 6109 ( 数据分割 ) Judge Status : Accepted//RunId : 21701687 Language : G++ Aut

2017-08-13 13:27:06 693

原创 HDU 6085 Rikka with Candies (暴力压位)

题目链接HDU 6085 Rikka with Candies分析其实很容易想到一个有技巧的暴力方法,我们可以这样办首先对于每一个 AiA_i 我们考虑比 AiA_i 大的部分 BB ，对于这部分 B,Ai%B=AiB,A_i\%B = A_i ,所以我们通过计算前缀和可以在 O(1)O(1) 计算出这部分值，因此我们先不考虑这部分值,对于比 AiA_i 小的部分 BB 我们反过来对每个 BiB_i

2017-08-09 14:54:50 415

原创 python语言学习:列表展开

问题提出这里遇到一个问题，比如我有两个列表 a=[1,2,3] b=[4,5,6]然后我们想做这样一个事情我们想把a中的元素加到b中，我们知道如果使用 b.append()a =[1,2,3]b = [4,5,6]b.append(a)bOut[37]: [4, 5, 6, [1, 2, 3]]加入的是嵌套链表,解决方法我们可以使用 b.extend(a)这样就可以将a的逐次迭代的

2017-08-07 21:42:03 595

原创 LA 3263 - That Nice Euler Circuit(计算几何基础题目　欧拉公式)

题目链接3263 - That Nice Euler Circuit分析这是紫书上的模板题，唯一有点思维的地方是用到欧拉公式,关于欧拉公式的简单证明可以看这篇blogAC code#include<bits/stdc++.h>#define pb push_back#define mp make_pair#define PI acos(-1)#define fi first#define

2017-08-07 16:45:07 485

原创 codeforces 833B. The Bakery(dp+线段树)

题目链接codeforces 833B. The Bakery分析这道题的dp方程不难，关键在于如何优化dp状态转移dp[i][j]=max1≤x≤j(dp[i−1][x]+c(x+1,j))dp[i][j] = max_{1 \le x \le j}(dp[i-1][x]+c(x+1,j))dp[i][j]dp[i][j] 表示 ii 个盒子装 1−j1-j 的时候的最优值c(i,j)c(i,

2017-08-06 14:51:38 509

原创 HDU 6069 Counting Divisors

题目描述HDU 6069 Counting Divisors分析设 i=pr11pr22…prkki = p_1^{r_1}p_2^{r_2}\dots p_k^{r_k},则 d(i)=∏ki=1(ri+1)d(i)=\prod_{i=1}^k(r_i+1),因此我们有d(ik)=∏ki=1(ri∗k+1)d(i^k)=\prod_{i=1}^k(r_i*k+1),显然用区间筛法，可以做到在 O(

2017-08-04 14:21:05 374

原创 HDU 6052 To my boyfriend

题目链接HDU 6052 To my boyfriend分析我这里用的是 n3logmn^3logm的算法。我们可以吧期望问题转化为计数问题，对于任意一个子矩阵来说，对于同种颜色的格子我们只计算从左到右，从上到下的第一个，那么我们反过来对于每个格子的贡献其实就是这个格子的颜色作为从左到右从上到下第一个的时候这样的矩阵有多少个。对于格子 (i,j)(i,j) 颜色为 cc 我们先求出从 1 to i1

2017-07-28 19:45:25 642

原创 HDU 6053 TrickGCD(莫比乌斯反演+分块)

题目链接TrickGCD分析题意其实很简单，反过来其实只要求得 gcd(bi)=1gcd(bi)=1　的个数,然后再用总的组合数减去就行了而对于求 gcd(bi)=1，bi<=Aigcd(b_i)=1，b_i <=A_i 由mobius公式 ∑x=1min{Ai}μ(x)∏i=1nfloor(Ai/x)\sum_{x =1}^{min\{A_i\}} \mu(x)\prod_{i=1}^nflo

2017-07-27 22:15:33 937

原创 <<python 核心编程>>第一章　正则表达式

基本使用方法match与searchimport rem = re.match('foo|bar','foo')print(m.group())footype(m.group())strm.group() 返回匹配的字符串pattern = 'foo|bar'm = re.match(pattern,'foobarfooooooook')print(m.group())foore.mat

2017-07-27 20:32:42 594

原创 HDU 6035 colorful tree

题目链接HDU 6035 colorful tree分析参考首先，我们单独考虑每种颜色，对于颜色 cc 　来说，他对答案的贡献就是包含这种颜色的路径条数，反过来想就是　n∗(n−1)/2−（不经过这种颜色的路径条数）n*(n-1)/2 - （不经过这种颜色的路径条数）,怎么求呢?我们以颜色　cc 来对树分块，那么　就是用总路径条数　n∗(n−1)/2n*(n-1)/2 剪掉　每一块的路径数目.不过我

2017-07-26 14:23:41 488

原创 codeforces D. Misha, Grisha and Underground(LCA)

题目链接D. Misha, Grisha and Underground分析LCA 模板题对于三个顶点　a,b,ca,b,c 因为，我们只需要找到他们的公共点，然后，公共点到各自距离子和的最大值就是ans u公共点怎么找呢？三个顶点会产生３个ＬＣＡ，一定有两个顶点的公共祖先是一样的，另外一个就是公共点，因此只需要对３个顶点异或一下就好了(还有这种操作~~~~)AC codehttp://cod

2017-07-25 22:16:28 548

原创 pandas 数据处理

Pandaspandas is a Python library for data analysis. It offers a number of data exploration, cleaning and transformation operations that are critical in working with data in Python.pandas build upon nu

2017-07-18 19:01:57 3428

原创 codeforces 830C Bamboo Partition

题目链接codeforces 830C Bamboo Partition分析首先我们需要知道什么样的d是满足的.对于任意的 ii ∑[a[i]+d−1d]−a[i]≤k\sum [\frac{a[i]+d-1}{d}] - a[i] \le k因此 d≤k+∑a[i]d \le k+\sum a[i]注意条件项中的 [a[i]+d−1d][\frac{a[i]+d-1}{d}] 在这中

2017-07-15 20:53:55 1123

原创 codeforces 830B. Cards Sorting

题目链接B. Cards Sorting分析这是一道很脑洞的玄学题.我们可以发现提取数字的顺序一定是从小到大的，每个数字需要用1次+转移到后面的次数。那问题的关键就是计算转移到后面的次数。我们可以发现每个数字转移到后面的次数等于在它后面比他小的次数+1,因此我们可以记录每个数字(设当前数字为 ii )出现的位置，从小到大排列，记录最后出现的位置为 pp ,用的次数为 (dis)那么 ii 后面的数字

2017-07-15 16:31:41 1471 3

原创 ISLR线性模型选择与正则化

ISLR线性模型选择与正则化ISLR线性模型选择与正则化zouzhitao2017年7月13日这一章主要讲进行特征选择(feature selection)的方法:子集选择压缩估计降维法子集选择方法最优子集选择最优子集选择法是一种穷搜( exhaustive search)策略，显然会发生维数灾难算法:记不含预测变量的模型为 $M_0$ .for k = 1:p拟合 \(\bin

2017-07-13 21:54:43 2249

原创《ISLR》交叉验证与自助法

注意本来是用r markdown写的所有代码放到Rstudio里面都可以运行，无奈csdn不支持这种东西理论所谓理论不过简单理解，并没有什么强力的证明验证集方法验证集方法，在Andrew的课里面涉及到过，就是将样本一部分作为训练集，一部分作为测试集，这种方法有两个缺点受训练集的影响波动性很大由于被训练的观测很少，统计方法表现不好留-交叉验证(LOOCV)这种方法用伪代码描述就是for i =

2017-07-12 19:51:13 1918

原创 R语言典型相关分析

参考资料《统计建模与R软件》典型相关的数学模型设 X=(X1,X2,…,Xp)T,Y=(Y1,Y2,…,Yq)TX=(X_1,X_2,\dots ,X_p)^T , Y=(Y_1,Y_2,\dots , Y_q)^T 为两条随机向量, 我们希望找到向量 a,ba,b 使得 U=aTX,V=aTY,ρ(U,V)达到最大U=a^TX,V=a^TY,\rho(U,V)达到最大，由于这样的向量有多组，我们加

2017-07-11 14:36:00 7846

原创 R语言主成分分析总结

简单总结R语言PCA相关函数这里是数据集 year X1 X2 X3 1951 1 -2.7 -4.3 1952 -5.3 -5.9 -3.5 1953 -2 -3.4 -0.8 1954 -5.7 -4.7 -1.1 1955 -0.9 -3.8 -3.1 1956 -5.7 -5.3

2017-07-10 16:56:28 28137 6

原创 R 统计推断

置信区间例有如下数据集求均值的置信区间 50,50,56,51,49,53,47,52,53,53,49,53,55,48,50,55> cof_interval <- t.test(x,conf.level = 0.95) #t.test 对总体均值的区间估计> cof_interval One Sample t-testdata: xt = 77.493, df = 15,

2017-07-08 22:57:55 613

原创 codeforces 819 B(区间加线性函数值)

问题描述首先我们考虑这样一个问题。给你一个操作 add [l,r] k(x-l)+b 对于 x∈[l,r],d[x]+=k(x−r)+bx \in [l,r],d[x]+=k(x-r)+b由于这个不是常数，我们没法用线段树或者树状数组来做，但是我们可以想到，如果是常数，我们一定可以这样 d[l]+=b,d[r+1]−=b;d[l]+=b,d[r+1]-=b; 然后都操作完成之后d

2017-07-01 00:50:36 663

原创 codeforce#420 E. Okabe and El Psy Kongroo(图论+矩阵快速幂)

题目链接E. Okabe and El Psy Kongroo分析首先对于坐标为 (x,y)(x,y) 的点,设 f(x,y)f(x,y) 表示从原点到 (x,y)(x,y) 的路径不难得到递推公式 f(x,y)=f(x−1,y−1)+f(x−1,y)+f(x−1,y+1)f(x,y) = f(x-1,y-1)+f(x-1,y)+f(x-1,y+1),很明显这和把每个点建图然后能达点建边构造的图是一

2017-06-27 16:00:19 559

原创 hihocoder 1305 区间求差(经典问题)

题目链接区间求差题目分析先来看一个经典问题，给你nn个区间，这些区间共覆盖了多少整点，很显然这些点将整条线段分为2n+12n+1段那么我们只需要计算每一段是否被覆盖就好了，怎么判断有没有被覆盖呢，很显然当当前段的左边的点中左端点个数比右端点个数多就行了.用cnt 作为计数器，遇到左端点+1，右端点-1,当cnt>0说明当前段会被计算，那么这道题就很简单了，对A与B分别计数，所要求得就是包

2017-05-29 22:25:58 672

原创 2017西安全国邀请赛总结

打完西安邀请赛感觉心有点被打垮了，最近几天有点颓废，加之期末不想刷题.思索再三感觉还是应该总结一下. 这一次拿了铜奖(参与奖)。心有不甘啊，感觉5题是怎么都有的. 一上来我就把两道签到题写了。然后看了一下榜，感觉A，C都是全场题，然后队友去开了C，我去搞了A推出来发现就是一个公式，上去交了6次，都是在边界的时候被卡了，这是已经过去3小时了等我过了以后，找了下队友做的C题发现了一个bug，交上

2017-05-26 11:07:41 1022

原创 hihocoder 1290 面试题(dp)

题目链接1290分析水dp 注意初始状态#include<bits/stdc++.h>#define pb push_back#define mp make_pair#define fi first#define se second#define INF 0x3f3f3f3fusing namespace std;typedef long long i64;typedef pair<

2017-05-25 19:57:03 425

原创 2017 江苏省赛JSCPC H(贪心,prim)

题目描述给你一棵树,两两节点间有一个距离cic_i,现在要重建地图，任意两个顶点间的花费是他们之间的最短路,建最大生成树求花费.分析这题比赛的时候没有做出来,赛后看题解原来很简单,想象prim的建树过程每次找剩余边集合中的最大距离.显然我们首先找到直径v0,v1v_0,v_1肯定作为第一条边然后照剩余点,对于剩下的每个点他要莫加到v0,要莫加到v1v_0,要莫加到v_1,因为可以证明他到任意点的距离

2017-05-15 22:46:04 846

原创 pandoc　md　文件转化pdf(中文字体)

教程在网上一搜就有不过有些已经陈旧了．先做出以下修改 1. 下载pandoc$ sudo apt-get install pandoc 下载xelatex引擎$ sudo apt-get install texlive-xetex texlive-latex-recommended texlive-latex-extra 配置字体模板\documentclass[$if(fonts

2017-05-12 11:21:44 3576 3

原创 UVA 11542 - Square(异或方程组消元经典题目)

题目链接11542 - Square分析每一个元素有选或则不选两种状态，因此设第i个元素为vi,xii个元素为v_i,x_i表示每个元素要不要选共nn个变元，由于我们可以先求出500500以内的素数,将选出来的数进行素因子分解后有pk1xj1+k2xj2+…,ki∗xji1p_1^{k_1x_{j_1}+k_2x_{j_2}+\dots,k_i*x_{j_i}}指数部分应该mod　2　为0的这样就会

2017-05-11 21:01:15 572

原创 UVA 10828 - Back to Kernighan-Ritchie(线性方程组jordan消元)

题目链接10828 - Back to Kernighan-Ritchie分析白书经典题目,每一个顶点用期望的线性性质列出方程.不过要注意特殊情况以及结合具体含义来分析,用gauss_jordan消元(转化为阶梯阵)之后.若有一行A[i][i]=0并且A[i][n]＝0,则说明此点不可被访问0次.继续若有矛盾方程A[i][i] = 0而A[i][n]不等于0则说明A[i][i]无解即不可能访问到终太

2017-05-10 22:58:56 383

原创 zoj 3954 Seven-Segment Display(暴力枚举)

题目链接Seven-Segment Display 分析预处理出所有解(29∗7!2^9*7!)，然后在线O(lgn)O(lgn)查询．具体来说,因为我们可以先枚举出7!7!种排列(交换情况)，对于每种交换由于查询的数码一定是9个数共292^9种情况，我们也可以枚举出来.最多需要7*9 = 63位，用一个long long保存就ｏｋ了.AC code#include <cstdio>#incl

2017-05-10 20:36:26 421

原创 zoj 3596Course Selection System(dp,01背包)

题目描述给你n(≤500)对数(0≤Hi()≤10000,0≤Ci≤100)n(\le 500)对数(0\le H_i()\le 10000,0\le C_i\le 100),让你选出ｍ个数最大化这个东西, 分析设x=∑mi=0Hxi,y=∑mi=0Cxi,则原等式可以化为(x−y)2−2∗y2x = \sum_{i = 0}^m H_{x_i},y = \sum_{i =0}^m C_{x_i}

2017-05-10 20:28:37 390

原创 HDU 4587 TWO NODES(割点模板题)

题目描述计算这个东西，即求删去两个点之后，的最大连通分量数. 3≤n,m≤5e33\le n,m\le 5e3分析其实就是一道模板题,枚举每一个要删去的点，然后计算剩余点的连通分量数，删掉一个对ans贡献最大的，即删掉它之后能够增加最大的连通分量数目.AC code#include <cstdio>#include <iostream>#include <vector>#include

2017-05-09 00:05:52 440

空空如也

空空如也