P5241 序列(滚动数组+前缀和优化dp)
挺神仙的一题
看看除了dp好像没什么其他办法了
想着怎么构个具体的图出来,然鹅不太现实。
于是我们想办法用几个参数来表示dp数组
加了几条边肯定要的吧,于是加个参数$i$表示已加了$i$条边
这显然是不够的。于是我们又想:强连通分量.....连通块.......
于是加个$j$表示还有$j$个强连通分量
于是dp数组为$f[i][j]$
这是我们发现一个问题,状态$f[i][j]$不一定是合法的。
那dp不就GG了吗
再次撕烤,我们发现每次加上的边无非就3种情况:
1.把2个强连通分量(或链)连成一条链
2.在某个强连通分量中瞎连(没啥用)
3.在1条链上的某点向回连,形成一个环,缩成一个新强连通分量(可以减少任意个强连通分量)
我们设$k-1$条边(dp数组下标$k$为正数较好处理)投入到第3种情况
要生成剩下$j$个强连通的情况,我们最少投入$n-j$条边用于第1种情况
所以$n-j+(k-1)<=i$
我们又发现,要生成剩下$j$个强连通的情况,我们最多共投入的边数$i$是有限制的
最多情况就是1个块有$n-j+1$个点,剩下$j-1$个块只有1个点,蓝后大块每个点连$n-1$条边,小块互相之间弱连通
那么最大边数为$(n-j+1)*(n-1)+(j-2+j-3+j-4+...+1)=(n-j+1)*(n-1)+(j-1)*(j-2)/2$
所以$i<=(n-j+1)*(n-1)+(j-1)*(j-2)/2$
总结一下,即设$f[i][j][k]$表示到第$i$条边,有$j$个强连通分量,$k-1$条边向回连的方案数
限制条件:
$n-j+(k-1)<=i$
$i<=(n-j+1)*(n-1)+(j-1)*(j-2)/2$
转移:
$f[i][j][k]+=f[i-1][j][k]$(第2种情况)
$f[i][j][k]+=\sum_{h=j+1}^{n}f[i-1][h][k-1]$
显然是可以滚动数组+前缀和优化的辣
然鹅复杂度还是太高,主要因为k很麻烦
仔细观察k,发现
$n-j+(k-1)<=i$
$k<=i+j-n+1$
发现$i>=2n$时k总是合法的
于是我们就可以愉快地缩成2维辣
#include<iostream>
#include<cstdio>
#include<cstring>
#define rint register int
using namespace std;
inline int Min(int a,int b){return a<b?a:b;}
const int mod=1e9+;
inline int Md(int x){return x<mod?x:x-mod;}
#define N 405
int n,f[][N][N],sf[][N][N],g[][N],sg[N][N],lim[N],ans[N*N];
int main(){
scanf("%d",&n); int tn=Min(n*(n-),n<<),w=;
for(rint j=;j<=n;++j) lim[j]=(n-j+)*(n-)+(j-)*(j-)/;
f[][n][]=ans[]=;
for(rint j=;j<=n;++j) sf[][n][]=;
for(rint i=;i<=tn;++i,w^=){
for(rint j=;j<=n;++j)
for(rint k=;k<=n;++k)
f[w][j][k]=;
for(rint j=;j<=n;++j) if(lim[j]>=i)
for(rint k=;k<=n;++k) if(i-(k-)>=n-j)
f[w][j][k]=Md(f[w^][j][k]+sf[w^][j+][k-]);
for(rint j=n;j;--j)
for(rint k=;k<=n;++k){
sf[w][j][k]=Md(sf[w][j+][k]+f[w][j][k]);
ans[i]=Md(ans[i]+f[w][j][k]);
}
}w=;
for(rint j=;j<=n;++j)
for(rint k=;k<=n;++k)
g[][j]=Md(g[][j]+f[][j][k]);
for(rint j=n;j;--j) sg[][j]=Md(sg[][j+]+g[][j]);//降维
for(rint i=tn+;i<=n*(n-);++i,w^=){
for(rint j=;j<=n;++j) g[w][j]=;
for(rint j=;j<=n;++j) if(lim[j]>=i)
g[w][j]=Md(g[w^][j]+sg[w^][j+]);
for(rint j=n;j;--j){
sg[w][j]=Md(sg[w][j+]+g[w][j]);
ans[i]=Md(ans[i]+g[w][j]);
}
}
for(rint i=;i<=n*(n-);++i) printf("%d ",ans[i]);
return ;
}
最新文章
- ASP.NET Core 中文文档 第二章 指南(4.2)添加 Controller
- zend studio 快捷键
- 在wex5平台grid显示问题
- LVS+Keepalived搭建MyCAT高可用负载均衡集群
- neXtep 安装过程整理
- PBOC金融IC卡,卡片与终端交互的13个步骤,简介-第三组
- uC/OS-II任务(OS_task)块
- CSAPP(前言)
- [leetcode]Combine Two Tables
- linux中$与()的一点使用疑惑解释
- [Java] 转:多线程 (并发)总结
- SQL中存储过程和自定义函数的区别
- Oracle 多版本控制
- Delphi 生成excel中的饼图
- hdu 4885 (n^2*log(n)推断三点共线建图)+最短路
- mybatis sql in not in的使用
- linux 软件安装篇
- js数组乱序输出 数组乱序排列
- 1个比较简单的使用java反射机制获取前台数据进行数据封装的例子
- express 调优的一个过程和心得,不错的文章
热门文章
- 对 data属性的使用之一
- 如何写一个优秀的GitHub项目README文档?
- [LeetCode] 78. Subsets tag: backtracking
- 配置完centos 6以后,大概需要安装的软件(主要是yum)
- java类的包装类
- json_decode 转数组
- 使用vscode编译sass
- BCB Access violateion at Address 0000 0003. Read of address 0000 0003
- jsp重新打开一个新的页面
- Unity shader学习之逐像素漫反射光照模型