2015-2016 Northwestern European Regional Contest (NWERC 2015) E

题外话：

vj上的题交到的网站不知道怎么回事数据都炸了。。。还是去交codeforces吧。。
ps：确切的说是交到icpc？网站的都炸了，，
题目链接：
http://codeforces.com/gym/101485

我是用网络流做的，后来听说官方题解是二分图匹配，不过都一样啦。
s连2500个点（容量为1），然后找2500个点分别去连他们能得到的三种数值，然后对于每种数值，连t（容量为1）。跑一遍最大流，如果不等于n（点的个数），那么无解。如果等了那么就去找对应的那条边就好了。

#include <string.h>#include <iostream>#include <cstdio>#include <map>using namespace std;typedef long long int lli;#define inf 0x3f3f3f3fconst lli maxn = 22000;struct edge{    lli to,next,v;}ed[maxn*10];lli d[maxn],cur[maxn],pre[maxn],gap[maxn],q[maxn+2000],cnte,head[maxn];void ae(lli x, lli y, lli v) {    ed[cnte].to = y;    ed[cnte].v = v;    ed[cnte].next = head[x];    head[x] = cnte++;    ed[cnte].to = x;    ed[cnte].v = 0;    ed[cnte].next = head[y];    head[y] = cnte++;}void rbfs (lli s,lli t) {    lli fi,se;    memset(gap,0,sizeof(gap));    memset(d,-1,sizeof(d));    d[t] = 0;    gap[0] = 1;    fi = se = 0;    q[se++] = t;    while (fi != se) {        lli u = q[fi++];        for (lli i=head[u];~i;i=ed[i].next) {            lli v = ed[i].to;            if (~d[v]) continue;            d[v] = d[u] + 1;            q[se++] = v;            gap[d[v]]++;        }    }}lli isap(lli s,lli t){    memcpy(cur,head,sizeof(head));    rbfs (s,t);    lli flow = 0,u = pre[s] = s,i;    while(d[s] < t+1) {        if(u==t) {            lli f = inf,neck;            for(i=s;i!=t;i=ed[cur[i]].to){                if(f > ed[cur[i]].v){                    f = ed[cur[i]].v;                    neck = i;                }            }            for(i=s;i!=t;i=ed[cur[i]].to){                ed[cur[i]].v -= f;                ed[cur[i]^1].v += f;            }            flow += f;            u = neck;        }        for(i=cur[u];~i;i=ed[i].next) if(d[ed[i].to]+1 == d[u] && ed[i].v) break;        if(~i) {            cur[u] = i;            pre[ed[i].to] = u;            u = ed[i].to;        }        else{            if(gap[d[u]]==0||0 == (--gap[d[u]]) ) break;            lli mind = t+1;            for(i=head[u];~i;i=ed[i].next){                if(ed[i].v && mind>d[ed[i].to]){                    cur[u] = i;                    mind = d[ed[i].to];                }            }            d[u] = mind + 1;            gap[d[u]]++;            u = pre[u];        }    }    return flow;}struct pp{    int a,b;}res[2600];lli fx[11000];void ini(){    memset(head,-1,sizeof(head));cnte = 0;}int main(){    int n,s = 0,t = 10200;    while(~scanf("%d",&n)){        ini();        map<lli,int> ma;int order = 2500;        for(int i = 1;i <= n;i++){           scanf("%d%d",&res[i].a,&res[i].b);           lli t1 = res[i].a+res[i].b, t2 = res[i].a-res[i].b, t3 = (lli)res[i].a*(lli)res[i].b;            if(!ma[t1]){                ma[t1] = ++order;fx[order] = t1;            }            if(!ma[t2]){                ma[t2] = ++order;fx[order] = t2;            }            if(!ma[t3]){                ma[t3] = ++order;fx[order] = t3;            }        }        for(int i = 1;i <= n;i++){            lli t1 = res[i].a+res[i].b, t2 = res[i].a-res[i].b, t3 = (lli)res[i].a*(lli)res[i].b;            ae(s,i,1);            ae(i,ma[t1],1);            if(t2 != t1)                ae(i,ma[t2],1);            if(t3 != t2 && t3 != t1)                ae(i,ma[t3],1);        }        for(int i = 2501;i <= order;i++){            ae(i,t,1);        }        if(isap(s,t) == n){            for(int i = 1;i <= n;i++){                for(int j = head[i];~j;j=ed[j].next){                    int to = ed[j].to;                    if(ed[j].v == 0){                        if(res[i].a+res[i].b == fx[to]){                            printf("%d + %d = %I64d\n",res[i].a,res[i].b,fx[to]);                        }                        else if(res[i].a-res[i].b == fx[to]){                            printf("%d - %d = %I64d\n",res[i].a,res[i].b,fx[to]);                        }                        else if((lli)res[i].a*(lli)res[i].b == fx[to]){                            printf("%d * %d = %I64d\n",res[i].a,res[i].b,fx[to]);                        }                        break;                    }                }            }        }        else{            puts("impossible");        }    }}