1049. 大盗阿福
每个店铺有两种状态,$f[i][0] $表示只选择前
i
i
i 个店铺,第
i
i
i 个店铺不洗劫,
f
[
i
]
[
1
]
f[i][1]
f[i][1] 表示选择前
i
i
i 个店铺,并洗劫第
i
i
i 个店铺。
f
[
i
]
[
0
]
=
m
a
x
(
f
[
i
−
1
]
[
1
]
,
f
[
i
−
1
]
[
0
]
)
f
[
i
]
[
1
]
=
f
[
i
−
1
]
[
0
]
+
w
[
i
]
f[i][0] = max(f[i-1][1],f[i-1][0])\\ f[i][1] = f[i-1][0]+w[i]
f[i][0]=max(f[i−1][1],f[i−1][0])f[i][1]=f[i−1][0]+w[i]
解题时存在的问题:
我发现有的时候我总是搞不懂dp到底应该怎么去初始化?
1A,发现这道题根本就不用初始化啊
y总的初始化,一开始状态机位置状态0。
f[0][0] = 0,f[0][1] = -inf;(不可到达)
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1e5+10;
int t,n,a[N];
int f[N][2];
int main(){
cin>>t;
while(t--){
for(int i=1;i<=n;i++){
f[i][0]=f[i][1]=0;
}
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i];
}
for(int i=1;i<=n;i++){
f[i][0]=max(f[i-1][1],f[i-1][0]);
f[i][1]=f[i-1][0]+a[i];
}
cout<< max(f[n][0],f[n][1]) <<endl;
}
return 0;
}
1057. 股票买卖 IV
一个做题时一定会想到的问题就是,为什么我做这道题的时候,要用dp来求解呢?
读完题之后,能够得到的一个非常直觉的东西就是,股票的买卖是有”后效性“的,也就是我当前的决策可能会受到之后的决策的影响。
这道题我是如何证明他是没有后效性的呢?
试着做一做。
f [ i ] [ j ] [ 0 ] f[i][j][0] f[i][j][0] 表示从前 i i i 天 中恰好进行 j j j 笔交易(定义卖出去一次才算进行了一笔交易),当前没有持股。 f [ i ] [ j ] [ 1 ] f[i][j][1] f[i][j][1] 表示从前 i i i 天 中恰好进行 j j j 笔交易(定义卖出去一次才算进行了一笔交易),当前有持股。
很明显我的状态定义没有问题
考虑一个问题,第二维度
j
j
j 是定义成卖出才算一笔交易还是买入是一笔交易呢?很明显卖出之后才能再买入,才算进行了一笔完整的交易,然后我TM就把状态转移方程给写错了
"
错
误
"
的
状
态
转
移
方
程
。
f
[
i
]
[
j
]
[
0
]
=
m
a
x
(
f
[
i
−
1
]
[
j
]
[
0
]
,
f
[
i
−
1
]
[
j
−
1
]
[
1
]
+
a
[
i
]
)
f
[
i
]
[
j
]
[
1
]
=
m
a
x
(
f
[
i
−
1
]
[
j
]
[
0
]
−
a
[
i
]
,
f
[
i
−
1
]
[
j
]
[
1
]
)
"
正
确
"
的
状
态
转
移
方
程
。
f
[
i
]
[
j
]
[
0
]
=
m
a
x
(
f
[
i
−
1
]
[
j
]
[
0
]
,
f
[
i
−
1
]
[
j
]
[
1
]
+
a
[
i
]
)
f
[
i
]
[
j
]
[
1
]
=
m
a
x
(
f
[
i
−
1
]
[
j
]
[
0
]
−
a
[
i
]
,
f
[
i
−
1
]
[
j
−
1
]
[
1
]
)
"错误"的状态转移方程。\\ f[i][j][0] = max(f[i-1][j][0],f[i-1][j-1][1]+a[i])\\ f[i][j][1] = max(f[i-1][j][0]-a[i],f[i-1][j][1])\\ "正确"的状态转移方程。\\ f[i][j][0] = max(f[i-1][j][0],f[i-1][j][1]+a[i])\\ f[i][j][1] = max(f[i-1][j][0]-a[i],f[i-1][j-1][1])\\
"错误"的状态转移方程。f[i][j][0]=max(f[i−1][j][0],f[i−1][j−1][1]+a[i])f[i][j][1]=max(f[i−1][j][0]−a[i],f[i−1][j][1])"正确"的状态转移方程。f[i][j][0]=max(f[i−1][j][0],f[i−1][j][1]+a[i])f[i][j][1]=max(f[i−1][j][0]−a[i],f[i−1][j−1][1])
看完其他大佬的博客才明白了自己哪"错"了,其实我状态定义和状态的转移都没错,但是我写代码的时候数组越界了。
a n s = m a x ( f [ n ] [ K ] [ 0 ] , f [ n ] [ K ] [ 1 ] ) , K ∈ [ 1 , k ] ans = max(f[n][K][0],f[n][K][1]),K \in [1,k] ans=max(f[n][K][0],f[n][K][1]),K∈[1,k]
初始化好像也是没有啊。
f[0][0][0]=0;
f[0][0][1]=-inf;
能感觉到递推公式有问题,但是说不出来。
Q1 你会打dfs吗?
没试过啊
#include <iostream>
#include <cstring>
#include <algorithm>
const int N = 1e5+10;
typedef long long ll;
using namespace std;
ll n,k,f[N][110][2],a[N];
int main(){
cin>>n>>k;
for(int i=1;i<=n;i++){
cin>>a[i];
}
for(int i=1;i<=n;i++){
for(int j=1;j<=k;j++){
// 应该从0开始dp,并且不能越界。
f[i][j][0] = max(f[i-1][j][0],f[i-1][j-1][1]+a[i]);
f[i][j][1] = max(f[i-1][j][0]-a[i],f[i][j][1]); // 感觉就是错的。。。
}
}
ll ans=0;
for(int i=0;i<=k;i++){
ans = max(ans,f[n][i][0]);
}
cout<<ans<<endl;
return 0;
}
看完了y总的讲解之后,发现我和y总的答案非常非常接近了,状态表示想对了,状态转移方程几乎是对的,答案输出也是对的,初始化没有想到。
#include <iostream>
#include <cstring>
#include <algorithm>
#define debug(x) cout<<#x<<":"<<x<<endl
const int N = 1e5+10;
typedef long long ll;
using namespace std;
ll n,k;
int f[N][110][2],a[N];
int main(){
cin>>n>>k;
for(int i=1;i<=n;i++){
cin>>a[i];
}
memset(f,-0x3f,sizeof f);
for(int i=0;i<=n;i++){
f[i][0][0]=0;
}
for(int i=1;i<=n;i++){
for(int j=0;j<=k;j++){//从0开始的
f[i][j][0] = f[i-1][j][0];
if(j)//防止j-1越界
f[i][j][0] = max(f[i][j][0],f[i-1][j-1][1]+a[i]);//卖出了算完成了整个交易。
f[i][j][1] = max(f[i-1][j][0]-a[i],f[i-1][j][1]); //
}
}
int ans=0;
for(int i=0;i<=k;i++){
ans = max(ans,f[n][i][0]);
}
cout<<ans<<endl;
return 0;
}
1058. 股票买卖 V
比上一题增加了一个约束条件,卖出股票之后有一天的冷冻期,相当于又多了一个维度的信息。
我直接考虑的是在上一题的前提下再增加一个维度的信息,表示进行到这个状态的上一个状态是不是卖出物品。
应该也能得到结果,但是感觉写起来很晕qwq
我TM才发现这题没有K这个限制了,我是SB
写出状态机机会,寻找状态机的入口和出口,然后写状态转移方程。
#include <iostream>
#include <cstring>
#include <algorithm>
#define debug(x) cout<<#x<<":"<<x<<endl
const int N = 1e5+10;
typedef long long ll;
using namespace std;
ll n;
int f[N][3],a[N];//最后一个维度记录上一次操作是卖出还是持仓
int main(){
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i];
}
memset(f,0xcf,sizeof f);
for(int i=0;i<=n;i++)f[i][2]=0;
for(int i=1;i<=n;i++){
f[i][0] = f[i-1][0];
if(i>=2){
// 麻了,明明是f[i-1],我怎么想的?
f[i][0] = max(f[i][0],f[i-2][2]-a[i]);
}
f[i][1]=f[i-1][0]+a[i];
f[i][2]=max(f[i-1][2],f[i-1][1]);
}
int ans=0;
for(int i=0;i<=n;i++){
ans = max(ans,max(f[i][1],f[i][2]));
}
cout<<ans<<endl;
return 0;
}
#include <iostream>
#include <cstring>
#include <algorithm>
#define debug(x) cout<<#x<<":"<<x<<endl
const int N = 1e5+10;
typedef long long ll;
using namespace std;
ll n;
int f[N][3],a[N];
int main(){
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i];
}
memset(f,0xcf,sizeof f);
//边界
f[0][2]=0;
for(int i=1;i<=n;i++){
f[i][0] = f[i-1][0];
f[i][0] = max(f[i][0],f[i-1][2]-a[i]);
f[i][1]=f[i-1][0]+a[i];
f[i][2]=max(f[i-1][2],f[i-1][1]);
}
cout<<max(f[n][1],f[n][2])<<endl;
return 0;
}
1052. 设计密码
计数类dp,发现时间复杂度完全不会分析啊。没有思路。
有串匹配,考虑KMP, f [ i ] [ j ] f[i][j] f[i][j] 表示已经填写到第 i i i 个字符,且 T 串和 S 串匹配到 j j j 的所有方案。
f [ i ] [ 0 ] f[i][0] f[i][0] 表示 T 串 和 S串没有匹配上。
状态机是什么呢?
递推时考虑用 f [ i ] f[i] f[i] 更新 f [ i + 1 ] f[i+1] f[i+1], 对于每个 i i i ,都有26个小写英文字母(26个状态),每个状态有唯一的对应 T T T 串中 j j j 的变化。
第一维循环 i i i , 第二维循环 j j j 枚举 i i i 匹配到的位置,第三维循环 i i i 所有可能的状态(26个小写字母), 对于每个小写字母,求出对应的 j j j 的状态。如果 j j j 不能匹配上 S S S 代表, S S S 串中没有 T T T 串作为子串,累加。
初始状态 f [ 0 ] [ 0 ] = 1 ( 计 数 类 d p ) f[0][0]=1(计数类dp) f[0][0]=1(计数类dp)
时间复杂度 O ( n ∗ m ∗ 26 ∗ 26 ) O(n*m*26*26) O(n∗m∗26∗26)
可以预处理出每个字母 a
的最后匹配到的
j
j
j 状态。
时间复杂度 O ( n ∗ m ∗ 26 ) O(n*m*26) O(n∗m∗26)
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 60, mod = 1e9+7;
char s[N],p[N];
int ne[N];
int n,m;
int f[N][N];
int main(){
cin>>n;
cin>>(p+1);
m=strlen(p+1);
for(int i=2,j=0;i<=m;i++){
while(j>0&&(p[i]!=p[j+1]))j=ne[j];
if(p[j+1]==p[i])j++;
ne[i]=j;
}
f[0][0]=1;
//我从第1位开始匹配的,感觉更自然一些。
//用以前推当前。
for(int i=1;i<=n;i++){
for(int k=0;k<m;k++){
for(char j='a';j<='z';j++){// i+1 位置的字符
int u=k;
while(u&&p[u+1]!=j)u=ne[u];
if(p[u+1]==j)u++;
if(u<m){
// cout<<"YES"<<endl;
f[i][u] = (f[i][u] + f[i-1][k]) %mod;
}
}
}
}
int res=0;
for(int i=0;i<m;i++){
res = (res+ f[n][i]) %mod;
}
cout<<res<<endl;
return 0;
}
如果用当前推以后。
f[0][0] = 1;
for (int i = 0; i < n; i ++ )
for (int j = 0; j < m; j ++ )
for (char k = 'a'; k <= 'z'; k ++ )
{
int u = j;
while (u && k != str[u + 1]) u = nxt[u];
if (k == str[u + 1]) u ++ ;
if (u < m) f[i + 1][u] = (f[i + 1][u] + f[i][j]) % mod;
}
预处理dp中KMP的过程,不预处理需要26ms,用了map花了80ms,qwq。unordered_map不知道为什么不行,应该是没重载比较符号
for(int k=0;k<m;k++){
for(char j='a';j<='z';j++){
int u=k;
while(u&&p[u+1]!=j)u=ne[u];
if(p[u+1]==j)u++;
pos[{k,j}]=u;
}
}
f[0][0]=1;
for(int i=1;i<=n;i++){
for(int k=0;k<m;k++){
for(char j='a';j<='z';j++){// i+1 位置的字符
int u = pos[{k,j}];
if(u<m){
f[i][u] = (f[i][u] + f[i-1][k]) %mod;
}
}
}
}
1053. 修复DNA
和上一道题很相似,唯一的不同就是用AC自动机维护字符串。