ZOJ3874 Permutation Graph(NTT&&cdq分治)
最近在看几道整体二分还有cdq分治的东西,突然间想起前几个礼拜的ZOJ题,然后看了一下代码,经过了一些深思熟虑之后,发现自己终于看懂了,下面就用别人的代码来剖析一下整个解题的思路吧,具体的内容我再看看。
首先要解决这个问题需要有一些卷积的知识,或者说是多项式乘法,一个很典型的多项式乘法的东东就是FFT,然后原来在数论意义下(即mod P)的意义下,也有相应的NTT(快速数论变换),思想是和FFT一致的,不过在这里原根稍微不一样,而且也不用去管复数,当然我是不会懂的,下面抄个链接:
NTT(快速数论变换) http://blog.csdn.net/zz_1215/article/details/40430041
现在假设自己懂了NTT,然后我们把它当作黑盒,然后就把它当作可以实现mod P意义下卷积的一个工具,然后去理解一下题目的做法。
经过一些理论推导,我们可以发现,实际上我们要求的东西是 dp[n]=n!-(dp[1]*(n-1)!+dp[2]*(n-2)!+... dp[n-1]*1!).
实际上dp[n]=n!- dp[i]和i!的卷积的第n项。这样的一个算法暴力做的话要算到n的话是O(n^2)的,下面看下cdq分治。
个人对cdq分治的理解是这样的:
T(n)=2T(n/2)+O(f(n)) 一个传统的典型的分治算法里,O(f(n))是指的将两个子问题合并的代价,非常典型的就是归并排序。而在cdq分治里,O(f(n))就不一定是合并的代价了,在归并排序里,左子问题对右子问题是没有影响的,而现实的分治里,可能会出现左子问题对右子问题产生影响的情况,及前面的操作是直接对右边操作产生影响的,我们必须先做了左子问题,然后把左子问题的影响加到右子问题,然后才能再递归右子问题。下面的框架给了两种分治的思想吧。
// traditional divide-and-conquer
void solve(l,r)
{
int mid=(l+r)>>1;
solve(l,mid);
solve(mid+1,r);
combine 2 subproblem.
} // cdq divide-and conquer void solve1(l,r)
{
int mid=(l+r)>>1;
solve(l,mid);
add the affect of [l,mid] to [mid+1,r]
solve(mid+1,r);
} void solve2(l,r)
{
int mid=(l+r)>>1;
add the affect of [l,mid] to [mid+1,r];
solve(l,mid);
solve(mid+1,r);
}
下面我们来看看别人的代码里是怎么做的,下面的代码抄自下面的链接,注释是自己给别人加上的,方便一下理解吧。
ZOJ3874:http://acm.hust.edu.cn/vjudge/problem/viewSource.action?id=3709562
// dp[l]要求的
// f[l] 是l的阶乘
void solve (int l, int r) {
// 递归边界,l==r时,说明所有比l小的卷积都算完了,所以dp[l]=f[l]-dp[l];
if (l == r) {
dp[l] = (f[l] - dp[l] + P) % P;
return;
}
int m = (l + r) >> 1;
// 递归左子问题,现在要做的是算出dp[l...mid]和f[]的卷积加到dp[mid+1...r]上
solve (l, m);
// 下面的部分就是将dp[l...mid]赋给a,将f的值赋给b,然后做NTT,然后算完之后再逆变换回来
// 做完逆变换后,a[x]存的就是dp和f的卷积的第x项,
int s = 1, n = m - l + 1;
while (s <= n * 2) s <<= 1;
a[0] = b[0] = 0;
for (int i = 1, j = l; i < s; i++, j++) a[i] = (j <= m ? dp[j] : 0);
for (int i = 1; i < s; i++) b[i] = f[i];
NTT (a, s); NTT (b, s);
for (int i = 0; i < s; i++) a[i] = a[i] * b[i] % P;
NTT (a, s, true);
// end
// 将影响加到dp[mid+1...r]上。
for (int i = m + 1, j = m - l + 2; i <= r; i++, j++) (dp[i] += a[j]) %= P;
solve (m + 1, r);
}
最后不难发现将左子问题的影响加到右子问题上其实就是一个O(nlogn)的过程。
T(n)=2T(n/2)+O(nlogn)
所以最后出来的复杂度应该是O(nlog^2n)的。
ZOJ3874 Permutation Graph(NTT&&cdq分治)的更多相关文章
- ZOJ3874 Permutation Graph
Time Limit: 2 Seconds Memory Limit: 65536 KB Edward has a permutation {a1, a2, … an}. He finds ...
- ZOJ3874 Permutation Graph 【分治NTT】
题目链接 ZOJ3874 题意简述: 在一个序列中,两点间如果有边,当且仅当两点为逆序对 给定一个序列的联通情况,求方案数对\(786433\)取模 题解 自己弄了一个晚上终于弄出来了 首先\(yy\ ...
- ZOJ 3874 Permutation Graph 分治NTT
Permutation Graph Time Limit: 2 Seconds Memory Limit: 65536 KB Edward has a permutation {a1, a2 ...
- HDU5322 Hope(DP + CDQ分治 + NTT)
题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5322 Description Hope is a good thing, which can ...
- hdu_5354_Bipartite Graph(cdq分治+并查集判二分图)
题目链接:hdu_5354_Bipartite Graph 题意: 给你一个由无向边连接的图,问对于每一个点来说,如果删除这个点,剩下的点能不能构成一个二分图. 题解: 如果每次排除一个点然后去DFS ...
- 【BZOJ-3456】城市规划 CDQ分治 + NTT
题目链接 http://www.lydsy.com/JudgeOnline/problem.php?id=3456 Solution 这个问题可以考虑dp,利用补集思想 N个点的简单图总数量为$2^{ ...
- BZOJ4555求和(cdq分治+NTT)
题意: 输出f(n)对998244353(7 × 17 × 223 + 1)取模的结果.1 ≤ n ≤ 100000 其中S(i,j)是第二类Stirling数,即有i个球,丢到j个盒子中,要求盒子不 ...
- BNUOJ 51279[组队活动 Large](cdq分治+FFT)
传送门 大意:ACM校队一共有n名队员,从1到n标号,现在n名队员要组成若干支队伍,每支队伍至多有m名队员,求一共有多少种不同的组队方案.两个组队方案被视为不同的,当且仅当存在至少一名队员在两种方案中 ...
- HDU5730 FFT+CDQ分治
题意:dp[n] = ∑ ( dp[n-i]*a[i] )+a[n], ( 1 <= i < n) cdq分治. 计算出dp[l ~ mid]后,dp[l ~ mid]与a[1 ~ r-l ...
随机推荐
- Java语言基础---转义符
转义符 转义符使用“\”表示.常用转义符如下: 1.‘\n’回车 2.‘\t’制表位字符,一个表示向右跳8-10个字符 3.‘\\’表示’\’ 4.‘\’’表示单引号 5.‘\’’’表示双引号 6.‘ ...
- poj2001Shortest Prefixes(trie)
Shortest Prefixes Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 18687 Accepted: 808 ...
- Spring mvc+hibernate+freemarker(实战)
Spring mvc+hibernate+freemarker(实战) 博客分类: Spring Spring mvchibernatefreemarkerwebjava 今天我为大家做了一个 sp ...
- java 解析/读取 种子/bt/torrent 内容
碰到不会的技术问题,我还是先度娘.能中文看懂,为什么非要看英文呢. java 解析/读取 种子/bt/torrent 内容,这个度娘给的满意答案并不是很多.GG之后的搜索结果出现了stackover ...
- 关于IOS下日期格式分隔符 - 、 /的问题
之前我们项目有一个低价日历,服务端下发的时间格式为: "2014-07-21 09:45:12" 然后一直出不了数据,后来发现. IOS下无论chrome.safari还是Uc如 ...
- HTML简易学习笔记
文字版地址 https://github.com/songzhenhua/github/blob/master/HTML简易学习笔记.txt
- 爬虫:Scrapy15 - 调试(Debugging)Spiders
考虑下面的 spider: import scrapy from myproject.items import MyItem class MySpider(scrapy.Spider): name = ...
- HDU 4109 Instrction Arrangement(DAG上的最长路)
把点编号改成1-N,加一点0,从0点到之前任意入度为0的点之间连一条边权为0的边,求0点到所有点的最长路. SPFA模板留底用 #include <cstdio> #include < ...
- Spring 学习笔记(六)—— AOP的简单理解
系统中的业务可以分为核心关注点和横切关注点. 核心关注点时业务处理的主要流程,而横切关注点是与核心业务无关但更为通用的业务. 各个横切关注点离散地穿插于核心业务之中,导致系统地每一个模块都与这些业务具 ...
- 获取任意网站的图标,标题栏logo,网站logo
https://www.hao123.com/favicon.ico 网站换成你想要的 大多数都可以