[NOIP2017模拟]流

样例数据1
输入
2
1 2 1 2 3
输出
1 4
样例数据2
输入
3
1 2 1 2 3
1 3 0 0 0
2 3 3 4 5
输出
-1 -1
样例数据3
输入
4
1 2 0 2 1
2 3 0 2 1
1 3 0 2 6
1 4 0 0 1
2 4 0 0 0
3 4 2 3 0
输出
2 15
样例数据4
输入
3
1 2 0 2 1
1 3 1 2 1
2 3 1 2 1
输出
2 6

分析：乍眼一看，怎么NOIP考了一道网络流，结果发现网络流不可做，因为网络流解决的是线性规划问题（即给出一堆不等式求最大值），而本题出现了$x^2$。看看数据范围，非常小，就直接暴力dfs。

100分代码：
（细节很多，我打了10遍都WA，确实不想打注释了，累orz）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int INF=0x3f3f3f3f;
int n,m,ans,bz,ansflow=INF,ansmax,l[7][7],h[7][7],low[7],high[7],a[7][7],in[7],out[7];

void dfs(int x,int y,int flow,int pro)
{
    if(flow>ansflow)return;
    if(x==n)
    {
        bz=1;
        if(flow<ansflow)ansflow=flow,ansmax=pro;
        else if(flow==ansflow)ansmax=max(ansmax,pro);
        return;
    }
    if(y>n)
    {
        dfs(x+1,x+2,flow,pro);
        return;
    }
    if(x>1&&y==n)
    {
        int s=in[x]-out[x];
        if(s>h[x][y]||s<l[x][y])return;
        dfs(x,y+1,flow,pro+s*s+(s>0)*(a[x][y]));
        return ;
    }
    for(int i=l[x][y];i<=h[x][y];i++)
    {
        out[x]+=i;in[y]+=i;
        if(x==1)dfs(x,y+1,flow+i,pro+i*i+(i>0)*a[x][y]);
        else dfs(x,y+1,flow,pro+i*i+(i>0)*a[x][y]);
        out[x]-=i;in[y]-=i;
    }
}

int main()
{
    freopen("flow.in","r",stdin);
    freopen("flow.out","w",stdout);
    int x,y;
    n=getint(),m=n*(n-1)/2;
    while(m--)
    {
        x=getint(),y=getint();
        l[x][y]=getint();
        h[x][y]=getint();
        a[x][y]=getint();
    }
    if(n==2)
    {
        if(l[1][2]>h[1][2])
        {
            cout<<-1<<" "<<-1;
            return 0;
        }
        cout<<l[1][2]<<" "<<a[1][2]+l[1][2]*l[1][2];
        return 0;
    }
    if(n==3)
    {
        low[1]=max(l[1][2],l[2][3])+l[1][3];
        high[1]=min(h[1][2],h[2][3])+h[1][3];
        if(l[1][2]>h[2][3]||h[1][2]<l[2][3])
        {
            cout<<-1<<" "<<-1;
            return 0;
        }
        l[1][2]=max(l[1][2],l[2][3]);
        if(l[1][2])ans+=a[1][2]+a[2][3]+2*l[1][2]*l[1][2];
        if(l[1][3])ans+=a[1][3]+l[1][3]*l[1][3];
        cout<<low[1]<<" "<<ans;
        return 0;
    }
    dfs(1,2,0,0);
    if(bz)cout<<ansflow<<" "<<ansmax;
    else cout<<-1<<" "<<-1;
    return 0;
}

本题结。