Codeforces Round #326 (Div. 1) B Duff in Beach

B.Duff in Beach
题意：有一个长为l的序列B是由一个长度为n的序列A循环得到的，即$b_i = a_{i \ mod \ n}$，现在定义一个新的序列C，C的生成条件有三个：
1.C中每个元素都在B中，并且C中的每个元素都在不同的循环节里（即C中两个相邻的元素在B中的距离不会超过2n，不会低于n）；
2.C的长度超过k；
3.C不递减。
问C有多少种，答案模1e9+7。
思路：
先把B分成两组，一组长度是$\lfloor {l \over n} \rfloor * n$, 另一组就是$l\ mod \ n$。
对于第一组，因为C是不递减序列而且元素不在同一循环节里，所以可以对A排序，这样只要$a_i \leq a_j$，那么在下个循环节里$a_j$就可以接在$a_i$后面。
用dp[i][j]表示C的长度为i，以不超过j结尾时所有可能的C序列的数量。
在$i \leq k$时，对于dp[i][j] = (dp[i][j - 1] + dp[i - 1][j]);
在$i > k$时，ans的增量是dp[k][$max\{a_i\}$]。
对于第二组，要判断加上这部分长度会不会超过k，ans的增量在dp里找就可以了。
时间复杂度$O(nk)$。

/****************************************************
  >Created Date: 2015-10-01-13.34.56
  >My Soul, Your Beats!
****************************************************/
#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

const int INF = 1 << 30;
const long long  LINF = 1LL << 50;
const int M = 1e6 + 10;
const double PI = acos(-1.0);
const double eps = 1e-6;
int seq[M], seq2[M];
LL dp[2][M];
const int MOD = 1e9 + 7;
int main(){
    #ifndef ONLINE_JUDGE
  //      freopen("in.in", "r", stdin);
    #endif // ONLINE_JUDGE
    LL n, l, k;
    scanf("%I64d %I64d %I64d", &n, &l, &k);
    for(int i = 0; i < n; i++) {
        scanf("%d", &seq[i]);
        seq2[i] = seq[i];
    }
    sort(seq2, seq2 + n);
    for(int i = 0; i < n; i++) seq[i] = lower_bound(seq2, seq2 + n, seq[i]) - seq2;
    int maxn = 0;
    for(int i = 0; i < n; i++) maxn = max(maxn, seq[i]);
    maxn++;
    LL c = l / n, d = l % n;
    for(int i = 0; i < d; i++) seq2[i] = seq[i];
    LL ans = 0;
    int id = 0;
    for(int i = 0; i < min(c, k); i++){
        id ^= 1;
        for(int j = 0; j <= maxn; j++) dp[id][j] = 0;
        for(int j = 0; j < n; j++)
            dp[id][seq[j]] = (dp[id][seq[j]] + dp[id ^ 1][seq[j]] + 1) % MOD;
        for(int j = 1; j < maxn; j++)
            dp[id][j] = (dp[id][j] + dp[id][j - 1]) % MOD;
        ans = (ans + dp[id][maxn - 1]) % MOD;
    }
    if(c > k)
        ans = (ans + dp[id][maxn - 1] * ((c - k) % MOD)) % MOD;
    int t = id;
    if(c >= k) t ^= 1;
    for(int i = 0; i < d; i++){
        ans += dp[t][seq2[i]] + 1;
        ans %= MOD;
    }
    cout << ans << endl;
    return 0;
}