一种神奇的计数dp

前言：
在做计数题时，我们会一类要求满足条件的排列个数或极值的题目。
这时有一个神奇的dp可以派上一些用场。

1. 适用范围

• 排列的条件适用于相邻的两数。
• 贡献可以拆开计算。
• 一般贡献因相邻两数的大小关系而发生结构变化。

2. 具体算法

• 设一个$dp$数组$f[N][N]$，从小到大插入。 （$fij$表示考虑到第$i$个数，形成了$j$个连续段。 ）
• 在转移的时候不用考虑很多，最多因题目的特殊限制（比如头尾必须是啥）进行一些特判。
• 可以抽象为一段一段的线段，除非特殊条件限制，其他时候可形成线段，或者和其他两个线段和在一起。（这段写的太抽象了可以直接看例题）

3. 复杂度

显然为$O(n^2)$

例题

1. CF704BAnt Man

这个一看就符合于我们前面说的场景，直接说一下怎么转移吧，这道题的特殊条件是它给出了开头$s$和结尾$e$。

• 如果我们选择形成一个新线段，那s和e一定不能已经和成一个线段了。

  if(i+1<max(e,s)||j>1)
  

• 如果我们选择和两个线段和并，段数肯定大于1。

  if(j>1)
  

• 还有一种如果我们选择到一个线段头，那一定不能只有一段并且包含着整个序列的头。

  if(j>1||(j==1&&i+1<e))
  

• 还有一种如果我们选择到一个线段尾，那一定不能只有一段并且包含着整个序列的尾。

  (j>1||(j==1&&i+1<s))
  

最后给个代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define N 5001
#define inf 1e15
ll n,s,e;
ll x[N],a[N],b[N],c[N],d[N];
ll f[N][N];
int main()
{
	//freopen("test1.in","r",stdin);
	//freopen("test1.out","w",stdout);
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>s>>e;
	for(int i=1;i<=n;i++)cin>>x[i];
	for(int i=1;i<=n;i++)cin>>a[i];
	for(int i=1;i<=n;i++)cin>>b[i];
	for(int i=1;i<=n;i++)cin>>c[i];
	for(int i=1;i<=n;i++)cin>>d[i];
	for(int i=1;i<=n;i++)
	{
		a[i]+=x[i];b[i]-=x[i];c[i]+=x[i];d[i]-=x[i];
	}
	memset(f,0x3f,sizeof(f));
	f[0][0]=0;
	for(int i=0;i<n;i++)
	{
		if(i+1!=s&&i+1!=e)
		{
			for(int j=0;j<=i;j++)
			{
				if(i+1<max(e,s)||j>1)
				{
					f[i+1][j+1]=min(f[i+1][j+1],f[i][j]+b[i+1]+d[i+1]);
				}
				if(j>1)
				{
					f[i+1][j-1]=min(f[i+1][j-1],f[i][j]+a[i+1]+c[i+1]);
				}
				if(j>1||(j==1&&i+1<e))
				{
					f[i+1][j]=min(f[i+1][j],f[i][j]+a[i+1]+d[i+1]);
				}
				if(j>1||(j==1&&i+1<s))
				{
					f[i+1][j]=min(f[i+1][j],f[i][j]+b[i+1]+c[i+1]);
				}
			}
		}
		if(i+1==s)
		{
			for(int j=0;j<=i;j++)
			{
				f[i+1][j+1]=min(f[i+1][j+1],f[i][j]+d[i+1]);
				if(j>0)f[i+1][j]=min(f[i+1][j],f[i][j]+c[i+1]);
			}
		}
		if(i+1==e)
		{
			for(int j=0;j<=i;j++)
			{
				f[i+1][j+1]=min(f[i+1][j+1],f[i][j]+b[i+1]);
				if(j>0)f[i+1][j]=min(f[i+1][j],f[i][j]+a[i+1]);
			}
		}
	}
	cout<<f[n][1]<<endl;
	return 0;
}

2.loj2743 JOI 摩天大楼
这个太经典了就不分析了。（而且和上面一个差不多）
放个代码吧，转移看代码。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define N 101
#define L 1001
#define yu 1000000007
ll n,l,f[N][N*2][L][3],a[N];
int main()
{
	//freopen("1.txt","r",stdin);
	cin>>n>>l;
	for(int i=1;i<=n;i++)cin>>a[i];
	sort(a+1,a+n+1);
	memset(f,0,sizeof(f));
	f[0][0][0][0]=1;
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<=2*i&&j<n;j++)
		{
			for(int k=0;k<=l;k++)
			{
				for(int p=0;p<=2;p++)
				{
					ll nk=k+(2*j-p)*(a[i+1]-a[i]);
					if(!f[i][j][k][p]||nk>l)continue;
					f[i+1][j+1][nk][p]=(f[i+1][j+1][nk][p]+f[i][j][k][p]*(j+1-p))%yu;
					f[i+1][j][nk][p]=(f[i+1][j][nk][p]+f[i][j][k][p]*(2*j-p))%yu;
					f[i+1][j-1][nk][p]=(f[i+1][j-1][nk][p]+f[i][j][k][p]*(j-1))%yu;
					if(p==0)
					{
						f[i+1][j+1][nk][1]=(f[i+1][j+1][nk][1]+2*f[i][j][k][p])%yu;
						if(j)f[i+1][j][nk][1]=(f[i+1][j][nk][1]+2*f[i][j][k][p])%yu;
						if(!j)f[i+1][j+1][nk][2]=(f[i+1][j+1][nk][2]+f[i][j][k][p])%yu;
					}
					if(p==1)
					{
						f[i+1][j+1][nk][2]=(f[i+1][j+1][nk][2]+f[i][j][k][p])%yu;
						if(j)f[i+1][j][nk][2]=(f[i+1][j][nk][2]+f[i][j][k][p])%yu;
					}
				}
			}
		}
	}
	ll ans=0;
	for(int i=0;i<=l;i++)
	{
		ans+=f[n][1][i][2];
		ans%=yu;
	}
	cout<<ans<<endl;
	return 0;
}