[APIO2016]划艇

题目描述

在首尔城中，汉江横贯东西。在汉江的北岸，从西向东星星点点地分布着 NNN 个划艇学校，编号依次为 111 到 NNN。每个学校都拥有若干艘划艇。同一所学校的所有划艇颜色相同，不同的学校的划艇颜色互不相同。颜色相同的划艇被认为是一样的。每个学校可以选择派出一些划艇参加节日的庆典，也可以选择不派出任何划艇参加。如果编号为 iii 的学校选择派出划艇参加庆典，那么，派出的划艇数量可以在 aia_ia​i​​ 至 bib_ib​i​​ 之间任意选择（ai≤bia_i leq b_ia​i​​≤b​i​​）。

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

输入所有学校的 ai,bia_i,b_ia​i​​,b​i​​ 的值，求出参加庆典的划艇有多少种可能的情况，必须有至少一艘划艇参加庆典。两种情况不同当且仅当有参加庆典的某种颜色的划艇数量不同。

输入输出格式

输入格式：

第一行包括一个整数 NNN，表示学校的数量。

接下来 NNN 行，每行包括两个正整数，用来描述一所学校。其中第 iii 行包括的两个正整数分别表示 ai,bia_i,b_ia​i​​,b​i​​（1≤ai≤bi≤1091 leq a_i leq b_i leq 10^91≤a​i​​≤b​i​​≤10​9​​）。

输出格式：

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

输入输出样例

输入样例#1：

2
1 2
2 3

输出样例#1：

7

说明

【样例解释】

在只有一所学校派出划艇的情况下有 444 种方案，两所学校都派出划艇的情况下有 333 种方案，所以答案为 777。

【数据范围】

子任务 111（999 分）：1≤N≤5001 leq N leq 5001≤N≤500 且对于所有的 1≤i≤N1 leq i leq N1≤i≤N，保证 ai=bia_i=b_ia​i​​=b​i​​。

子任务 222（222222 分）：1≤N≤5001 leq N leq 5001≤N≤500 且 ∑i=1N(bi−ai)≤106sum_{i=1}^N (b_i-a_i) leq 10^6∑​i=1​N​​(b​i​​−a​i​​)≤10​6​​。

子任务 333（272727 分）：1≤N≤1001 leq N leq 1001≤N≤100。

子任务 444（424242 分）：1≤N≤5001 leq N leq 5001≤N≤500。

http://m.blog.csdn.net/qq_22541499/article/details/51674707这个博客说的很清楚

组合数部分很难理解，建议手推一下

很玄学的是我用数组排序判重离散WA

改成STL set就AC了，回头还要改一下

已改完，为二分错误

以防万一，以后遇到一些奇怪的二分还是用upper_bound

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int Mod=1000000007;
typedef long long lol;
lol A[5001],f[501][2001],len[2001];
lol a[1001],b[1001],ge[2001],n,cnt;
int main()
{lol i,j,k;
  cin>>n;
  A[1]=1;
  for (i=2;i<=n;i++)
    {
      A[i]=(Mod-Mod/i)*A[Mod%i]%Mod;
    }
  for (i=1;i<=n;i++)
    {
      scanf("%lld%lld",&a[i],&b[i]);
      ge[i]=a[i];ge[n+i]=b[i]+1;
    }
  sort(ge+1,ge+2*n+1);
  for (i=1;i<2*n;i++)
    {
      if (ge[i]==ge[i+1])
       ge[i]=2e9;
    }
  sort(ge+1,ge+2*n+1);
  cnt=2*n;
  while (ge[cnt]==2e9) cnt--;
  for (i=1;i<=n;i++)
    {
      a[i]=upper_bound(ge,ge+cnt+1,a[i])-ge;
      b[i]=upper_bound(ge,ge+cnt+1,b[i])-ge;
      //cout<<a[i]<<' '<<b[i]<<endl;
    }
  for (i=0;i<=cnt;i++)
    {
      f[0][i]=1;
      if (i)
      len[i]=ge[i]-ge[i-1];
    }
  for (i=1;i<=n;i++)
    {
      f[i][0]=1;
      for (j=a[i];j<=b[i];j++)
    {
      f[i][j]=(f[i-1][j-1]*len[j])%Mod;
      lol c=len[j]-1,now=1;
      for (k=i-1;k;k--)
        {
          if (j>=a[k]&&j<=b[k])
        {
          now++;
          c=((c*(len[j]+now-2)%Mod)*A[now])%Mod;
          if (!c) break;
          f[i][j]+=(f[k-1][j-1]*c)%Mod;
          f[i][j]%=Mod;
        }
        }
    }
      for (j=1;j<=cnt;j++)
    {
      f[i][j]=(f[i][j]+f[i-1][j]+f[i][j-1]-f[i-1][j-1]+Mod)%Mod;
    }
    }
  cout<<(f[n][cnt]-1+Mod)%Mod;
}