hdu5760

区间Dp，经典问题最长上升公共子序列的变式；
计数需要去重，思想在于对于本质相同的序列均将其表示到最大的哪个；
然后dp时依据包含关系就可以去重啦；

#include<cstdio>
#include<algorithm>
#include<cstring>
#define rep(i,k,n) for(int i=k;i<=n;i++)
#define rep2(i,k,n) for(int i=k;i>=n;i--)
#define cnt second
#define len first
#define mk make_pair
#define pa pair<int,int>
using namespace std;
const int N=5e3+7;
const int mod=1e9+7;
void upd(int& x,int y){x+=y;if(x>=mod)x-=mod;}
int n,m,a[N],p[N],pre[N][N],nxt[N][N];
pa d[N][N],now;
void init(){
    rep(i,1,n)scanf("%d",&a[i]),p[i]=a[i];
    sort(p+1,p+n+1);m=unique(p+1,p+n+1)-p;
    rep(i,1,n)a[i]=lower_bound(p+1,p+m,a[i])-p;
    memset(p,0,sizeof(p));rep(i,0,n+1){rep(j,1,m)pre[i][j]=p[j];p[a[i]]=i;}
    memset(p,0,sizeof(p));rep2(i,n+1,0){rep(j,1,m)nxt[i][j]=p[j];p[a[i]]=i;}
}
int main(){
    while(scanf("%d",&n)==1){
        init();
        rep2(i,n,1){
            d[i][i]=mk(1,1);
            now=mk(0,1);
            rep(j,i+1,n){
                d[i][j]=mk(0,0);
                if(a[i]==a[j])d[i][j]=mk(now.len+2,now.cnt);
                if(a[i]>=a[j]){
                int l=nxt[i][a[j]],r=pre[j][a[j]];
                    if(d[l][j].len>now.len)now=d[l][j];
                    else if(d[l][j].len==now.len){
                        if(l<=r && d[l][r].len==d[l][j].len){
                            upd(now.cnt,(mod-d[l][r].cnt)%mod);
                        }
                        upd(now.cnt,d[l][j].cnt);
                    }
                }
            }
        }now=mk(0,0);
        rep(i,1,m){
            int l=nxt[0][i],r=pre[n+1][i];
            if(!l || !r)break;
            if(d[l][r].len>now.len)now=d[l][r];
            else if(d[l][r].len==now.len){
                upd(now.cnt,d[l][r].cnt);
            }
        }printf("%d %d\n",now.len,now.cnt);
    }
}