【BZOJ3456】城市规划

最新推荐文章于 2019-07-07 20:03:40 发布
最新推荐文章于 2019-07-07 20:03:40 发布
阅读量414 收藏
点赞数
分类专栏： 【OJ】BZOJ 【类型】做题记录 【算法】FFT与NTT 【算法】CDQ分治与整体二分 【算法】生成函数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39972971/article/details/80364039
版权

【题目链接】

【思路要点】

  • 令\(f_i\)为当\(N=i\)时，原问题的答案。
  • 显然有\(f_1=1\)。
  • 当\(i>1\)，考虑用\(i\)个点的简单无向图的个数减去\(i\)个点的简单无向图中不连通的个数。
  • 考虑经典思路：枚举1号节点所在连通块的大小\(j\)，则有\(f_i=2^{\binom{i}{2}}-\sum_{j=1}^{i-1}f_j*\binom{i-1}{j-1}*2^{\binom{i-j}{2}}\)。
  • 展开组合数，整理得\(f_i=2^{\binom{i}{2}}-(i-1)!*\sum_{j=1}^{i-1}\frac{f_j}{(j-1)!}*\frac{2^{\binom{i-j}{2}}}{(i-j)!}\)。
  • 直接使用分治NTT求解，时间复杂度\(O(NLog^2N)\)，常数较好的实现可以通过本题。
  • 观察上式，若将\(f_i\)移项，则有：\(2^{\binom{i}{2}}-(i-1)!*\sum_{j=1}^{i}\frac{f_j}{(j-1)!}*\frac{2^{\binom{i-j}{2}}}{(i-j)!}=0\)。
  • 也即\(\frac{2^{\binom{i}{2}}}{(i-1)!}=\sum_{j=1}^{i}\frac{f_j}{(j-1)!}*\frac{2^{\binom{i-j}{2}}}{(i-j)!}\)。
  • 令\(C_i=\frac{2^{\binom{i}{2}}}{(i-1)!}\)，\(A_i=\frac{f_i}{(i-1)!}\)，\(B_i=\frac{2^{\binom{i}{2}}}{i!}\)，则上式等价于\(C=A*B\)。
  • 其中\(C\)和\(B\)均为已知多项式，通过对\(B\)求逆倒推用\(A=C*B^{-1}\)求得\(A\)即可。
  • 时间复杂度\(O(NLogN)\)，可以通过本题。

【代码】

/*Inverse of Polynomial O(NLogN)*/
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 524288;
const int P = 1004535809;
template <typename T> void chkmax(T &x, T y) {x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {x = min(x, y); } 
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
namespace NTT {
	const int G = 3;
	int power(int x, int y) {
		if (y == 0) return 1;
		int tmp = power(x, y / 2);
		if (y % 2 == 0) return 1ll * tmp * tmp % P;
		else return 1ll * tmp * tmp % P * x % P;
	}
	int powerll(int x, long long y) {
		if (y == 0) return 1;
		int tmp = powerll(x, y / 2);
		if (y % 2 == 0) return 1ll * tmp * tmp % P;
		else return 1ll * tmp * tmp % P * x % P;
	}
	int N, Log, home[MAXN];
	void NTTinit() {
		for (int i = 0; i < N; i++) {
			int ans = 0, tmp = i;
			for (int j = 1; j <= Log; j++) {
				ans <<= 1;
				ans += tmp & 1;
				tmp >>= 1;
			}
			home[i] = ans;
		}
	}
	void NTT(int *a, int mode) {
		for (int i = 0; i < N; i++)
			if (home[i] < i) swap(a[i], a[home[i]]);
		for (int len = 2; len <= N; len <<= 1) {
			int delta;
			if (mode == 1) delta = power(G, (P - 1) / len);
			else delta = power(G, P - 1 - (P - 1) / len);
			for (int i = 0; i < N; i += len) {
				int now = 1;
				for (int j = i, k = i + len / 2; k < i + len; j++, k++) {
					int tmp = a[j];
					int tnp = 1ll * a[k] * now % P;
					a[j] = (tmp + tnp) % P;
					a[k] = (tmp - tnp + P) % P;
					now = 1ll * now * delta % P;
				}
			}
		}
		if (mode == -1) {
			int inv = power(N, P - 2);
			for (int i = 0; i < N; i++)
				a[i] = 1ll * a[i] * inv % P;
		}
	}
	void times(int *a, int *b, int *c, int limit) {
		N = 1, Log = 0;
		while (N < 2 * limit) {
			N <<= 1;
			Log++;
		}
		for (int i = limit; i < N; i++)
			a[i] = b[i] = 0;
		NTTinit();
		NTT(a, 1);
		NTT(b, 1);
		for (int i = 0; i < N; i++)
			c[i] = 1ll * a[i] * b[i] % P;
		NTT(c, -1);
	}
	void timesabb(int *a, int *b, int *c, int limit) {
		N = 1, Log = 0;
		while (N < 2 * limit) {
			N <<= 1;
			Log++;
		}
		for (int i = limit; i < N; i++)
			a[i] = 0;
		for (int i = limit / 2; i <= N; i++)
			b[i] = 0;
		NTTinit();
		NTT(a, 1);
		NTT(b, 1);
		for (int i = 0; i < N; i++)
			c[i] = 1ll * a[i] * b[i] % P * b[i] % P;
		NTT(c, -1);
	}
	void inverse(int *a, int *b, int limit) {
		b[0] = power(a[0], P - 2);
		for (int now = 1; now < limit; now <<= 1) {
			static int c[MAXN], d[MAXN];
			for (int i = 0; i < now * 2; i++)
				c[i] = a[i], d[i] = b[i];
			timesabb(c, d, d, now * 2);
			for (int i = 0; i < now * 2; i++)
				b[i] = (2ll * b[i] - d[i] + P) % P;
		}
	}
}
int n, fac[MAXN], inv[MAXN], two[MAXN];
int invb[MAXN], b[MAXN], c[MAXN], a[MAXN];
int main() {
	read(n);
	fac[0] = 1;
	for (int i = 1; i <= n; i++)
		fac[i] = 1ll * fac[i - 1] * i % P;
	inv[n] = NTT::power(fac[n], P - 2);
	for (int i = n - 1; i >= 0; i--)
		inv[i] = 1ll * inv[i + 1] * (i + 1) % P;
	for (int i = 1; i <= n; i++)
		two[i] = NTT::powerll(2, i * (i - 1ll) / 2);
	b[0] = 1; c[0] = 0;
	for (int i = 1; i <= n; i++) {
		b[i] = 1ll * two[i] * inv[i] % P;
		c[i] = 1ll * two[i] * inv[i - 1] % P;
	}
	NTT::inverse(b, invb, n + 1);
	NTT::times(invb, c, a, n + 1);
	writeln(1ll * a[n] * fac[n - 1] % P);
	return 0;
}
/*DivideAndConquer + NTT O(NLog^2N)*/
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 262144;
const int P = 1004535809;
template <typename T> void chkmax(T &x, T y) {x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {x = min(x, y); } 
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
namespace NTT {
	const int G = 3;
	int power(int x, int y) {
		if (y == 0) return 1;
		int tmp = power(x, y / 2);
		if (y % 2 == 0) return 1ll * tmp * tmp % P;
		else return 1ll * tmp * tmp % P * x % P;
	}
	int powerll(int x, long long y) {
		if (y == 0) return 1;
		int tmp = powerll(x, y / 2);
		if (y % 2 == 0) return 1ll * tmp * tmp % P;
		else return 1ll * tmp * tmp % P * x % P;
	}
	int N, Log, home[MAXN];
	void NTTinit() {
		for (int i = 0; i < N; i++) {
			int ans = 0, tmp = i;
			for (int j = 1; j <= Log; j++) {
				ans <<= 1;
				ans += tmp & 1;
				tmp >>= 1;
			}
			home[i] = ans;
		}
	}
	void NTT(int *a, int mode) {
		for (int i = 0; i < N; i++)
			if (home[i] < i) swap(a[i], a[home[i]]);
		for (int len = 2; len <= N; len <<= 1) {
			int delta;
			if (mode == 1) delta = power(G, (P - 1) / len);
			else delta = power(G, P - 1 - (P - 1) / len);
			for (int i = 0; i < N; i += len) {
				int now = 1;
				for (int j = i, k = i + len / 2; k < i + len; j++, k++) {
					int tmp = a[j];
					int tnp = 1ll * a[k] * now % P;
					a[j] = (tmp + tnp) % P;
					a[k] = (tmp - tnp + P) % P;
					now = 1ll * now * delta % P;
				}
			}
		}
		if (mode == -1) {
			int inv = power(N, P - 2);
			for (int i = 0; i < N; i++)
				a[i] = 1ll * a[i] * inv % P;
		}
	}
	void times(int *a, int *b, int *c, int limit) {
		N = 1, Log = 0;
		while (N < 2 * limit) {
			N <<= 1;
			Log++;
		}
		for (int i = limit; i < N; i++)
			a[i] = b[i] = 0;
		NTTinit();
		NTT(a, 1);
		NTT(b, 1);
		for (int i = 0; i < N; i++)
			c[i] = 1ll * a[i] * b[i] % P;
		NTT(c, -1);
	}
}
int n, f[MAXN], fac[MAXN], inv[MAXN], two[MAXN];
void work(int l, int r) {
	if (l == r) return;
	int mid = (l + r) / 2;
	work(l, mid);
	int limit = r - l + 1;
	static int a[MAXN], b[MAXN], c[MAXN];
	for (int i = l; i <= mid; i++)
		a[i - l] = 1ll * f[i] * inv[i - 1] % P;
	for (int i = mid + 1; i <= r; i++)
		a[i - l] = 0;
	b[0] = 0;
	for (int i = 1; i < limit; i++)
		b[i] = 1ll * two[i] * inv[i] % P;
	NTT::times(a, b, c, limit);
	for (int i = mid + 1; i <= r; i++)
		f[i] = (f[i] - 1ll * fac[i - 1] * c[i - l] % P + P) % P;
	work(mid + 1, r);
}
int main() {
	read(n);
	fac[0] = 1;
	for (int i = 1; i <= n; i++)
		fac[i] = 1ll * fac[i - 1] * i % P;
	inv[n] = NTT::power(fac[n], P - 2);
	for (int i = n - 1; i >= 0; i--)
		inv[i] = 1ll * inv[i + 1] * (i + 1) % P;
	for (int i = 1; i <= n; i++)
		f[i] = two[i] = NTT::powerll(2, i * (i - 1ll) / 2);
	work(1, n);
	writeln(f[n]);
	return 0;
}
cz_xuyixuan
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
bzoj 3456 城市规划
我要吃熊猫
06-29 439
多项式求逆，NTT
BZOJ 3456 城市规划
Fuxey
05-11 1489
FNT
bzoj3456 城市规划
wennitao的博客
06-10 159
题面：bzoj3456 题意：求nnn个点的简单无向连通图的方案数 题解：设g(n)g(n)g(n)为简单无向图的方案数，f(n)f(n)f(n)为简单无向连通图的方案数 g(n)=2Cn2g(n) = 2^{C_n^2}g(n)=2Cn2​ 枚举111号点的连通块大小：g(n)=∑i=1nCn−1i−1f(i)g(n−i)g(n) = \sum _{i=1} ^n C_{n-1} ^{i-1} ...
BZOJ3456 城市规划
aoye9670的博客
11-30 73
一句话题意 求\(n\)个带标号点能构成的简单无向图数目 考虑令\(f[i]\)表示\(i\)个带标号点的简单无向图数目. \(f[n]=2^{C_n^2}-\sum_{i=1}^{n-1}C_{n-1}^{i-1}*f[i]*2^{C_{n-i}^2}\) 两边同除\((n-1)!\)得 \(\frac{f[n]}{(n-1)!}=\frac{2^{C_n^2}}{(n-1)!}-\...
Bzoj3456 城市规划
dezhen7015的博客
05-01 103
Time Limit: 40 SecMemory Limit: 256 MBSubmit: 681Solved: 383 Description 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通.为了省钱, 每两个城市之间最...
BZOJ3456 城市规划（多项式求逆）
weixin_30839881的博客
01-18 71
　　设f[i]为连通图的数量，g[i]为不连通图的数量，显然有f[i]=2i*(i-1)/2-g[i]，g[i]通过枚举1所在连通块大小转移，有g[i]=Σf[j]*C(i-1,j-1)·2(i-j)*(i-j-1)/2，也即f[i]=2i*(i-1)/2-(i-1)!·Σf[j]·2(i-j)*(i-j-1)/2/(j-1)!/(i-j)!。显然是一个卷积形式，可以分治NTT。 　　进一步将式...
BZOJ3456城市规划
weixin_30776273的博客
02-15 57
题目描述 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了.刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通.为了省钱, 每两个城市之间最多只能有一条直接的贸易路径. 对于两个建立路线的方案, 如果存在一个城市对, 在两个方案中是否建立路线不一样, 那么这两个方案就是不同的, 否则就是相同的. 现在你需要求出一...
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
BZOJ题目镜像
09-14
通过这个镜像，你可以查看到BZOJ自2012年02月29日以来的所有题目，涵盖了各种难度和类型的算法问题，如动态规划、图论、排序算法、搜索算法等。这对于编程爱好者和准备参加ACM/ICPC（国际大学生程序设计竞赛）或其他...
bzoj problem
08-21
这些题目通常涵盖从基础的逻辑思维到复杂的算法设计，包括但不限于排序、搜索、图论、动态规划、数学逻辑等多个领域。用户可以通过解题来锻炼自己的编程思维和问题解决能力。 【标签】"bzoj" 是这个平台的标识，...
BZOJ第二部分
11-04
【BZOJ第二部分】是针对BZOJ（北京邮电大学在线评测系统）的一个专题，这个专题主要涵盖了BZOJ平台上的P2001到P3000之间的编程题目。对于想要深入学习算法、提升编程能力的IT从业者或学生来说，这是一个宝贵的资源...
【LOJ2320】「清华集训 2017」生成树计数
cz_xuyixuan的博客
06-24 1814
【题目链接】点击打开链接【思路要点】由于场。时间复杂度\(O(DC^2+D^2+TD^2LogN)\)，其中\(C=165,D=146\)。【代码】#include&lt;bits/stdc++.h&gt; using namespace std; const int MAXM = 65; const int MAXN = 262144; const int MAXLOG = 20; const ...
【2019 江苏省队集训】Day1 解题报告
cz_xuyixuan的博客
07-02 1343
【T1】 光影交错 【思路要点】 NNN 轮后仪式仍然进行的概率为 (1−p)N(1-p)^N(1−p)N ，当 p=10−5p=10^{-5}p=10−5 时，取 N&gt;107N&gt;10^7N>107 可以保证该概率在 10−2010^{-20}10−20 以内，因此忽略这部分情况不会导致精度要求不能接受的误差产生。 记 f(i)f(i)f(i) 表示出现 ii...
【LOJ575】「LibreOJ NOI Round #2」不等关系
cz_xuyixuan的博客
07-07 1283
【题目链接】 点击打开链接 【思路要点】 见 官方题解 。（反正也是我写的） 时间复杂度 O(NLog2N)O(NLog^2N)O(NLog2N) 。 【代码】 #include<bits/stdc++.h> using namespace std; const int MAXN = 262144; const int P = 998244353; typedef l...
【CodeForces】Mail.Ru Cup 2018 Round 1 (Div. 1 + Div. 2) 题解
cz_xuyixuan的博客
10-21 774
【比赛链接】 点击打开连接 【题解链接】 点击打开链接 **【A】**Elevator or Stairs? 【思路要点】 按照题意计算两种方式的用时，取较优的方案采纳。 时间复杂度 O(1)O(1)O(1) 。 【代码】 #include&lt;bits/stdc++.h&gt; using namespace std; const int MAXN = 2e5 + 5...
【LOJ2554】「CTSC2018」青蕈领主
cz_xuyixuan的博客
04-19 739
【题目链接】 点击打开链接 【思路要点】 值域连续的区间的交和并均为值域连续的区间，因此一个合法输入应当能够构成一个树形结构。 考虑一个节点的子树，其对应的区间值域连续，且在后续的考虑中，我们会避免使其出现在更大的连续区间中，因此我们可以等效地将其看做一个数。 因此，记 fif_ifi​ 表示长度为 i+1i+1i+1 的，去掉最后一个元素后不存在非平凡连续区间的排列个数， sis_i...
【LOJ2983】「WC2019」数树
cz_xuyixuan的博客
03-31 623
【题目链接】 点击打开链接 【思路要点】 op=0op=0op=0 ，算一算两棵树的公共边数即可。 时间复杂度 O(N)O(N)O(N) 或 O(NLogN)O(NLogN)O(NLogN) 。 op=1op=1op=1 ，考虑枚举蓝树上的一个边集 SSS ，强制红树上同样存在这些边，计算将剩余 N−∣S∣N-|S|N−∣S∣ 个联通块连成一棵树的方案数，更新答案。 上述算法中，枚举蓝...
【CodeForces】CodeForces Round #485 (Div. 1 + Div. 2) 题解
cz_xuyixuan的博客
06-05 561
【比赛链接】Div. 1Div. 2【题解链接】点击打开链接【Div.2 A】Infinity Gauntlet【思路要点】按照题意模拟即可。时间复杂度\(O(NLogN)\)。【代码】#include&lt;bits/stdc++.h&gt; using namespace std; const int MAXN = 100005; template &lt;typename T&gt; voi...
bzoj.org俄罗斯方块
最新发布
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。在该游戏中，玩家需要控制方块不断向下移动，直到方块落到底部或者与其他方块重叠。游戏过程中，玩家可以通过旋转和移动方块来调整方块的位置，使其能够合理地落下。当一行方块被填满时，该行会被清除，并获得相应的分数。 bzoj.org俄罗斯方块不仅仅是一款游戏，更是一个很好的算法练习平台，你可以通过该游戏来学习和研究游戏算法问题。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值