spoj periodni

题解：

dp

方程弄出来就好做了

代码：

#include<bits/stdc++.h>
const int N=,M=;
typedef int arr[N];
typedef long long ll;
int n,k,c[N][N],h[N];
arr g;
void add(int &x,int y){x+=y;if (x>=M)x-=M;}
void dfs(int l,int r,int H,arr &ans)
{
    memset(ans,,sizeof(ans));
    ans[]=;
    int p=l,m=r-l+;
    for (int i=l;i<=r;i++)
     if (h[i]<h[p]) p=i;
    arr L,R;
    L[]=R[]=;
    int h0=H;H=h[p]-H;
    if (l<p)dfs(l,p-,h[p],L);
    if (p<r)dfs(p+,r,h[p],R);
    memset(g,,sizeof(g));
    for (int i=;i<=p-l;i++)
     for (int j=m;j>=i;j--)
      (g[j]+=(ll)ans[j-i]*L[i]%M)%=M;
    memcpy(ans,g,sizeof(ans));
    memset(g,,sizeof(g));
    for (int i=;i<=r-p;i++)
     for (int j=m;j>=i;j--)
      (g[j]+=(ll)ans[j-i]*R[i]%M)%=M;
    memcpy(ans,g,sizeof(ans));
    memset(g,,sizeof(g));
    ll P=;
    for (int i=;i<=m&&i<=H;i++)
     {
        for (int j=m;j>=i;j--)
         (g[j]+=(ll)ans[j-i]*c[m-j+i][i]%M*P%M)%=M;
        P=P*(H-i)%M;
     }
    memcpy(ans,g,sizeof(ans));
}
int main()
{
    scanf("%d%d",&n,&k);
    for (int i=;i<=n;i++)scanf("%d",&h[i]);
    c[][]=;
    for (int i=;i<=;i++)
     {
        c[i][]=;
        for (int j=;j<=i;j++)
         c[i][j]=(c[i-][j]+c[i-][j-])%M;
     }
    arr ans;
    dfs(,n,,ans);
    printf("%d\n",ans[k]);
    return ;
}