#include <cstdio>
#include <queue>

using namespace std;

struct edge
{
    int nod,cost;
    bool operator <(const edge aux) const
    {
        return cost>aux.cost;
    }
};

const int inf=1e9;
vector<edge> v[360];
priority_queue<edge> q;
queue<int> q1;
int n,sursa,dest,dist[360],distt[360],cap_flow[360][360],d[360],vaz[360],tata[360];

void bellman_ford()
{
    for(int i=1;i<=n;i++) dist[i]=inf;
    dist[sursa]=0;
    vaz[sursa]=1;
    q1.push(sursa);
    while(!q1.empty())
    {
        int nod=q1.front();
        q1.pop();
        vaz[nod]=0;
        for(int i=0;i<v[nod].size();i++)
        {
            int vec=v[nod][i].nod,c=v[nod][i].cost;
            if(cap_flow[nod][vec]>0 && dist[nod]+c<dist[vec])
            {
                dist[vec]=c+dist[nod];
                if(vec!=dest && vaz[vec]==0) {vaz[vec]=1;q1.push(vec);}
            }
        }
    }
}

int dijkstra()
{
    for(int i=1;i<=n;i++) d[i]=distt[i]=inf;
    d[sursa]=distt[sursa]=0;
    q.push({sursa,0});
    while(!q.empty())
    {
        int nod=q.top().nod,c=q.top().cost;
        q.pop();
        if(d[nod]!=c) continue;
        for(int i=0;i<v[nod].size();i++)
        {
            int vec=v[nod][i].nod,cost=v[nod][i].cost;
            if(cap_flow[nod][vec]>0 && d[nod]+dist[nod]+cost-dist[vec]<d[vec])
            {
                d[vec]=d[nod]+dist[nod]+cost-dist[vec];
                distt[vec]=distt[nod]+cost;
                tata[vec]=nod;
                if(vec!=dest) q.push({vec,d[vec]});
            }
        }
    }
    for(int i=1;i<=n;i++) dist[i]=distt[i];
    return d[dest]<inf;
}

int main()
{
    freopen("fmcm.in","r",stdin);
    freopen("fmcm.out","w",stdout);
    int m,x,y,z,c,sol=0;
    scanf("%d%d%d%d",&n,&m,&sursa,&dest);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&x,&y,&z,&c);
        v[x].push_back({y,c});
        v[y].push_back({x,-c});
        cap_flow[x][y]=z;
    }
    bellman_ford();
    while(dijkstra())
    {
        int s=inf;
        for(int nod=dest;nod!=sursa;nod=tata[nod]) s=min(s,cap_flow[tata[nod]][nod]);
        for(int nod=dest;nod!=sursa;nod=tata[nod])
        {
            cap_flow[tata[nod]][nod]-=s;
            cap_flow[nod][tata[nod]]+=s;
        }
        sol+=s*dist[dest];
    }
    printf("%d",sol);
    return 0;
}