【PA2014】【BZOJ3711】Druzyny

Description

体育课上，n个小朋友排成一行（从1到n编号），老师想把他们分成若干组，每一组都包含编号连续的一段小朋友，每个小朋友属于且仅属于一个组。
第i个小朋友希望它所在的组的人数不多于d[i]，不少于c[i]，否则他就会不满意。
在所有小朋友都满意的前提下，求可以分成的组的数目的最大值，以及有多少种分组方案能达到最大值。

Input

第一行一个整数n(1<=n<=1000000)，表示小朋友的数目。
接下来n行，每行两个整数c[i],di，表示i所在组的人数的最小值和最大值。

Output

如果不存在这样的方案，仅输出一行NIE。
否则输出一行包含两个整数，组的数目的最大值、方案数量。（方案数量对1000000007取模）

Sample Input

样例输入1：

9

1 4

2 5

3 4

1 5

1 1

2 5

3 5

1 3

1 1

样例输入2：

2

1 1

2 2
Sample Output

样例输出1:

5 2

样例输出2：

NIE

HINT

Source

鸣谢Jcvb
这题有毒
智商不足不会做
↓这是Claris的课件

我连这个课件的代码实现都做不好..又去翻了Claris的博客..
(╯‵□′)╯︵┻━┻

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define P 1000000007#define GET (ch>='0'&&ch<='9')#define MAXN 1000010#define lchild rt<<1,l,mid#define rchild rt<<1|1,mid+1,r#define ln rt<<1#define rn rt<<1|1#define MAXINT 1000000using namespace std;int n;int c[MAXN],d[MAXN],g[MAXN];struct node{    int x,y;    node()  {   x=y=0;  }    node(int tx,int ty) {   x=tx;y=ty;  }    inline node operator + (const node& b)    {        if (x<b.x)  return  b;        if (x>b.x)  return  (node){x,y};        return (node){x,(y+b.y)%P};    }    inline node operator + (int tx) {   return (node){x+tx,y};  }    inline void operator +=(const node& b)  {   *this=*this+b;  }}f[MAXN];struct seg{    int l,r;    int pos,vald;    node f,flag;}tree[(MAXN<<1)+(int)(1e5)];inline void in(int &x){    char ch=getchar();x=0;    while (!GET)    ch=getchar();    while (GET) x=x*10+ch-'0',ch=getchar();}inline void build1(int rt=1,int l=0,int r=n){    tree[rt].l=l;tree[rt].r=r;    if (l==r)   {   tree[rt].pos=l;tree[rt].vald=d[l];return;   }    int mid=(l+r)>>1;build1(lchild);build1(rchild);    tree[rt].pos=c[tree[ln].pos]>c[tree[rn].pos]?tree[ln].pos:tree[rn].pos;    tree[rt].vald=min(tree[ln].vald,tree[rn].vald);}inline void build2(int rt=1,int l=0,int r=n){    tree[rt].flag=(node){-MAXINT,0};    if (l==r)   {   tree[rt].vald=d[l];tree[rt].f=f[l];return;  }    int mid=(l+r)>>1;build2(lchild);build2(rchild);    tree[rt].f=tree[ln].f+tree[rn].f;    tree[rt].vald=min(tree[ln].vald,tree[rn].vald);}inline void modify_add(int rt,int l,int r,node delta){    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    if (l<=L&&R<=r) {   tree[rt].flag+=delta;return;    }    if (r<=mid) modify_add(ln,l,r,delta);    else    if (l>mid)  modify_add(rn,l,r,delta);    else    modify_add(ln,l,mid,delta),modify_add(rn,mid+1,r,delta);}inline int query_c(int rt,int l,int r){    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    if (l<=L&&r>=R) return tree[rt].pos;    if (r<=mid) return query_c(ln,l,r);    if (l>mid)  return query_c(rn,l,r);    int t1=query_c(ln,l,mid),t2=query_c(rn,mid+1,r);    return c[t1]>c[t2]?t1:t2;}inline int query_d(int rt,int l,int r){    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    if (l<=L&&r>=R) return tree[rt].vald;    if (r<=mid) return query_d(ln,l,r);    if (l>mid)  return query_d(rn,l,r);    return min(query_d(ln,l,mid),query_d(rn,mid+1,r));}inline node query_f(int rt,int l,int r){    if (l>r)    return node(-MAXINT,0);    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    if (l<=L&&r>=R) return tree[rt].f;    if (r<=mid) return query_f(ln,l,r);    if (l>mid)  return query_f(rn,l,r);    return query_f(ln,l,mid)+query_f(rn,mid+1,r);}inline void modify_cov(int rt,int x,node delta){    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    if (L==R)   {   tree[rt].f=delta;return;    }    if (x<=mid) modify_cov(ln,x,delta); else    modify_cov(rn,x,delta);    tree[rt].f=tree[ln].f+tree[rn].f;}node ask(int x){    int rt=1,L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;    node t=node(-MAXINT,0);    while (L!=R)    {        t+=tree[rt].flag;        mid=(L+R)>>1;rt<<=1;        if (x>mid)  rt|=1;        L=tree[rt].l,R=tree[rt].r;    }    t+=tree[rt].flag;    return t;}inline void update(int l,int mid,int r){    int i=max(c[mid]+l,mid);    if (i>r||g[i]>=mid) return;    int newl=max(l,g[i]),newr=i-c[mid];    node tmp=query_f(1,newl,newr)+1;    for (;i<=mid-1+c[mid]&&i<=r;i++)    {        if (g[i]>newl)        {            if (g[i]>=mid)  return;            newl=g[i];tmp=query_f(1,newl,newr)+1;        }        f[i]+=tmp;newr++;        if (newr>=newl) tmp+=f[newr]+1;    }    while (i<=r)    {        if (g[i]>newl)        {            if (g[i]>=mid)  return;            newl=g[i];        }        tmp=query_f(1,newl,mid-1)+1;        int t=query_d(1,newl+1,n);        if (t>r)    {   modify_add(1,i,r,tmp);return;   }        modify_add(1,i,t-1,tmp);i=t;    }}void solve(int l,int r){    if (l==r)    {        if (l)  modify_cov(1,l,f[l]=f[l]+ask(l));        return;    }    int mid=query_c(1,l+1,r);    solve(l,mid-1);update(l,mid,r);solve(mid,r);}int main(){    in(n);int i=0,j=0;    for (i=1;i<=n;i++)  in(c[i]),in(d[i]);    build1();    for (i=0;i<=n;i++)  d[i]=n+1,f[i]=node(-MAXINT,0);    f[0]=node(0,1);    for (i=0;i<=n;i++)    {        while (j<i&&i-j>query_d(1,j+1,i))   j++;        g[i]=j;        if (d[g[i]]>n)  d[g[i]]=i;    }    build2();solve(0,n);printf(f[n].x>0?"%d %d\n":"NIE\n",f[n].x,f[n].y);}