Even Three is Odd

题意：

问题是对于所有的长度为n，且$1<=ai<=n$的整数序列求 $\prod_{i=1}^{n-2}{max \{w_i,w_{i+1},w_{i+2}}\}$ 之和。

解法：

首先设dp状态为 $f(i,j,k)$ ，长度为$i+3$的，最大值为k，且最大值出现的位置集合为j的序列的乘积和。

显然可以由 $f(i-1,j2,k2)$ 转移到 $f(i,j,k)$，做前缀和优化，总效率$O(n^2 * 2^6)$

重新设计dp状态，改变j的定义，j表示最大值最后出现的位置。

这样对于状态 $f(i,j,k)$，我们确定了长度为$i+j$的序列的值，并且确定了$a(i+j+1)...a(i+2)<k$ 。

假设之前的三个数字最大值为$k2$，之后的最大值为k，这样的的话只要分为 $k>k2, k<k2, k=k2$ 讨论即可得出答案。

再加以前缀和优化，总效率$O(n^2)$。

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4
 5 #define LL long long
 6 #define N 2010
 7 #define P 1000000007LL
 8
 9 using namespace std;
10
11 int n;
12 LL w[N],S[N],S2[N],f[N][3][N];
13
14 LL sum(LL S[],int l,int r)
15 {
16     if(l>r) return 0LL;
17     LL ans = S[r]+P-S[l-1];
18     if(ans>=P) ans-=P;
19     return ans;
20 }
21
22 int main()
23 {
24 //    freopen("test.txt","r",stdin);
25     while(~scanf("%d",&n))
26     {
27         for(int i=1;i<=n;i++) scanf("%lld",&w[i]);
28         for(int i=0;i<=n-2;i++)
29             for(int k=1;k<=n;k++)
30                 f[i][0][k]=0, f[i][1][k]=0, f[i][2][k]=0;
31         for(int x1=1;x1<=n;x1++)
32             for(int x2=x1;x2<=n;x2++) f[0][2][x2]++;
33         for(int x1=1;x1<=n;x1++) f[0][1][x1]=1;
34         for(int i=1;i<=n-2;i++)
35         {
36             for(int k=1;k<=n;k++)
37             {
38                 S[k]  = S[k-1] +f[i-1][0][k];
39                 S2[k] = S2[k-1]+f[i-1][0][k]*(k-1)*(k-1);
40                 S2[k] += f[i-1][1][k]*(k-1);
41                 S2[k] += f[i-1][2][k];
42             }
43             for(int k=1;k<=n;k++)
44             {
45                 f[i][2][k] += sum(S2,1,k-1);
46                 f[i][1][k] += f[i-1][2][k];
47                 f[i][0][k] += f[i-1][1][k];
48                 f[i][2][k] += f[i-1][0][k]*(k-1)*(k-1);
49                 f[i][2][k] += f[i-1][1][k]*(k-1);
50                 f[i][2][k] += f[i-1][2][k];
51                 f[i][0][k] += sum(S,k+1,n);
52                 f[i][1][k] += sum(S,k+1,n)*k;
53                 f[i][2][k] += sum(S,k+1,n)*k*k;
54                 f[i][0][k] = f[i][0][k]%P * w[k]%P;
55                 f[i][1][k] = f[i][1][k]%P * w[k]%P;
56                 f[i][2][k] = f[i][2][k]%P * w[k]%P;
57             }
58         }
59         LL ans=0;
60         for(int k=1;k<=n;k++)
61         {
62             ans += f[n-2][0][k]*(k-1)*(k-1);
63             ans += f[n-2][1][k]*(k-1);
64             ans += f[n-2][2][k];
65             ans %= P;
66         }
67         cout << ans << endl;
68     }
69     return 0;
70 }