【GDOI2017第二轮模拟day1】最长路径(性质题,容斥,组合数学)
来源:互联网 发布:vue.js 动态绑定事件 编辑:程序博客网 时间:2024/06/05 09:51
Description
在Byteland 一共有n 个城市,编号依次为1 到n,它们之间计划修建n(n-1)/2条单向道路,对于任意两个不同的点i 和j,在它们之间有且仅有一条单向道路,方向要么是i 到j,要么是j 到i。换句话说,这是一个n 个点的竞赛图。
Byteasar 居住在1 号城市,他希望从1 号城市出发,沿着单向道路不重复地访问一些城市,使得访问的城市数尽可能多。
请写一个程序,帮助Byteasar 计算有多少种道路修建方式,使得从1 号点出发的最长简单路径经过点数恰好为k,由于答案可能很大,请对P 取模输出。
Solution
这是一道性质题。考场的时候肯定是打表了,60分!!!
首先我们知道一些性质:
1、竞赛图如果是一个强连通分量,那么这个图肯定有哈密顿回路(可用数学归纳法证明,假设有一个k个点的满足这个性质,多加一个点……)
2、这个竞赛图经过强连通分量缩点之后会形成一条链(有一个贡献答案的主链和很多小链构成)
3、答案为1号节点在所在的强连通分量的点数+链上1的强连通点后面的点总数(这些是缩点之前的)
那么现在我们DP一下就好了。
我们设
g可以用容斥来做。
然后根据性质3,我们枚举i为1号所在的强连通分量的节点数量和之后的节点数量j
因为i枚举的时候包含了1,所以组合数要注意一下。
Code
#include<iostream> #include<stdio.h>#include<string.h>#include<algorithm>#include<math.h>#define fo(i,a,b) for(i=a;i<=b;i++)typedef long long ll;using namespace std;const int maxn=2e3+7;ll i,j,k,l,t,n,m,mo,ans[maxn];ll c[maxn][maxn],f[maxn],g[maxn],er[maxn*maxn];int main(){// freopen("path.in","r",stdin);// freopen("path.out","w",stdout); scanf("%d%d",&n,&mo); fo(i,1,n)c[i][0]=c[0][i]=1; fo(i,1,n)fo(j,1,i)c[i][j]=(c[i-1][j]+c[i-1][j-1])%mo;c[0][0]=1; er[0]=1;fo(i,1,n*n-1)er[i]=er[i-1]*2%mo; fo(i,1,n)f[i]=er[(i-1)*i/2]; f[0]=g[0]=1; fo(i,1,n){ g[i]=f[i]; fo(j,1,i-1) g[i]=(g[i]-c[i][j]*g[j]%mo*f[i-j]%mo)%mo; } fo(i,1,n){ fo(j,0,n-i){ ans[i+j]=(ans[i+j]+c[n-1][i-1]*c[n-i][j]%mo*g[i]%mo*f[n-i-j]%mo*f[j]%mo)%mo; } } fo(i,1,n){ ans[i]=(ans[i]+mo)%mo; printf("%lld\n",ans[i]); }}
- 【GDOI2017第二轮模拟day1】最长路径(性质题,容斥,组合数学)
- 【GDOI2017第二轮模拟day1】最长路径
- 【JZOJ5061】【GDOI2017第二轮模拟day1】最长路径
- jzoj 5061. 【GDOI2017第二轮模拟day1】最长路径 动态规划
- 【GDOI2017第二轮模拟day1】公路建设
- 【JZOJ5060】【GDOI2017第二轮模拟day1】公路建设
- GDOI2017第二轮模拟
- 【jzoj5060】【GDOI2017第二轮模拟day1】【公路建设】【数据结构】
- 【GDOI2017第二轮模拟day1】公路建设(克鲁斯卡尔最小生成树+线段树+归并)
- GDOI2017第二轮模拟总结
- GDOI2017模拟第二轮总结
- GDOI2017第二轮模拟总结
- jzoj 5060. 【GDOI2017第二轮模拟day1】公路建设 线段树+最小生成树
- GDOI2017模拟第二轮 4.15-4.17
- jzoj【GDOI2017第二轮模拟day2】开房间
- 【GDOI2017第二轮模拟day2】中位数
- 【GDOI2017第二轮模拟day2】开房间
- 【GDOI2017第二轮模拟day2】中位数
- 【数据压缩】Huffman编解码器
- SQLserver统计某个字段不重复的个数以及怎么以时间为单位取出一段时间的数据
- bzoj 4811: [Ynoi2017]由乃的OJ (树链剖分)
- 深入分析Java ClassLoader原理
- Mac终端切换root用户
- 【GDOI2017第二轮模拟day1】最长路径(性质题,容斥,组合数学)
- FreeBSD install KDE plasma5 and Flash
- 循环语句 for
- (转)诺亚财富巨亏5亿,“嗜血模式”走向末路?
- win10 uwp 改变鼠标
- Can not find the tag library descriptor for “http://java.sun.com/jsp/jstl/core”等
- open live writer's plugins
- 【动态规划】数字三角形最大和(二)(递推循环)
- Valid Palindrome