BZOJ 4798: [Ceoi2015]Calvinball championship

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
const int N = 10050;
const int p = 1000007;

inline void Add(int &x,LL y) { x=(x+y)%p; }

int n,mx,cur;
int a[N];
int f[2][N][2];

int main() {
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	cur=mx=f[1][1][1]=1;
	for(int i=2;i<=n;i++) {
		cur^=1;
		for(int j=1;j<=i;j++) {
			f[cur][j][0]=f[cur][j][1]=0;
			Add(f[cur][j][0],(LL)f[cur^1][j][0]*j%p+f[cur^1][j-1][0]);
			if(j==mx) Add(f[cur][j][0],f[cur^1][j][1]*(a[i]-1));
		}mx=max(mx,a[i]),f[cur][mx][1]=1;
//		for(int j=1;j<=n;j++) printf("%d,%d%c",f[cur][j][0],f[cur][j][1]," \n"[j==n]);
	}
	int ans=0;
	for(int i=1;i<=n;i++) Add(ans,f[cur][i][0]);
	printf("%d\n",(ans+1)%p);
	return 0;
}