// #define _GLIBCXX_DEBUG // for STL debug (optional)
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <cstring>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <utility>
#include <algorithm>
#include <map>
#include <set>
#include <complex>
#include <cmath>
#include <limits>
#include <cfloat>
#include <climits>
#include <ctime>
#include <cassert>
#include <numeric>
#include <fstream>
#include <functional>
#include <bitset>
using namespace std;
using ll = long long int;
using int64 = long long int;
 
template<typename T> void chmax(T &a, T b) {a = max(a, b);}
template<typename T> void chmin(T &a, T b) {a = min(a, b);}
template<typename T> void chadd(T &a, T b) {a = a + b;}
 
int dx[] = {0, 0, 1, -1};
int dy[] = {1, -1, 0, 0};
const int INF = 1LL << 29;
const ll LONGINF = 1LL << 60;
const ll MOD = 1000000007LL;

// Ford-Fulkerson 法による 最大流 O( F |E| )
// Bellman-Ford 法による 最小費用流 O( F |V| |E| )
// [条件に注意] Dijkstra 法による 最小費用流 O( F |E| log |V| )

template <typename CapTp=int, typename CostTp=int>
struct Edge {
    int to, rev;
    CapTp cap; CostTp cost;
    bool is_rev;
    Edge(int t, bool f, int r, CapTp ca, CostTp co=0)
        : to(t), rev(r), cap(ca), cost(co), is_rev(f) {}
};

template <typename CapTp=int, typename CostTp=int>
struct Flow {
    using Graph = vector< vector< Edge<CapTp, CostTp> > >;
    Graph G; const CapTp IA; const CostTp IO;
    vector< pair<int, int> > r_edges;
    Flow(int N_, CapTp IA_=1<<29, CostTp IO_=1<<29)
        : G(N_), IA(IA_), IO(IO_), r_edges() {}
    // 辺を追加 (from -> to に流量 ca, コスト co)
    void add_edge(int from, int to, CapTp ca, CostTp co=0) {
        G[from].emplace_back(to, false, G[to].size(), ca, co);
        G[to].emplace_back(from, true, G[from].size() - 1, 0, -co);
        r_edges.emplace_back(to, G[to].size() - 1);
    }
    // k 番目の辺にいくつ流れたか
    CapTp get_flowed_cap(size_t k) {
        if(r_edges.size() <= k) return -1;
        int v, i; tie(v, i) = r_edges[k];
        return G[v][i].cap;
    }
    // s -> t 最大流
    CapTp max_flow(int s, int t) {
        vector<bool> used(G.size());
        auto dfs = [&](auto &&func, int v, int t, CapTp f) -> CapTp {
            if(v == t) return f;
            used[v] = true;
            for(auto &e : G[v]) {
                if(used[e.to] or e.cap == 0) continue;
                CapTp d = func(func, e.to, t, min(f, e.cap));
                if(d == 0) continue;
                e.cap -= d; G[e.to][e.rev].cap += d;
                return d;
            }
            return 0;
        };

        CapTp res(0);
        while(true) {
            fill(used.begin(), used.end(), false);
            CapTp delta = dfs(dfs, s, t, IA);
            if(delta == 0) return res;
            res += delta;
        }
    }
    // ベルマンフォードをつかって最小費用流
    pair<CostTp, bool> mincost_flow(int s, int t, CapTp f) {
        vector<CostTp> dist(G.size()); CostTp res(0);
        vector<int> prevv(G.size()), preve(G.size());
        while(f > 0) {
            fill(dist.begin(), dist.end(), IO);
            dist[s] = 0;
            while(1) {
                bool upd = false;
                for(int v=0; v<(int)G.size(); v++) {
                    if(dist[v] == IO) continue;
                    for(size_t i=0; i<G[v].size(); i++) {
                        auto &e = G[v][i];
                        if(e.cap == 0 or dist[e.to] <= dist[v] + e.cost) continue;
                        dist[e.to] = dist[v] + e.cost;
                        prevv[e.to] = v, preve[e.to] = i;
                        upd = true;
                    }
                }
                if(!upd) break;
            }

            if(dist[t] == IO) return make_pair(res, false);
            CapTp d = f;
            for(int v=t; v!=s; v=prevv[v]) d = min(d, G[prevv[v]][preve[v]].cap);
            f -= d; res += d * dist[t];
            for(int v=t; v!=s; v=prevv[v]) {
                auto &e = G[prevv[v]][preve[v]];
                e.cap -= d, G[v][e.rev].cap += d;
            }
        }
        return make_pair(res, true);
    }
    // ポテンシャルの導入により、ダイクストラ法で最小費用流を解く
    // [仮定している条件]
    //     1. グラフに負の閉路が存在しない (流量の 0 初期化のため)
    //        もし存在するならベルマンフォードで負の閉路を見つけ
    //        そこに流せるだけ流してスタート
    //     2. グラフに負の辺が存在しない (pot_0 の計算可能性)
    //        もし存在する場合は最初のみベルマンフォードを使う必要あり
    pair<CostTp, bool> fast_mincost_flow(int s, int t, CapTp f) {
        CostTp res = 0;
        vector<CostTp> dist(G.size()), pot(G.size());
        vector<int> prevv(G.size()), preve(G.size());
        while(f > 0) {
            using PT = pair<CostTp, int>;
            priority_queue< PT, vector<PT>, greater<PT> > que;
            fill(dist.begin(), dist.end(), IO);

            dist[s] = 0;
            que.push(make_pair(0, s));
            while(!que.empty()) {
                PT cur = que.top(); que.pop();
                int v = cur.second;
                if(dist[v] < cur.first) continue;
                for(size_t i=0; i<G[v].size(); i++) {
                    auto& e = G[v][i];
                    if(e.cap > 0 and dist[e.to] > dist[v] + e.cost + pot[v] - pot[e.to]) {
                        dist[e.to] = dist[v] + e.cost + pot[v] - pot[e.to];
                        prevv[e.to] = v;
                        preve[e.to] = i;
                        que.push(make_pair(dist[e.to], e.to));
                    }
                }
            }
            if(dist[t] == IO) {
                return make_pair(res, false);
            }
            for(int v=0; v<(int)G.size(); v++) pot[v] += dist[v];

            CapTp d = f;
            for(int v=t; v!=s; v=prevv[v]) {
                d = min(d, G[prevv[v]][preve[v]].cap);
            }
            f -= d;
            res += d * pot[t];
            for(int v=t; v!=s; v=prevv[v]) {
                auto& e = G[prevv[v]][preve[v]];
                e.cap -= d;
                G[v][e.rev].cap += d;
            }
        }
        return make_pair(res, true);
    }    
};

template <typename cap_tp=int, typename cost_tp=int,
          template<typename, typename> class flow_fw=Flow>
struct GeneralizedMincostFlow {
private:
    flow_fw<cap_tp, cost_tp> fl;
    vector<cap_tp> D;
    int N, source, sink;
    cap_tp Z_cap;
    cost_tp Z_cost, ofs;

public:
    GeneralizedMincostFlow(int N_,
                           cap_tp zero_cap=0, cost_tp zero_cost=0)
        : fl(N_+2), D(N_), N(N_), source(N_), sink(N_+1),
          Z_cap(zero_cap), Z_cost(zero_cost), ofs(zero_cost) {}

    void add_edge(int u, int v, cap_tp low, cap_tp high, cost_tp cost) {
        assert(low <= high);
        // コストが負の場合
        if(cost < Z_cost) {
            // 流すだけ得なので high 流れたことにする
            ofs += high * cost;
            D[u] -= high; D[v] += high;
            // 目一杯流したのを high - low 分だけキャンセル可能
            if(high - low > Z_cap) {
                fl.add_edge(v, u, high - low, -cost);
            }
        }
        else {
            // v 側に low 流れたことにする
            ofs += low * cost;
            D[u] -= low; D[v] += low;
            // high - low 分の辺を張る
            if(high - low > Z_cap) {
                fl.add_edge(u, v, high - low, cost);
            }
        }
    }
    void add_edge(int u, int v, cap_tp cap, cost_tp cost) {
        add_edge(u, v, Z_cap, cap, cost);
    }
    void change_d(int u, cap_tp delta) {
        D[u] += delta;
    }
    
    pair<cost_tp, bool> mincost_flow() {
        cap_tp sum_d = 0;
        for(int i=0; i<N; i++) {
            if(D[i] > 0) {
                fl.add_edge(source, i, D[i], 0);
                sum_d += D[i];
            }
            else {
                fl.add_edge(i, sink, -D[i], 0);
            }
        }

        pair<cost_tp, bool> res = fl.mincost_flow(source, sink, sum_d);
        res.first += ofs;
        return res;
    }
};

int main() {
    int N, M, s, t; scanf("%d%d%d%d", &N, &M, &s, &t);
    s--; t--;
    GeneralizedMincostFlow<ll, ll> fl(N);
    vector< vector<int> > adj(N, vector<int>(N, INF));
    for(int i=0; i<N; i++) adj[i][i] = 0;
    for(int i=0; i<M; i++) {
        int u, v, d, c; scanf("%d%d%d%d", &u, &v, &d, &c);
        u--; v--;
        adj[u][v] = d;
        fl.add_edge(u, v, c, d);
    }
    for(int k=0; k<N; k++) {
        for(int i=0; i<N; i++) {
            for(int j=0; j<N; j++) {
                chmin(adj[i][j], adj[i][k] + adj[k][j]);
            }
        }
    }
    fl.add_edge(t, s, INF, -adj[s][t] - 1);

    ll res, ok; tie(res, ok) = fl.mincost_flow();
    assert(ok);
    printf("%lld\n", -res);
    return 0;
}



Measured memory_bytes: 3396