dp[i]:以第i个字符为max,左边的可以构成的上升子序列的个数 。
dp2[i];以第i个字符为max,右边的可以构成的下降子序列的个数 。
num[i][j];在第i个位置以字符j为结尾的上升子序列的个数
num2[i][j];在第i个位置以字符j为开头的下降子序列的个数
很容易的推出状态转移:
int x=str[i-1]-'a';
for(j=0;j<26;j++)num[i][j]=num[i-1][j]%mod;
for(j=0;j<x;j++)num[i][x]+=num[i-1][j];
num[i][x]++;
num[i][x]%=mod;
dp[i]=0;
for(j=0;j<x;j++)dp[i]+=num[i-1][j];
dp[i]%=mod;
以此类推dp2[i],num2[i][j];
则ans=sigma((dp[i])*(dp2[i]));(1<=i<=n)
很显然,这道题目并不是很难。但是一定要注意,取余的时候对每一步都要取余。
不然很容易造成dp数组和num数组的越界。
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<stack>
using namespace std;
#define LL long long
#define maxn 110000
#define mod 2012
char str[maxn];
LL dp[maxn];
LL dp2[maxn];
LL num[maxn][27];
LL num2[maxn][27];
int main()
{// freopen("data1.in","r",stdin);int i,j,n;while(~scanf("%d",&n)){scanf("%s",str);memset(dp,0,sizeof(dp));memset(dp2,0,sizeof(dp2));memset(num,0,sizeof(num));memset(num2,0,sizeof(num2));for(i=1;i<=n;i++){int x=str[i-1]-'a';for(j=0;j<26;j++)num[i][j]=num[i-1][j]%mod;for(j=0;j<x;j++)num[i][x]+=num[i-1][j];num[i][x]++;num[i][x]%=mod;dp[i]=0;for(j=0;j<x;j++)dp[i]+=num[i-1][j];dp[i]%=mod;}for(i=n;i>=1;i--){int x=str[i-1]-'a';for(j=0;j<26;j++)num2[i][j]=num2[i+1][j]%mod;for(j=0;j<x;j++)num2[i][x]+=num2[i+1][j];num2[i][x]++;num2[i][x]%=mod;dp2[i]=0;for(j=0;j<x;j++)dp2[i]+=num2[i+1][j];dp[i]%=mod;}// cout<<endl;LL maxx=0;for(i=1;i<=n;i++){maxx+=(dp[i])*(dp2[i]);// maxx=maxx%mod;}cout<<maxx%mod<<endl;}return 0;
}