dp学习笔记 5 序列删除(区间dp)

题目
题意: 给定n个数，把除a1和an的其余数全部删除。删除第i个数的代价为a[i]*a[l]*a[r]. a[l]和a[r]是a[i]左右尚未被删除的数。
思路: 区间dp，一直在想先删除哪边另一边就会变，绕不明白。只能说对区间dp的理解不是很到位。
对于区间[i,j]枚举中间点k，k即区间[i,j]中最后删除的数，他的代价自然容易计算。
f[i][j] = min(f[i][j],f[i][k-1]+f[k+1][j]+a[k]a[i-1]a[j+1])
！这里要注意，直接把数组赋值为INF会寄，因为在转移式中存在非法区间，非法区间代价是0，但是却记为INF。可以特判边界，因为只会在k=i或者k=j时存在非法边界，也可以把所有非法区间赋值为0。
时间复杂度: O(nnn)
代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int a[502];
int f[502][502];
void solve()
{
   memset(f,0x3f,sizeof(f));
   cin>>n;
   for(int i=1;i<=n;++i) cin>>a[i];
   for(int i=2;i<=n;++i)
   for(int j=0;j<i;++j) f[i][j] = 0;
   
   for(int len=1;len<=n;++len)
   {
   	   for(int i=2,j=i+len-1;j<=n-1;++i,j=i+len-1)
   	   {
   	   	   // if(len==1) f[i][i] = a[i]*a[i-1]*a[i+1];
   	   	   for(int k=i;k<=j;++k) //枚举最后抽的哪个数
   	   	   {
   	   	   	  f[i][j] = min(f[i][j],f[i][k-1]+f[k+1][j]+a[k]*a[i-1]*a[j+1]);
   	   	   }
   	   }
   }
   cout<<f[2][n-1];
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}