Problem:
A: Attack the Moles

ICPC で良い成績を収めるには修行が欠かせない．うさぎは ICPC で勝ちたいので，今日も修行をすることにした．

今日の修行は，もぐらたたきを何回も行って，反射神経と記憶力を高めようというものである．出てくるもぐらを次々に叩き，出来るだけ多くのポイントを獲得したい．

もぐらが出てくる可能性のある場所は直線状に並んでいて，基準点からの距離によって座標が定まっている．うさぎはしばらく修行を続けるうちに，もぐらの出現する場所と時間が常に一緒であることに気が付いた．うさぎは，その情報をすべて記録し，コンピュータで解析を行うことにした．

もぐらを叩くには，もぐらの出現位置に手を動かした後，もぐらの出てくるタイミングにぴったり合わせてもぐらを叩かなければならない．もぐらをうまく叩けると，そのもぐらに応じてポイントを得ることが出来る．もぐらを叩く動作は一瞬で行うことが出来るが，手を移動させる速さには限界がある．うさぎはもぐらを叩くにあたって左右両方の手を用いることができる．左手と右手は独立に動かすことが可能であるが，左手は常に右手より座標が小さい位置に存在しなければならない．このような条件下で，最大でどれだけのポイントが得られるかを調べたい．

Input

N V XLeft XRight
X1 T1 P1
...
XN TN PN

N は出てくるもぐらの数，V は手を移動させられる最大の速さ，XLeft, XRight はそれぞれ，左手，右手の初期位置の座標である．Xi, Ti, Pi はそれぞれ，i 番目のもぐらの，出現する位置の座標，ゲーム開始から出現までの時間，叩けた際に得られるポイントである．

1 ≤ N ≤ 3,000，1 ≤ V ≤ 10,000，1 ≤ XLeft < XRight ≤ 100,000，1 ≤ X1 ≤ X2 ≤ ... ≤ XN ≤ 100,000，1 ≤ Ti ≤ 100,000，1 ≤ Pi ≤ 100,000 を満たす．(Xi, Ti) として同一の組は複数回現れない．

Output

うさぎが得られる最大のポイントを 1 行に出力せよ．

Sample Input 1

3 10 150 250
100 20 123
201 10 67
202 10 45

Sample Output 1

190

Sample Input 2

1 7 20 90
55 5 73

Sample Output 2

73

Sample Input 3

10 2 1000 2000
400 300 1
600 200 1
700 800 1
700 500 1
900 600 1
1000 700 1
1300 900 1
1400 400 1
1500 1000 1
2000 100 1

Sample Output 3

10

Solution in C++14 (judge verdict: Accepted)

#include <bits/stdc++.h>
#define ll long long
#define INF 1000000005
#define MOD 1000000007
#define EPS 1e-10
#define rep(i,n) for(int i=0;i<(int)(n);++i)
#define rrep(i,n) for(int i=(int)(n)-1;i>=0;--i)
#define srep(i,s,t) for(int i=(int)(s);i<(int)(t);++i)
#define each(a,b) for(auto (a): (b))
#define all(v) (v).begin(),(v).end()
#define len(v) (int)(v).size()
#define zip(v) sort(all(v)),v.erase(unique(all(v)),v.end())
#define cmx(x,y) x=max(x,y)
#define cmn(x,y) x=min(x,y)
#define fi first
#define se second
#define pb push_back
#define show(x) cout<<#x<<" = "<<(x)<<endl
#define spair(p) cout<<#p<<": "<<p.fi<<" "<<p.se<<endl
#define svec(v) cout<<#v<<":";rep(kbrni,v.size())cout<<" "<<v[kbrni];cout<<endl
#define sset(s) cout<<#s<<":";each(kbrni,s)cout<<" "<<kbrni;cout<<endl
#define smap(m) cout<<#m<<":";each(kbrni,m)cout<<" {"<<kbrni.first<<":"<<kbrni.second<<"}";cout<<endl

using namespace std;

typedef pair<int,int> P;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<double> vd;
typedef vector<P> vp;
typedef vector<string> vs;

const int MAX_N = 100005;

struct Edge{
    Edge(const short to,const short reverse,const char capacity,const int cost,bool is_reverse = false):to{to},reverse{reverse},capacity{capacity},flow{0},cost{cost},is_reverse{is_reverse}{}
    short to;
    short reverse;
    char capacity;
    char flow;
    int cost;
    bool is_reverse;
};

class tsort {
public:
	vector<int> deg;
    vector<short> order;
	short V;
	tsort(const vector<vector<Edge>>& g)
    {
		V = len(g);
		deg.resize(V,0);
        rep(i,len(g)){
            for(const Edge& e : g[i]){
                if(!e.is_reverse){
                    deg[e.to]++;
                }
            }
        }
        make(g);
	}
	void make(const vector<vector<Edge>>& G) {
		queue<short> que;
		rep(i,V){
	        if(deg[i] == 0){
	            que.push(i);
	        }
	    }
		while(!que.empty()){
			short p = que.front();
	        que.pop();
			order.push_back(p);
			rep(i,G[p].size()){
	            if(--deg[G[p][i].to] == 0){
	                que.push(G[p][i].to);
	            }
	        }
		}
	}
};

class CostFlow
{
public:
    using T = int;

    CostFlow(const short v):V{v}
    {
        edge.resize(v);
        m_dist.resize(v);
        m_potential.resize(v);
        m_prev_v.resize(v);
        m_prev_e.resize(v);
    }

    void addEdge(const short from, const short to, const char capacity,const T cost)
    {
        edge[from].pb(Edge{to,(short)edge[to].size(),capacity,cost,false});
        edge[to].pb(Edge{from,(short)(edge[from].size()-1),0,-cost,true});
    }

    T minCostFlow(const short s,const short t, char f)
    {
    //    const T INF = numeric_limits<T>::max()/100;
        using P = pair<T,short>;
        T res = 0;
        fill(all(m_potential),0);
        tsort ts(edge);
// potential
        for(short i =0; i < V; i++){
            m_potential[i] = INF;
        }
        m_potential[s] = 0;
//        assert(s == order[0]);
    //    svec(order);

        for(const short i:ts.order){
            if(m_potential[i] != INF){
                for(const auto& e: edge[i]){
                    if(e.is_reverse){
                        continue;
                    }
                    m_potential[e.to] = min(m_potential[e.to],m_potential[i]+e.cost);
                }
            }
        }

//        svec(m_potential);

        while(f > 0)
        {
            priority_queue<P,vector<P>,greater<P>> q;
            fill(all(m_dist),INF);
            m_dist[s]=0;
            q.push(make_pair(0,s));
            while(not q.empty()){
                const P p = q.top();
                q.pop();
                const short v = p.second;
                if(m_dist[v] < p.first){
                    continue;
                }
                for(short i = 0; i < edge[v].size();i++){
                    const auto& e = edge[v][i];
                    if(e.capacity > e.flow and m_dist[e.to] > m_dist[v]+e.cost+m_potential[v]-m_potential[e.to]){
                        m_dist[e.to] = m_dist[v]+e.cost+m_potential[v]-m_potential[e.to];
                        m_prev_v[e.to] = v;
                        m_prev_e[e.to] = i;
                        q.push(make_pair(m_dist[e.to],e.to));
                    }
                }
            }

            if(m_dist[t] == INF){
                return -1;
            }
            for(short v = 0; v < V; v++){
                m_potential[v] += m_dist[v];
            }
            char d = f;
            for(short v = t; v != s; v = m_prev_v[v]){
                const auto& e = edge[m_prev_v[v]][m_prev_e[v]];
                d = min(d,(char)(e.capacity-e.flow));
            }
            f-= d;
            res += d*m_potential[t];
            for(short v = t; v != s;v=m_prev_v[v]){
                auto& e = edge[m_prev_v[v]][m_prev_e[v]];
                e.flow += d;
                edge[v][e.reverse].flow-=d;
            }
        }
        return res;
    }

    const short V;
    vector<vector<Edge>> edge;
private:
    vector<T> m_dist;
    vector<T> m_potential;
    vector<short> m_prev_v;
    vector<short> m_prev_e;
};

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);
    int n,v,xr,xl;
    cin >> n >> v >> xr >> xl;
    vvl vec(n,vl(3));
    rep(i,n){
        rep(j,3){
            cin >> vec[i][j];
        }
    }
    CostFlow mc(2*n+4);
    rep(i,n){
        mc.addEdge(i,n+i,1,0);
    }
    rep(i,2){
        mc.addEdge(2*n,2*n+1+i,1,0);
    }
    rep(i,n){
        if(abs(xl-vec[i][0]) <= v*vec[i][1]){
            mc.addEdge(2*n+1,i,1,0);
        }
    }
    rep(i,n){
        if(abs(xr-vec[i][0]) <= v*vec[i][1]){
            mc.addEdge(2*n+2,i,1,0);
        }
    }
    rep(i,n){
        rep(j,n){
            if(i != j){
                if(abs(vec[i][0]-vec[j][0]) <= v*(vec[j][1]-vec[i][1])){
                    mc.addEdge(n+i,j,1,-vec[i][2]);
                }
            }
        }
    }
    rep(i,n){
        mc.addEdge(n+i,2*n+3,1,-vec[i][2]);
    }
    rep(i,2){
        mc.addEdge(2*n+1+i,2*n+3,1,0);
    }
    cout << -mc.minCostFlow(2*n,2*n+3,2) << "\n";
    return 0;
}