#include<cstdio>
#include<deque>
#include<cstring>
#include<queue>
#include<algorithm>
#include<vector>
const int nmax=350;
const int mmax=12500;
const int inf=2000000000;
using namespace std;
int n,m,s,t;
int flux[nmax+5][nmax+5];
int ca[nmax+5][nmax+5],co[nmax+5][nmax+5],ta[nmax+5],cv[nmax+5];
deque<int> que;
bool be[nmax+5],inq[nmax+5];
vector<int > ve[nmax+5];
priority_queue<pair<int,int> ,vector<pair<int,int> >,greater<pair<int, int> > > pq;
void belman()
{
    int x,i;
     memset(cv,inf,sizeof(cv));
    cv[s]=0;
    que.push_back(s);
    inq[s]=1;
    while(!que.empty())
    {
        x=que.back();
        que.pop_back();
        for(i=ve[x].size()-1; i>=0; i--)
            if(ca[x][ve[x][i]]>0 && cv[x]+co[x][ve[x][i]]<cv[ve[x][i]])
            {
                cv[ve[x][i]]=cv[x]+co[x][ve[x][i]];
                if(inq[ve[x][i]]==0)
                {
                    inq[ve[x][i]]=1;
                    que.push_back(ve[x][i]);
                }
            }
        inq[x]=0;
    }

}
int oldd[nmax+5],newd[nmax+5],reald[nmax+5];
bool djkstra()
{
    int i,x,y,nd,cs;
   for(i=1;i<=n;i++)
   newd[i]=inf;
    reald[s]=newd[s]=0;
    pq.push(make_pair(newd[s],s));
    while(!pq.empty())
    {
        cs=pq.top().first;
        y=pq.top().second;
        pq.pop();
        if(cs==newd[y])
        {
        for(i=ve[y].size()-1;i>=0;i--)
        {
            nd=ve[y][i];
        if(ca[y][nd]-flux[y][nd]>0)
        {
            x=newd[y]+co[y][nd]-(oldd[nd]-oldd[y]);
            if(x<newd[nd])
            {
                newd[nd]=x;
                ta[nd]=y;
                reald[nd]=reald[y]+co[y][ve[y][i]];
                pq.push(make_pair(newd[nd],nd));
            }
        }
        }
        }
    }
memcpy(oldd,reald,sizeof(oldd));
    if(newd[t]!=inf)
    return 1;
    else
    return 0;
}
void da()
{
    int i,cmin,z;
    belman();
    for(i=1;i<=n;i++)
    oldd[i]=cv[i];
    while(djkstra())
    {
        cmin=-1;
        z=t;
        while(z!=s)
        {
            if(ca[ta[z]][z]-flux[ta[z]][z]<cmin || cmin==-1)
                cmin=ca[ta[z]][z]-flux[ta[z]][z];
            z=ta[z];
        }
        z=t;
        while(z!=s)
        {
            flux[ta[z]][z]+=cmin;
            flux[z][ta[z]]-=cmin;
            z=ta[z];
        }
    }
}
int main()
{
    freopen("fmcm.in","r",stdin);
    freopen("fmcm.out","w",stdout);
    int i,j,x,y,cmin;
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(i=1; i<=m; i++)
    {
        scanf("%d%d",&x,&y);
        scanf("%d%d",&ca[x][y],&co[x][y]);
        co[y][x]=-co[x][y];
        ve[x].push_back(y);
        ve[y].push_back(x);
    }
    da();
    cmin=0;
    for(i=1; i<=n; i++)
        for(j=1; j<=n; j++)
            if(flux[i][j]>0)
            {
                cmin=cmin+flux[i][j]*co[i][j];
            }
    printf("%d\n",cmin);
    return 0;
}