2017年09月_孤鸿子_

原创 HDU 5528 Count a * b 2015 长春现场赛(积性函数)

题目链接Count a * b分析这是很有意思的积性函数问题反过来定义 h(m)=m2−f(m)=∑a,b[a∗b%m=0]=∑a=1,bmgcd(a,m)|b=∑a=1mgcd(a,m)=∑d|md∗ϕ(m/d)\begin{align}h(m) &= m^2 - f(m)\\&=\sum_{a,b}[a*b\%m=0]\\&=\sum_{a=1,b}^m gcd(a,m)|b\\&=

2017-09-29 12:37:06 343

原创 HDU 5517 2015沈阳现场赛(二维树状数组)

题目链接Triple题目分析首先分析集合C中的不同点不会超过10510^5个,因为每个ee 只会对应唯一的一个aa,因为我们只需要对每个bb存一个最大的a就行. 然后将CC中的点集排序,重大到小排序，重大到小遍历，用二维树状数组判断　≥(c,d)\ge (c,d)　的点的个数,若多于１个则说明这个点不符合情况,否则符合条件,因为在二维树状数组中的点的a都比当前点大,所以一次少过去就好.注意计算重点

2017-09-25 22:55:29 307

原创 2017北京网络赛hihocoder 1580 matrix(dp)

题目链接1580 : Matrix分析和最大子矩阵和类似，扩展一维就好了. 详细分析请见这位老歌写的blog 注意特判Ac 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

2017-09-25 21:29:29 327

原创 HDU The Closest M Points（KD-tree模板题）

题目链接HDU 4347 The Closest M Points分析kd-tree模板题AC 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#defi

2017-09-25 16:13:06 651

原创数论经典问题:与x互质的数的和

问题xx已知，求 ∑gcd(a,x)=1a\sum_{gcd(a,x)=1}a 我们知道这样的　aa 只有ϕ(x)\phi(x) 个，我们想能否用经典的gauss求等差数列的方法来解决这个问题呢? s设　gcd(a,x)=1gcd(a,x)=1，那麽有　gcd(x−a,x)=1gcd(x-a,x)=1因此对于任意的数　a,gcd(a,x)=1a,gcd(a,x)=1 则有　对应的x−a，满

2017-09-23 11:42:09 1312

原创 2015沈阳现场赛F (HDU 5514)(经典问题数论phi函数)

题目链接Frogs分析首先我们可以发现能跳到的数一定满足存在i,gcd(ai,m)|x i,gcd(a_i,m) | x,即存在　ii gcd(ai,m)|gcd(x,m)gcd(a_i,m) | gcd(x,m) 将 x,按照　gcd(x,m)=dx,按照　gcd(x,m)=d 分类，最多有　mm的因子个集合,那麽每个集合的和为 ∑gcd(x,m)=dx=d∗∑gcd(x/d,m/d)

2017-09-23 11:33:51 300

原创 2017数学建模总结

2017年数学建模第一天晚上选题，看了看A，看了看B，B题分析性太强，属于应用统计，数据分析类，和队友商量着战略性放弃，然后就开始查资料，搞东西了，突然回忆起以前姜起源的书上貌似讲过CT成像这种东西,翻出来看果然有这玩意儿。然后看了看题目发现精度是4位小数，似乎可以暴力枚举答案

2017-09-19 01:47:30 8111 1

原创 python爬虫自动提交HDU并获取AC状态(p3+request+Beatifulsoup)

说明今天下午搞了好一会儿，源于上周日的沈阳网络赛,竟然可以枚举答案交1024发……想到如果自己能写个爬虫自动提交那该多好啊,源于这总冲动自己今天花了一个下午实现了这个小玩意儿.效果图实现过程其实主要是3个步骤:login 登录submit 提交代码report status 范围状态(AC没有?)首先是登录，就是要把你的id 和　密码post 到HDU ，我们可以打开HDU 的主页然后分析提

2017-09-13 22:36:24 897

原创 HDU 4109 Instrction Arrangement(差分约束)

题目链接Instrction Arrangement分析这是一道简单的差分约束的题目，关于差分约束，这个blog写的很好设tit_i 表示第i条指令完成的最短时间，那麽对于一个约束x,y,zx,y,z有ty−txtx−ty≥z≤z\begin{align}t_y -t_x&\ge z\\t_x - t_y &\le z\end{align} 这不就是差分约束的标准形式吗,要求完成的最

2017-09-13 21:56:28 345

原创 poj 3317 Stake Your Claim(极大极小搜索经典 a-b+剪枝+记忆化dp)

poj 3317 Stake Your Claim分析哇这题也算把我卡得够死了，怪我没能正确理解所谓a-b剪枝的姿势. 其实alpha-beta剪枝不用也是可以过的，非得我不用才找到了错误首先这道题只有10种选择，每种选择有3种状态，共310<600003^{10} <60000,因此直接记忆话搜就行了.用3进制状压一下.然后就是枚举空点的时候也可以状压一下这样枚举快.AC code//Source Cod

2017-09-12 00:34:53 1964

原创 kd-tree : k近邻查询和范围查询

想象一下我们有如下两个任务:我现在想骑一辆小黄车，我想查找离我最近的k辆小黄车.找到百度地图中显示在屏幕上区域中的所有酒店这两个任务均可以用kd-tree来解决 kd-tree 主要两个用途:查询离某个点的最近的kk 个邻居,搜索某个区域内的所有点.后者在计算几何中称为范围查询,例如查询某个平面区域内的点的个数.kd-tree是什么玩意儿kd-tree就是高维平衡树…… kd-tre

2017-09-11 16:07:13 11702 4

原创 STL nth_element神器

这周在看kd-tree,苦于找中位数写起来太过麻烦，去观摩了一下其他大佬的代码，竟然还有nth_element这种东西。 http://www.cplusplus.com/reference/algorithm/nth_element/ nth_element(first,nth,last) first，last 第一个和最后一个迭代器，也可以直接用数组的位置。 nth，要定位的第n

2017-09-08 13:52:19 6845

原创 ACM概率期望dp刷题总结

这个周刷了很多概率期望有关的dp题目，缘起2016青岛D题和取log的神操作题HDU 5988 2016青岛区域赛 (最小费用流)这类题目没有固定的模板，而且概率可以很容易插入一些经典模型，比如下面的 TSP。最小费用流。我刷的这部分题目，无一例外均可用dp解决，主要是找到状态，很多题目都可以抽象成马尔科夫链下面记录一下刷的经典题目B - Discovering Gold 一排1到n的格子，每个

2017-09-07 22:10:48 3628

原创 codeforces 854 D. Jury Meeting(dp水题)

题目链接D. Jury Meeting分析昨晚打正式赛的时候太急了，代码写的有点毒，下了课后一看，立马发现bug…..将去的航班表和返回航班表分别按照时间最小和最大排序，对于去的航班来说，挨个扫描过去，对每一个人维护当前的最小花费，若时间总的航班满足等于　nn 说明这一天可以作为起始时间，更新这一天的花费dp[d]=now_ansdp[d] = now\_ans,now_ansnow\_ans 记录

2017-09-07 12:10:53 586

原创线性规划:单纯形算法

作者 dylanFrank(滔滔)转载请联系作者原文链接：http://blog.csdn.net/Dylan_Frank/article/details/77876006这里简要总结一下线性规划的单纯形算法，做如下几个方面的总结，其余以后再来填坑.几何表示标准型的代数解法其他情形与标准型的互换退化情形几何表示先看这样一个问题我们很容易用下面的数学语言来

2017-09-07 00:18:58 9667 1

原创 HDU 5988 2016青岛区域赛 (最小费用流)

题目链接HDU 5988分析这题一看就是一个网络流的模板题，不过需要注意建边的费用!! 首先很容易想到，s与每个区域连一条边，费用为0，容量为 sis_i,区域与 tt 连一条边容量为 bib_i 费用为0。然后就是考虑区域与区域之间的连边了。我们想象一下最终的分配方案，假设对于边 eije_{ij} 来说，我们让 kijk_{ij} 个人通过，那么eije_{ij} 不被破坏的概率是 (1−p

2017-09-06 16:46:46 2303

原创 codeforces 851 D. Arpa and a list of numbers(前缀和+bruteforce)

题目链接D. Arpa and a list of numbers分析是我太天真了，暴力枚举素数还是计算还是有问题的，复杂度太高，不过你可以发现只需要再处理一个前缀和就可以了….下面进入正题首先显然可以想到枚举素数，那么最终就是对于每个素数的计算问题了，如果暴力计算的化肯定是要超时的，其实我们可以在O(1)O(1) 计算花费.对于素数 kk 来说，不为kk 的倍数的数都要加到kk 的倍数,所以只需要

2017-09-05 09:11:21 772

taotao 的大学墓志