2017年05月_Nicetomeetu-

10月 09月 08月 07月 06月 05月 04月 03月 02月 01月

原创 HDU 2089 不要62 初探数位dp

第一次写传说中的数位dp，久仰大名终于得见真容，无比激动！数位dp就是求在给定的区间[l, r]内满足条件C的数字的个数，这类区间统计问题往往可以用数学上的递推来描述，也就是dp了。关于这道题，定义状态dp[i][j]表示以数字j开头的i位数字中不含62和4的个数。什么意思呢？dp[2][6]就表示以数字6开头的两位数字中不包含62和4的个数。也就是61,63,65,66,67,68,69

2017-05-27 19:31:53 606

原创 NJUST 1722 所有的平方差 Pollard_Rho大数分解

原题链接：所有的平方差代码如下：#include #include #include #include #include #include #include using namespace std;typedef long long int ll;const int MAX_SIZE = 3000;const ll mod = 10000000019;ll

2017-05-25 21:56:43 676

原创 HDU 4344 Mark the Rope Pollard_Rho大数分解

PS：图片由作者用wps制作，使用请注明链接，O(∩_∩)O谢谢！注意如果n的质因数分解只用到一个质数时，要多除一个质数。如对Miller-Rabin有疑问，请参考作者其他博客：Miller-Rabin算法详解代码如下：#include #include #include #include #include #include using namespace s

2017-05-24 19:03:36 490

原创 POJ 2429 GCD & LCM Inverse Pollard_Rho大数分解+Miller_Rabin

PS：图片由作者用wps制作，使用请注明链接，O(∩_∩)O谢谢！代码如下：#include #include #include #include #include #include using namespace std;typedef long long int ll;const int MAX_SIZE = 1000;const ll INF = 0x

2017-05-24 17:46:45 556

原创 Miller-Rabin素性测试算法详解

看了一些别人的博客，发现里面涉及到的公式没有证明，于是就打算自己写一篇比较详细的讲解。先看两个引理及其证明(建议把证明搞懂)。PS：以下图片均为作者用wps制作，如想使用请附上作者博客链接，谢谢O(∩_∩)O。看完了上面的引理，那就可以正式开始Miller-Rabin算法的讲解了。背景：素性测试(即测试给定的数是否为素数)是近代密码学中的一个非常重要的课题。虽然Wilson

2017-05-24 14:54:09 26521 7

原创 POJ 1930 Dead Fraction 数论

别人说这是小学奥数水平的题目，orz，作者受到了极大的打击。如果是用笔写的话相信很多人都能将循环小数化成分数。来个例子吧：0.333.....怎么化成分数呢？设 x = 0.333......10x = 3.33......10x - x = 39x = 3x = 1 / 3和这个思路类似，对于一般情况，我们可以推导出一个公式(推导出了公式才方便编程)。引用别人

2017-05-21 15:48:50 466

原创 POJ 1679 The Unique MST(算法导论23-1次优最小生成树)

只需要讲解算法导论的题即可。23-1次优最小生成树a. 最小生成树唯一性证明：已知当前构造的边集A是最小生成树的子集。令无向图G的一个切割是，显然该切割是尊重A的。已知跨越该切割的轻量级边对于A是安全的，又因为该无向图G的每条边的权值都不相同，所以对于当前A而言，安全边有且只有一条，即对于每个状态下的A，构造最小生成树的方式是唯一的。所以最小生成树是唯一的。次优最小

2017-05-13 00:25:09 2056 1

原创 HDU 1495 非常可乐 bfs状态空间搜索

状态空间搜索就是将问题求解过程表现为从初始状态到目标状态寻找路径的过程。普通的搜索是在显式图中进行，从一个节点搜索到其他节点。状态空间搜索是在一个由问题状态构成的隐式图中进行，从一个状态搜索到另一个状态。状态空间搜索和普通搜索的区别无非是，这里搜索的不再是一个个实实在在的节点，而是一个个问题的状态，不过二者思路和方法都是相同的。本题问题的状态可以用四个值来描述，第一个杯子的

2017-05-06 21:19:48 461

原创 POJ 3356 AGTC（算法导论15-5编辑距离）经典dp

POJ 3356 为算法导论原题简化版，这里只讲算法导论版编辑距离。定义状态dp[i][j]表示将x[i……m]变为y[j……n]所需要的最小代价。对于每一个状态，有6种操作，具体如下。1、复制功能：i + 1，j + 1代价：cost12、替换功能：i + 1，j + 1代价：cost23、删除功能：i + 1，j不变代价：cost34

2017-05-05 16:11:44 788

原创 POJ 2677 （算法导论15-3）双调欧几里得旅行商问题 dp

首先按横坐标递增给所有点排序。定义状态dp[i][j]表示从点i向n走一条路L1，从点j向n走另一条路L2（如下图，两条路互不相交，并且L1在L2上面），L1 + L2的最小值。程序中用distance(i, j)表示点i到点j的距离。如何计算dp[i, j]呢?我们考虑k = max(i, j) + 1这个点，这个点肯定在L1或者L2上。k在L1上时，如下图

2017-05-05 16:00:40 930 1

原创 UVA 1103 Ancient Messages dfs

按照《算法竞赛入门经典》的思路，我们只需要判断出当前文字有几个“洞”，就可以确定到底是哪个文字了。W：零个洞A：一个洞K：两个洞J：三个洞S：四个洞D：五个洞本题重点就是如何判断一个文字有几个洞。我们知道每一个文字都是由很多个 "1" 组成的，背景是 “0”。请看下面这个文字：0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

2017-05-01 17:12:41 534