题目
思路
一看到不能相交的线段,显然是区间 d p \tt dp dp 无疑,因为连接的点一定是区间。
只用 f ( l , r ) f(l,r) f(l,r) 显然转移不动,因为你没法把这个区间切割成一个更小的区间。所以我们加点佐料!
用 f ( l , r , 0 / 1 ) f(l,r,0/1) f(l,r,0/1) 表示 l , r l,r l,r 是否必须用边连接 情况下的方案数。若直接连接,则考虑断开这条边之后,形成的两个树。枚举这个划分点,方案为
f ( l , r , 1 ) = ∑ i = l + 1 r − 1 [ f ( l , i , 0 ) + f ( l , i , 1 ) ] [ f ( i + 1 , r , 0 ) + f ( i + 1 , r , 1 ) ] f(l,r,1)=\sum_{i=l+1}^{r-1}\Big[f(l,i,0)+f(l,i,1)\Big]\Big[f(i+1,r,0)+f(i+1,r,1)\Big] f(l,r,1)=i=l+1∑r−1[f(l,i,0)+f(l,i,1)][f(i+1,r,0)+f(i+1,r,1)]
若不直接连接,则考虑 l l l 向最右边的那个点连的边。方案为
f ( l , r , 0 ) = ∑ i = l + 1 r − 1 f ( l , i , 1 ) [ f ( i , r , 0 ) + f ( i , r , 1 ) ] f(l,r,0)=\sum_{i=l+1}^{r-1}f(l,i,1)\Big[f(i,r,0)+f(i,r,1)\Big] f(l,r,0)=i=l+1∑r−1f(l,i,1)[f(i,r,0)+f(i,r,1)]
代码
#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
inline int readint(){
int a = 0; char c = getchar(), f = 1;
for(; c<'0'||c>'9'; c=getchar())
if(c == '-') f = -f;
for(; '0'<=c&&c<='9'; c=getchar())
a = (a<<3)+(a<<1)+(c^48);
return a*f;
}
template < class T >
void getMax(T&a,T b){ if(a < b) a = b; }
template < class T >
void getMin(T&a,T b){ if(b < a) a = b; }
const int MaxN = 502;
const int Mod = 1000000007;
// 0: l,r不直接连接,连成树
// 1: l,r直接相连
int dp[MaxN][MaxN][2];
int a[MaxN][MaxN];
int main(){
int n = readint();
for(int i=1; i<=n; ++i)
for(int j=1; j<=n; ++j)
a[i][j] = readint();
for(int i=1; i<=n; ++i)
dp[i][i][0] = 1;
for(int len=1; len<n; ++len)
for(int i=1; i+len<=n; ++i){
dp[i][i+len][1] = // l没了
(dp[i+1][i+len][0]+
dp[i+1][i+len][1])%Mod;
for(int j=i+1; j<i+len; ++j){
dp[i][i+len][1] += 1ll*
(dp[i][j][1]+\
dp[i][j][0])*\
(dp[j+1][i+len][0]+\
dp[j+1][i+len][1])%Mod;
dp[i][i+len][0] +=
1ll*dp[i][j][1]*\
(dp[j][i+len][0]+\
dp[j][i+len][1])%Mod;
for(int t=0; t<2; ++t)
dp[i][i+len][t] %= Mod;
}
dp[i][i+len][1] *= a[i][i+len];
}
int ans = (dp[1][n][0]+dp[1][n][1])%Mod;
printf("%d\n",ans);
return 0;
}