BZOJ 1996 [Hnoi2010]chorus 合唱队

最新推荐文章于 2019-03-01 13:34:00 发布

原创最新推荐文章于 2019-03-01 13:34:00 发布

· 297 阅读

1 ·

版权

动态规划专栏收录该内容

40 篇文章

订阅专栏

本文介绍了一种使用区间动态规划解决特定取数方案数问题的方法。通过定义f[i][j]和g[i][j]分别表示在[i,j]区间内，最后一个取数为i或j的方案数，实现了一个高效的算法。代码示例展示了如何通过状态转移方程进行计算，并考虑了模运算以防止整数溢出。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接

https://www.lydsy.com/JudgeOnline/problem.php?id=1996

题解

区间dp，记f[i][j]表示 $[i, j]$ 区间，最后一个取 $i$ 的方案数，g[i][j]表示 $[i, j]$ 区间最后一个取 $j$ 的方案数，状态转移很显然。

代码

#include <cstdio>

int read()
{
  int x=0,f=1;
  char ch=getchar();
  while((ch<'0')||(ch>'9'))
    {
      if(ch=='-')
        {
          f=-f;
        }
      ch=getchar();
    }
  while((ch>='0')&&(ch<='9'))
    {
      x=x*10+ch-'0';
      ch=getchar();
    }
  return x*f;
}

const int maxn=1000;
const int mod=19650827;

int n,h[maxn+10],f[maxn+10][maxn+10],g[maxn+10][maxn+10];

int main()
{
  n=read();
  for(int i=1; i<=n; ++i)
    {
      h[i]=read();
    }
  for(int i=1; i<=n; ++i)
    {
      f[i][i]=1;
    }
  for(int len=2; len<=n; ++len)
    {
      for(int i=1; i+len-1<=n; ++i)
        {
          int j=i+len-1;
          f[i][j]=(h[i]<h[i+1])*f[i+1][j]+(h[i]<h[j])*g[i+1][j];
          if(f[i][j]>=mod)
            {
              f[i][j]-=mod;
            }
          g[i][j]=(h[j]>h[j-1])*g[i][j-1]+(h[j]>h[i])*f[i][j-1];
          if(g[i][j]>=mod)
            {
              g[i][j]-=mod;
            }
        }
    }
  int ans=f[1][n]+g[1][n];
  if(ans>=mod)
    {
      ans-=mod;
    }
  printf("%d\n",ans);
  return 0;
}