#include <bits/stdc++.h>
using namespace std;
template<typename T>
struct MaxFlowMinCost {
struct Edge {
T cap, flow, cost;
int to;
};
vector<Edge> Edges;
vector<vector<int>> G;
int src, dest;
vector<int> Parent, ParentEdge, InQ;
vector<T> Dist;
MaxFlowMinCost& Initialize(int n, int m = 0) {
G.clear();
G.resize(n);
Edges.clear();
Edges.reserve(m);
Parent.resize(n);
ParentEdge.resize(n);
Dist.resize(n);
InQ.resize(n);
return *this;
}
void _addEdge(int from, int to, T cap, T cost) {
int ei = Edges.size();
Edges.push_back(Edge {cap, 0, cost, to});
G[from].push_back(ei);
}
MaxFlowMinCost& AddEdge(int from, int to, T cap, T cost) {
_addEdge(from, to, cap, cost);
_addEdge(to, from, 0, -cost);
return *this;
}
MaxFlowMinCost& SetSourceSink(int src, int dest) {
this->src = src; this->dest = dest;
return *this;
}
bool Bellman() {
static queue<int> Q;
fill(Dist.begin(), Dist.end(), numeric_limits<T>::max());
fill(Parent.begin(), Parent.end(), -1);
fill(InQ.begin(), InQ.end(), 0);
Dist[src] = 0;
Q.push(src);
while(!Q.empty()) {
int node = Q.front();
Q.pop();
InQ[node] = 0;
for(auto ei : G[node]) {
auto &e = Edges[ei];
if(e.flow < e.cap && Dist[e.to] > Dist[node] + e.cost) {
Dist[e.to] = Dist[node] + e.cost;
Parent[e.to] = node;
ParentEdge[e.to] = ei;
if(!InQ[e.to]) {
Q.push(e.to);
InQ[e.to] = 1;
}
}
}
}
return Parent[dest] != -1;
}
pair<T, T> Compute() {
T flow = 0, cost = 0;
while(Bellman()) {
T m = numeric_limits<T>::max();
for(int node = dest; node != src; node = Parent[node]) {
int ei = ParentEdge[node];
m = min(m, Edges[ei].cap - Edges[ei].flow);
}
for(int node = dest; node != src; node = Parent[node]) {
int ei = ParentEdge[node];
Edges[ei].flow += m;
Edges[ei ^ 1].flow -= m;
cost += Edges[ei].cost * m;
}
flow += m;
}
return {flow, cost};
}
};
int main() {
freopen("fmcm.in", "r", stdin);
freopen("fmcm.out", "w", stdout);
int n, m, s, d;
cin >> n >> m >> s >> d;
auto mfmc = MaxFlowMinCost<int>()
.Initialize(n, m)
.SetSourceSink(s - 1, d - 1);
while(m--) {
int a, b, c, k;
cin >> a >> b >> c >> k;
mfmc.AddEdge(a - 1, b - 1, c, k);
}
cout << mfmc.Compute().second << endl;
return 0;
}