[洛谷1250] 种树

---

题目描述

一条街的一边有几座房子。因为环保原因居民想要在路边种些树。路边的地区被分割成块，并被编号为1..n。每个块大小为一个单位尺寸并最多可种一棵树。每个居民想在门前种些树并指定了三个号码b,e,t。这三个数表示该居民想在b和e之间最少种t棵树。当然，b<=e,居民必须保证在指定地区不能种多于地区被分割成块数的树，即要求t<=e-b+1。允许居民想种树的各自区域可以交叉。出于资金短缺的原因，请你求出满足所有居民要求前提下，需要种树的最少数量。

---

输入格式

第一行为n，表示区域的个数；
第二行为m，表示房子的数目；
接下来m行描述居民的需要：b e t (0 < b <= e <= 30000,r <= e-b+1) 分别用一个空格分开。

---

输出格式

输出仅一个整数，即满足所有需要最少种树的数量。

---

样例数据

样例输入

9
4
1 4 2
4 6 2
8 9 2
3 5 2

样例输出

5

---

数据范围

30%的数据满足1<=n<=1000；1<=m<=500
100%的数据满足n<=30000；m<=5000

---

题目分析

贪心或差分约束
同[bsoj1271] 序列长度
不解释了

---

源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
const int maxn=50005;
struct Edge {
    int from,to,dist;
};
struct Difference_Constraints { //差分约束系统（最长路） 
    int n,m;
    vector<Edge>edges;
    vector<int>G[maxn];
    bool inque[maxn];
    int dist[maxn],used[maxn];
    int ans[maxn]; //ans表示不等式组答案 
    void init(int n) {
        this->n=n;
        edges.clear();
        for(int i=0; i<=n; i++)G[i].clear();
    }
    void AddEdge(int from,int to,int dist) {
        edges.push_back((Edge) {
            from,to,dist
        });
        m=edges.size();
        G[from].push_back(m-1);
    }
    void insert(int x,int y,int v,bool flag) { //x-y<=v (flag=0),x-y>=v (flag=1)
        v*=-1;
        if(flag==0)AddEdge(y,x,v);
        if(flag==1)AddEdge(x,y,-v);
    }
    bool spfa(int s) {
        for(int i=0; i<=n; i++)dist[i]=0x7fffffff/2;
        memset(inque,0,sizeof(inque));
        memset(used,0,sizeof(used));
        deque<int>Q;
        Q.push_back(s);
        dist[s]=0;
        inque[s]=1;
        used[s]++;
        while(!Q.empty()) {
            int Now=Q.front();
            Q.pop_front();
            inque[Now]=0;
            for(int i=0; i<G[Now].size(); i++) {
                Edge& e=edges[G[Now][i]];
                int Next=e.to;
                if(dist[Next]>dist[Now]+e.dist) {
                    dist[Next]=dist[Now]+e.dist;
                    if(!inque[Next]) {
                        used[Next]++;
                        if(used[Next]>n+1)return false; //负权回环
                        if(!Q.empty()&&dist[Next]<dist[Q.front()])Q.push_front(Next); //SLF优化
                        else Q.push_back(Next);
                        inque[Next]=1;
                    } 
                }
            }
        }
        return true;
    }
    bool main(int Start) { //求出正整数解 
        if(spfa(Start)==0)return false;
        for(int i=0; i<=n; i++)dist[i]*=-1;
        int Min=0x7fffffff/2;
        for(int i=0; i<=n; i++)Min=min(Min,dist[i]); //转为正数 
        for(int i=0; i<=n; i++)ans[i]=dist[i]-Min;
        return true;
    }
    void Output() {
        for(int i=0; i<=n; i++)printf("%d\n",ans[i]);
    }
};
int n,m;
Difference_Constraints dc;
int main() {
    n=Get_Int(); 
    m=Get_Int();
    dc.init(n+1);
    for(int i=1; i<=m; i++) {
        int x=Get_Int(),y=Get_Int(),v=Get_Int();
        dc.insert(y,x-1,v,0);
    }
    for(int i=1; i<=n; i++) {
        dc.insert(i,i-1,0,0);
        dc.insert(i-1,i,-1,0);
        dc.AddEdge(n+1,i,0);
    }
    dc.AddEdge(n+1,0,0);
    if(dc.main(n+1)==0) {
        puts("-1");
        return 0;
    }
    printf("%d\n",dc.ans[n]);
    return 0;
}

---