TopCoder SRM 489 Div1 Lev3:AppleTree
挺优秀的一道题,想出做法时有些惊艳。
题意:
数轴上有\(D\)个连续整数刻度,有\(N\)棵树要种在这些刻度上,其中第\(i\)棵与两旁(如果有的话)相邻的树至少要相距\(R_i\),问方法数。
\(1 \leq N , R_i \leq 40\)
思路:
首先,如果确定了种树的顺序,就确定了相邻树的最小间距。把\(D\)减掉最小间距之和,所得的就是“冗余刻度”的数量。
把这个数量分配给\(N+1\)段间隙,用插板法可以求出方法数。
所以问题在于,对于每一个\(L\),求出1到\(N\)的排列\(P\)的数量,满足:
\]
注意到,对于使\(R_i\)最大的\(i\),它的两侧种的是什么树,不影响这两段间隙的最小长度。
根据套路,这个时候我们应该在\(i\)的位置把排列割开并分别处理。
对于一个1到\(N\)的子集的长度为\(l\)的排列\(P\),定义其代价为:\(\sum_{i=1}^{l-1} \mathrm{max}(R_{P_i}, R_{P_{i+1}})\)。
我们先把\(R\)数组从小到大排序,接着DP:\(dp[i][j][k]\)表示,1到\(i\)这\(i\)个数,组成了\(j\)个不相交排列,排列的代价总和为\(k\)的方法数。
转移时考虑第\(i\)个数在一个排列的两端还是中间,删除之,加上第\(i\)个数两侧的空隙长度并转移即可。把树排序的意义在于,排序后第\(i\)棵树一定是\([1,i]\)中\(R\)值最大的那个。
注意到\(i,j \leq 40\),\(k \leq 1600\),故复杂度没有问题。
等等,正确性有问题……
因为长度为1的排列只有一个端点,而在先前的状态中很难区分长度为1和大于1的排列,故无法求出转移时乘的系数。
所以不妨直接考虑从长度为1的那些排列入手,发现……可以用它们合并成更长的排列?
于是想到可以DP状态不变,转移做如下修改:
首先,初始状态\(dp[0][0][0]=1\)。
其次,允许通过以下方式转移:
1、把\(i+1\)单独组成一个长为1的排列
2、任选两个排列,把\(i\)放在中间将其连接
3、任选一个排列,把\(i\)吸附在其端点旁边
这样把\(i+1\)添加进当前的排列集合之后,就可以转移到\(i'=i+1\)的情况。
对于每一个\(L\),求出1到\(N\)的排列\(P\)的数量,满足\(\sum_{i=1}^{N-1} \mathrm{max}(R_{P_i}, R_{P_{i+1}})=L\)
而对于每一个\(L\),要求的值就是\(dp[N][1][L]\)。Bingo!
代码:
// BEGIN CUT HERE
// END CUT HERE
#line 5 "AppleTrees.cpp"
#include <bits/stdc++.h>
using namespace std;
#define iinf 2000000000
#define linf 1000000000000000000LL
#define ulinf 10000000000000000000ull
#define MOD 1000000007LL
#define lson(v) ((v)<<1)
#define rson(v) (((v)<<1)^1)
#define mpr make_pair
typedef long long LL;
typedef unsigned long long ULL;
typedef unsigned long UL;
typedef unsigned short US;
typedef pair < int , int > pii;
clock_t __stt;
inline void TStart(){__stt=clock();}
inline void TReport(){printf("\nTaken Time : %.3lf sec\n",(double)(clock()-__stt)/CLOCKS_PER_SEC);}
template < typename T > T MIN(T a,T b){return a<b?a:b;}
template < typename T > T MAX(T a,T b){return a>b?a:b;}
template < typename T > T ABS(T a){return a>0?a:(-a);}
template < typename T > void UMIN(T &a,T b){if(b<a) a=b;}
template < typename T > void UMAX(T &a,T b){if(b>a) a=b;}
int n,dp[42][42][1655],fac[100005];
int inv(int val){
int ret=1,tms=MOD-2;
while(tms){
if(tms&1) ret=((LL)ret*(LL)val)%MOD;
val=((LL)val*(LL)val)%MOD;
tms>>=1;
}
return ret;
}
int C(int n,int m){
if(n<m) return 0;
if(!n) return 1;
int u=fac[n],d=((LL)fac[m]*(LL)fac[n-m])%MOD;
return ((LL)u*(LL)inv(d))%MOD;
}
class AppleTrees {
public:
int theCount(int D, vector <int> r){
sort(r.begin(),r.end());
n=(int)r.size();
int i,j,k,res=0;
fac[0]=1;
for(i=1;i<=D;++i){
fac[i]=((LL)fac[i-1]*(LL)i)%MOD;
}
dp[0][0][0]=1;
for(i=0;i<n;++i){
for(j=0;j<=n;++j){
for(k=0;k<=1620 && k<=D;++k){
if(!dp[i][j][k]) continue;
dp[i+1][j+1][k]+=dp[i][j][k];
dp[i+1][j+1][k]%=MOD;
dp[i+1][j][k+r[i]]+=((LL)dp[i][j][k]*2LL*(LL)j)%MOD;
dp[i+1][j][k+r[i]]%=MOD;
if(j>1){
dp[i+1][j-1][k+2*r[i]]+=((LL)dp[i][j][k]*(LL)j*(LL)(j-1))%MOD;
dp[i+1][j-1][k+2*r[i]]%=MOD;
}
}
}
}
for(i=1;i<=1620 && i<=D;++i){
res+=((LL)dp[n][1][i-1]*(LL)C(D-i+n,n))%MOD;
res%=MOD;
}
return res;
}
};
TopCoder SRM 489 Div1 Lev3:AppleTree的更多相关文章
- Topcoder Srm 726 Div1 Hard
Topcoder Srm 726 Div1 Hard 解题思路: 问题可以看做一个二分图,左边一个点向右边一段区间连边,匹配了左边一个点就能获得对应的权值,最大化所得到的权值的和. 然后可以证明一个结 ...
- Topcoder SRM 643 Div1 250<peter_pan>
Topcoder SRM 643 Div1 250 Problem 给一个整数N,再给一个vector<long long>v; N可以表示成若干个素数的乘积,N=p0*p1*p2*... ...
- topcoder srm 714 div1
problem1 link 倒着想.每次添加一个右括号再添加一个左括号,直到还原.那么每次的右括号的选择范围为当前左括号后面的右括号减去后面已经使用的右括号. problem2 link 令$h(x) ...
- Topcoder SRM 602 div1题解
打卡- Easy(250pts): 题目大意:rating2200及以上和2200以下的颜色是不一样的(我就是属于那个颜色比较菜的),有个人初始rating为X,然后每一场比赛他的rating如果增加 ...
- Topcoder SRM 627 div1 HappyLettersDiv1 : 字符串
Problem Statement The Happy Letter game is played as follows: At the beginning, several players ...
- Topcoder SRM 584 DIV1 600
思路太繁琐了 ,实在不想解释了 代码: #include<iostream> #include<cstdio> #include<string> #include& ...
- TopCoder SRM 605 DIV1
604的题解还没有写出来呢.先上605的. 代码去practice房间找. 说思路. A: 贪心,对于每个类型的正值求和,如果没有正值就取最大值,按着求出的值排序,枚举选多少个类型. B: 很明显是d ...
- topcoder srm 575 div1
problem1 link 如果$k$是先手必胜那么$f(k)=1$否则$f(k)=0$ 通过对前面小的数字的计算可以发现:(1)$f(2k+1)=0$,(2)$f(2^{2k+1})=0$,(3)其 ...
- topcoder srm 635 div1
problem1 link 首先枚举长度$L$.然后计算每一段长度$L$的差值最大公约数,然后差值除以最大公约数的结果可以作为当前段的关键字.然后不同段就可以比较他们的关键字,一样就是可以转化的. p ...
随机推荐
- 11GR2 双节点RAC 配置单节点DG
只记录主要步骤,供大家参考: RAC 搭建单节点 DG 1 修改源数据库开启归档和force loggingalter system set shared_servers=0; alter datab ...
- vue 单页应用中微信支付的坑
vue 单页应用中微信支付的坑 标签(空格分隔): 微信 支付 坑 vue 场景 在微信H5页面(使用 vue-router2 控制路由的 vue2 单页应用项目)中使用微信 jssdk 进行微信支付 ...
- 关于git的使用
一.关于GIT Git --- The stupid content tracker, 傻瓜内容跟踪器.Linus Torvalds 是这样给我们介绍 Git 的. Git 是用于 Linux内核 ...
- table中设置tr行间距
CSS border-collapse 属性设置表格的边框是否被合并为一个单一的边框 值 描述 separate 默认值.边框会被分开.不会忽略 border-spacing 和 empty-cell ...
- 异常:org.hibernate.id.IdentifierGenerationException
在有关联关系的数据表中保存数据时,先保存一端,再保存多端的抛出的异常(此时不管一端,还是多端的对象都没有设置id,属性,也就是要保存的两个对象的id 属性为空.) org.hibernate.id.I ...
- FFT抄袭笔记
你看我都不好意思说是学习笔记了,毕竟\(FFT\)我怎么可能学得会 那就写一篇抄袭笔记吧ctrl+c真舒服 先从多项式说起吧 1.多项式 我们定义一个多项式 \[F(x)=\sum_{i=0}^{n- ...
- 【[AHOI2012]树屋阶梯】
卡特兰数! 至于为什么是卡特兰数,就稍微说那么一两句吧 对于一个高度为\(i\)的阶梯,我们可以在左上角填一个高度为\(k\)的阶梯,右下角填一个高度为\(i-1-k\)的阶梯剩下的我们用一个大的长方 ...
- 【洛谷P3205】[HNOI2010]CHORUS 合唱队
合唱队 区间DP f[l][r][0/1]表示生成目标序列中的区间[l,r],最后一个数是a[l]/a[r] 的方案数 边界: f[i][i][0]=1 转移: f[l][r][0]=(a[l]< ...
- 前端静态文件如何应对HTTPS的到来
近几年,越来越多的网站开始支持https,我们可以看到国外的比如github.谷歌.facebook:国内的有百度.淘宝.博客园.coding.net.worktile等一系列的网站. 我再最近的开发 ...
- insertAdjacentHTML与innerHTML
insertAdjacentHTML:insertAdjacentHTML() 将指定的文本解析为HTML或XML,并将结果节点插入到DOM树中的指定位置.它不会重新解析它正在使用的元素,因此它不会破 ...