写在前面
从接触到DP的第一天起就一直想知道很多转移方程为何这么写状态,每次看到都觉得很奇妙的话那一定是自己太垃圾了。我想是因为自己理论知识不到位所以不能站在一定的高度看待问题,但每次找论文有时竟然深入到最优化理论这种程度,又实在是静不下心来研究那些数学公式理论,所以渐渐明白做题和搞理论不一样,理论是搞那些总的规律,而做题是从具体问题中提取特殊的信息,具体问题具体分析,虽然也有规律,但绝不会是一劳永逸的规律。
这几篇DP试图探讨一些DP的规律和思想
hdu6578 Blank【2019多校赛第一场】
传送带
【思考】
这部分是我的思考过程可以跳过
分析问题,总结一句话就是:区间[l,r]中不同数恰好为x个的方案总数
直接归约问题的话这题就很像一个区间相关的题【什么线段树、莫队,总之自己学的不精就会想一大堆】,但这题有很多附加信息可以使题目变得特殊化,梳理一下:
1)这题并不是给定序列再进行区间上的增删改查,而是一个计数问题
2)这个区间很小,最多是100
3)能填的数很少,只能是1,2,3,4
4)要满足m个区间的条件交集,m最多为100
第一个信息就很显然说明解法不是线段树之类的,而且虽然看到第二个信息第一反应是区间dp,但这题不是求最优解类型的问题,而是计数问题,有点要构造的意思,第三第四两个信息为构造又营造了良好的时空复杂度,所以可以通过高维(n小,就像在对你说:你看我这么小你怎么能不多加几维呢)的计数dp来求解。
怎么构造呢,首先由于dp有很多种优化方式,所以我们能接受的空间范围大概在(空间上限*某个题目数据)之内。本题可以从能填的这4个数入手,也能从m个条件入手,还能从区间的长度入手,如果从m个条件入手,也就是满足前i个条件的数量,会用到容斥原理之类的,但这个条件不好表示,不像概率期望那样加加减减就行了,不好想,遂放弃;再看从4个数入手,有区间的限制条件,但也有些信息是冗余的并不需要知道,举个栗子,1 1 2 1和1 3 1 1,对于区间[2,3]都有贡献,对于区间[3,4]只有前者有贡献,所以数的排列和数的种类不重要,重要的是随便哪个数要至少出现一次,将这个问题归约为随便哪个数最后一次出现的位置在l到r之间,而这“随便哪个数”最多也只有四种情况,所以愉快地开个四位数组吧!
我才不会告诉你我是看了题解之后才凹的思路,不过也大差不差吧
【题解】
记录最后一次出现位置i j k cur,i<j<k<cur(限制顺序可以优化时间),并且cur可以通过滚动数组实现
考虑第当前cur位置放哪一个数f[i][j][k][cur-1]会转移到哪
放第一个数可以转移到f[j][k][cur-1][cur]
放第二个数可以转移到f[i][k][cur-1][cur]
放第三个数可以转移到f[i][j][cur-1][cur]
放第四个数可以转移到f[i][j][k][cur]
由于本题的区间操作没有修改操作,所以对于区间的限制可以将循环变量作为右端点,只对左端点进行判断
【代码】
注意取模的写法,容易写T
#include <bits/stdc++.h>
using namespace std;
#define mp make_pair
const int MAXN=101;
const int P=998244353;
int T,n,m,l,r,x;
int f[MAXN][MAXN][MAXN][3],ans;
vector<pair<int,int> >d[MAXN];
void init(){
ans=0;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
d[i].clear();
d[i].push_back(mp(i,1));
}
for(int i=0;i<m;i++){
scanf("%d%d%d",&l,&r,&x);
d[r].push_back(mp(l,x));
}
memset(f,0,sizeof(f));
f[0][0][0][0]=1;
for(int cur=1;cur<=n;cur++){
int c=cur&1;
for(int i=0;i<=cur;i++)
for(int j=i;j<=cur;j++)
for(int k=j;k<=cur;k++)
f[i][j][k][c]=0;
for(int i=0;i<=cur;i++)
for(int j=i;j<=cur;j++)
for(int k=j;k<=cur;k++){
(f[j][k][cur-1][c]+=f[i][j][k][c^1])%=P;
(f[i][k][cur-1][c]+=f[i][j][k][c^1])%=P;
(f[i][j][cur-1][c]+=f[i][j][k][c^1])%=P;
(f[i][j][k][c]+=f[i][j][k][c^1])%=P;
}
for(int i=0;i<=cur;i++)
for(int j=i;j<=cur;j++)
for(int k=j;k<=cur;k++)
for(auto pi:d[cur]){
l=pi.first;
r=cur;
x=pi.second;
if((i>=l)+(j>=l)+(k>=l)+(cur>=l)!=x) f[i][j][k][c]=0;
}
}
for(int i=0;i<=n;i++)
for(int j=i;j<=n;j++)
for(int k=j;k<=n;k++)
(ans+=f[i][j][k][n&1])%=P;
printf("%d\n",ans);
}
int main(){
//freopen("hdu6578.in","r",stdin);
//freopen("hdu6578.out","w",stdout);
scanf("%d",&T);
while(T--){
init();
}
return 0;
}