APIO2016赛艇

• 首先每个学校的边界范围是\(1e9\)，肯定不能直接\(dp[i][j]\)表示前i所学校，第\(i\)所学校派出\(j\)艘船，但\(b<=500\)所有考虑把\(a,b\)离散，第二维改为当前这个学校派出数量在那个区间里
• 转移的时候，如果前面所有的学校都不在此区间内，那贡献为\(\sum_{i‘=1}^{i}\sum_{j‘=1}^{j}dp[i‘][j‘]\)
• 但是如果有学校要选在该区间内，问题便转化为在一段长度为len的区间内选出k个数，使他们递增，当然，有些数是可以不选的，该方案数为\(\sum_{i=0}^{k}C(k,k-i)*C(len,i)\),这个东西看起来一点都不好算,我们把这个式子化为\(C(len+k,k)\)
• 所以如果现在是第i个学校到第k个学校在此范围内，并且i,k学校强制选择，则方案数为\(C(i-k-1+len,i-k-1)\),贡献为\(\sum_{i‘=i-1}^{1}C(len+i-i‘-1,i-i‘-1)\sum_{k=1}^{i‘-1}\sum_{j‘=1}^{j-1}dp[k][j‘]\)其中&j&是当前枚举的区间
• 对于每一段区间，\(len\)相同，所以我们在枚举区间是处理一下组合数，后面的两重循环用前缀和处理一下，时间复杂度\(O(n^3)\)

  #include<bits/stdc++.h>
  using namespace std;
  typedef int sign;
  typedef long long ll;
  #define For(i,a,b) for(register sign i=(sign)a;i<=(sign)b;++i)
  #define Fordown(i,a,b) for(register sign i=(sign)a;i>=(sign)b;--i)
  const int N=500+5;
  bool cmax(sign &a,sign b){return (a<b)?a=b,1:0;}
  bool cmin(sign &a,sign b){return (a>b)?a=b,1:0;}
  template<typename T>T read()
  {
    T ans=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)&&ch!=‘-‘)ch=getchar();
    if(ch==‘-‘)f=-1,ch=getchar();
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch-‘0‘),ch=getchar();
    return ans*f;
  }
  template<typename T>void write(T x,char y)
  {
    if(x==0)
    {
    putchar(‘0‘);putchar(y);
    return;
    }
    if(x<0)
    {
    putchar(‘-‘);
    x=-x;
    }
    static char wr[20];
    int top=0;
    for(;x;x/=10)wr[++top]=x%10+‘0‘;
    while(top)putchar(wr[top--]);
    putchar(y);
  }
  void file()
  {
    #ifndef ONLINE_JUDGE
    freopen("3643.in","r",stdin);
    freopen("3643.out","w",stdout);
    #endif
  }
  int n,a[N],b[N],l[N<<1],tot;
  void input()
  {
  n=read<int>();
  For(i,1,n)
  {
      a[i]=read<int>(),b[i]=read<int>();
      l[++tot]=a[i],l[++tot]=b[i]+1;
  }
  }
  void init()
  {
  sort(l+1,l+tot+1);
  tot=unique(l+1,l+tot+1)-l-1;
  For(i,1,n)
  {
      a[i]=lower_bound(l+1,l+tot+1,a[i])-l;
      b[i]=lower_bound(l+1,l+tot+1,b[i]+1)-l;
  }
  }
  const int mo=1e9+7;
  int inv[N],dp[N],C[N];
  void work()
  {
  int len;
  inv[1]=1;For(i,2,n)inv[i]=1ll*(mo-mo/i)*inv[mo%i]%mo;
  C[0]=dp[0]=1;
  For(i,1,tot-1)
  {
      len=l[i+1]-l[i];
      For(i,1,n)C[i]=1ll*C[i-1]*(len+i-1)%mo*inv[i]%mo;
      Fordown(j,n,1)if(a[j]<=i&&i+1<=b[j])
      {
          int f=0,pos=1,c=len;
          Fordown(k,j-1,0)
          {
              (f+=1ll*dp[k]*c%mo)%=mo;
              if(a[k]<=i&&i+1<=b[k])c=C[++pos];
          }
          (dp[j]+=f)%=mo;
      }
  }
  int ans=0;
  For(i,1,n)(ans+=dp[i])%=mo;
  write(ans,‘\n‘);
  }
  int main()
  {
  file();
  input();
  init();
  work();
  return 0;
  }

原文地址：https://www.cnblogs.com/dengyixuan/p/8881136.html