Edu 69 E Culture Code (线段树建图+拓扑序求DAG最短路 / dp）

大致题意

给n个俄罗斯套娃娃，每个有个外体积out[i]，内体积 in[i]，外体积out[j]<=in[i]的情况下，j可以放在i里面，产生的空隙是in[i]-out[j]。现规定一个合格的套娃娃集合是在该套娃娃集合中，或者外面，无法再套上或者嵌套任何的剩余套娃。现要求产生最小空隙的套娃娃方式共有多少种。（注意一个没有套别的物品，其空隙就是in[i])

思路

这是一道我调了一下午的题…最后终于发现是最短路终点没赋初值…
codeforces好题 首先这个思路比较好想，每个套娃娃向in[i] >=out[j]的套娃娃连边，那么从各个起点出发，最短路即为最小空隙，顺便更新种树。但是这样建图会超时超内存。由于一个物品i，会向所有out[j]<=in[i]的物品连边，那么考虑按out递增排序。对于每个in[i]，二分查找到可连边的范围 [1,pos]，即区间 [1,pos] 都需要被连边，考虑建立一颗线段树，父亲都像儿子连权值为0的边，当l==r时，线段树叶子节点 x 向 id = l 连权值为 -out[l] 的边。在需要向 [1,pos] 连边的时候，将区间 [1,pos] 覆盖的区间节点找到，然后 i 向当前节点 x 连权值为 a[i].in 的边。
这样总体增加了nlogn个节点，和nlogn条边，但是建图只要 nlogn 的复杂度。
这样转化的局限和方法在于，要将加边的对象转化在一段连续的区间上，才能如此操作。在未排序之前是不行的。
然后DAG上求最短路，只要拓扑序更新可以线性求解最短路。

好像网上有博客可以 dp在O(n)或者O(nlogn)的时间范围内求解，蒟蒻以后再补。

代码

调的有点乱。

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000005
#define maxm 500006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define mod 1000000007
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int n,len,s,all,te;
struct pp{int in,out;
    bool operator < (const pp &rid)const{
        if(out!=rid.out)return out < rid.out;
        return in < rid.in;
    }
}a[maxn];
struct node{int v,w,nxt;}ed[maxn<<2];
int head[maxn],tot,du[maxn];
void add(int x,int y,int z){ed[++tot]={y,z,head[x]};head[x]=tot;du[y]++;}
queue<int>q;
ll dis[maxn],dp[maxn];
int calc(int x){
    int l=1,r=n,ans=-1;
    while(l<=r){
        int mid=(l+r)>>1;
        if(a[mid].out<=x){
            ans=mid;
            l=mid+1;
        }
        else r=mid-1;
    }
    return ans;
}
int t[maxn<<2];
void build(int x,int l,int r){
    all=max(all,x);
    if(l==r){add(n+x,l,-a[l].out);return ;}
    int mid=(l+r)>>1;
    int ls=x<<1,rs=x<<1|1;
    add(n+x,n+ls,0);add(n+x,n+rs,0);
    build(x<<1,l,mid);build(x<<1|1,mid+1,r);
}
void update(int x,int l,int r,int lp,int rp,int now){
    if(lp<=l&&r<=rp){add(now,n+x,a[now].in);return ;}
    int mid=(l+r)>>1;
    if(lp<=mid)update(x<<1,l,mid,lp,rp,now);
    if(mid+1<=rp)update(x<<1|1,mid+1,r,lp,rp,now);
}
int c[maxn];
int main()
{
    n=read();
    inc(i,1,n){a[i].out=read();a[i].in=read();}
    sort(a+1,a+n+1);
    build(1,1,n);
    s=n+all+1;te=s+1;
    inc(i,1,n){
        int pos=calc(a[i].in);
        if(pos>0){
            update(1,1,n,1,pos,i);
            c[1]++;c[pos+1]--;
        }
        else add(i,te,a[i].in);
    }
    inc(i,1,n){
        c[i]+=c[i-1];
        if(c[i]==0)add(s,i,0);
    }
    n+=all;
    inc(i,1,te){dp[i]=0;dis[i]=inf;}
    dp[s]=1;dis[s]=0;
    q.push(s);
    inc(i,1,n)if(du[i]==0)q.push(i);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=ed[i].nxt){
            int v=ed[i].v;
            du[v]--;
            if(dis[v]>dis[u]+ed[i].w){
                dis[v]=dis[u]+ed[i].w;
                dp[v]=dp[u];
            }
            else if(dis[v]==dis[u]+ed[i].w){
                dp[v]=(dp[v]+dp[u])%mod;
            }
            if(du[v]<=0)q.push(v);
        }
    }
    printf("%lld\n",dp[te]);
    return 0;
}