BZOJ 4584 luogu P3643: [Apio2016]赛艇

4584: [Apio2016]赛艇

Time Limit: 70 Sec  Memory Limit: 256 MB
[Submit][Status][Discuss]

Description

在首尔城中，汉江横贯东西。在汉江的北岸，从西向东星星点点地分布着个划艇学校，编号依次为到。每个学校都

拥有若干艘划艇。同一所学校的所有划艇颜色相同，不同的学校的划艇颜色互不相同。颜色相同的划艇被认为是一

样的。每个学校可以选择派出一些划艇参加节日的庆典，也可以选择不派出任何划艇参加。如果编号为的学校选择

派出划艇参加庆典，那么，派出的划艇数量可以在Ai至Bi之间任意选择（Ai<=Bi）。值得注意的是，编号为i的学

校如果选择派出划艇参加庆典，那么它派出的划艇数量必须大于任意一所编号小于它的学校派出的划艇数量。输入

所有学校的Ai、Bi的值，求出参加庆典的划艇有多少种可能的情况，必须有至少一艘划艇参加庆典。两种情况不同

当且仅当有参加庆典的某种颜色的划艇数量不同

 

Input

 第一行包括一个整数N，表示学校的数量。接下来N行，每行包括两个正整数，用来描述一所学校。其中第行包括的

两个正整数分别表示Ai,Bi（1<=Ai<=Bi<=10^9），N<=500

Output

 输出一行，一个整数，表示所有可能的派出划艇的方案数除以1,000,000,007得到的余数

Sample Input

2

1 2

2 3

Sample Output

7

HINT

 

Source

 

BZOJ链接

UOJ链接

看dalao博客弄懂的 链接

#include<bits/stdc++.h>

using namespace std;

const long long INF = 1000000007;

const int MAXN = 510;
long long n, p, top, ans;
long long inv[MAXN];
long long a[MAXN], b[MAXN], l[MAXN], r[MAXN]; 
long long h[MAXN + MAXN];
long long C[MAXN], g[MAXN];
map<long long, int> k;

template <typename tn> void read (tn & a) {
    tn x = 0, f = 1; char c = getchar();
    while (c < '0' || c > '9') { if (c == '-') f = -1; c = getchar(); }
    while (c >= '0' && c <= '9') { x = x * 10 + c - '0'; c = getchar(); }
    a = f == 1 ? x : -x;
}

void get_inv() {
    inv[1] = 1;
    for (int i = 2; i <= n; ++i) {
        inv[i] = - ( inv[INF % i] * (INF / i) ) % INF;
        while (inv[i] < 0) inv[i] += INF;
    }
}

int main() {
    read(n);
    get_inv();
    for (int i = 1; i <= n; ++i) {
        read(a[i]);
        read(b[i]);
        ++b[i];
        h[i + i - 1] = a[i];
        h[i + i] = b[i];
    }
    sort(h + 1, h + 1 + n + n);
    top = unique(h + 1, h + n + n + 1) - h - 1;
    for (int i = 1; i <= top; ++i) {
        k[h[i]] = i;
    }
    for (int i = 1; i <= n; ++i) {
        l[i] = k[a[i]];
        r[i] = k[b[i]];
    }
//    cout << top << "\n";
//    for (int i = 1; i <= top; ++i) cout << h[i]<<" ";
//    for (int i = 1; i <= n; ++i) cout << l[i]<<" "<<r[i]<<"\n";
    g[0] = 1;
    for (int j = 1; j <= top; ++j) {
        int L = h[j + 1] - h[j];
        C[0] = 1;
//        cerr<<"j = "<<j<<"\n";
        for (int i = 1; i <= n; ++i) {
            C[i] = C[i - 1] * (L + i - 1) % INF * inv[i] % INF;
//            cerr<<i<<" "<<C[i]<<"\n";
        }
        for (int i = n; i > 0; --i) {
            if (l[i] <= j && r[i] > j) {
                long long f = 0, m = 1, c = L;
                for (int p = i - 1; p >= 0; --p) {
                    f = (f + c * g[p]) % INF;
                    if (l[p] <= j && r[p] >= j + 1) {
                        c = C[++m];
                    }
                }
                g[i] = (g[i] + f) % INF;
            }
        }
    }
    for (int i = 1; i <= n; ++i) ans = (ans + g[i]) % INF;
    printf("%lld\n", ans);
    return 0;
}