AGC 004 F Namori

神仙题

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<stack>
#include<assert.h>
#define gc getchar()
#define N 100010
#define clr(a,n) memset(a,0,sizeof(int)*((n)+1))
#define build_edge(x,y) add_edge(x,y),add_edge(y,x)
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
inline int inn()
{
    int x,ch;while((ch=gc)<'0'||ch>'9');
    x=ch^'0';while((ch=gc)>='0'&&ch<='9')
        x=(x<<1)+(x<<3)+(ch^'0');return x;
}
inline int gabs(int x) { return x<0?-x:x; }
struct edges{
    int to,pre;
}e[N<<1];int h[N],etop,find_cir,s[N],onc[N],vis[N];
int val[N],u[N],v[N],a[N],d[N],k[N];stack<int> stc;
inline int add_edge(int u,int v) { return e[++etop].to=v,e[etop].pre=h[u],h[u]=etop; }
inline int getdpt(int x,int dpt=0)
{
    vis[x]=1;
    for(int i=h[x],y;i;i=e[i].pre)
        if(!vis[y=e[i].to]) getdpt(y,dpt^1);
    return val[x]=((d[x]=dpt)?1:-1);
}
inline int getans(int x,int ans=0)
{
    vis[x]=1,s[x]=val[x];
    for(int i=h[x],y;i;i=e[i].pre)
        if(!vis[y=e[i].to]) ans+=getans(y),s[x]+=s[y],k[x]+=k[y];
    return ans+=gabs(s[x]);
}
int dfs(int x,int fa=0)
{
    if(find_cir) return 0;
    stc.push(x),vis[x]=1;
    for(int i=h[x],y;i;i=e[i].pre)
        if((y=e[i].to)^fa)
        {
            if(find_cir) return 0;
            if(!vis[y]) dfs(y,x);
            else{
                find_cir=1,onc[x]=onc[y]=1;
                while(stc.top()!=y) onc[stc.top()]=1,stc.pop();
            }
        }
    return (stc.top()==x?stc.pop(),0:0),0;
}
int main()
{
    int n=inn(),m=inn(),ans=0,c=0,x=0;
    for(int i=1;i<=m;i++)
        u[i]=inn(),v[i]=inn(),build_edge(u[i],v[i]);

    getdpt(1),clr(vis,n),ans=getans(1),clr(vis,n);
    if(m<n) return !printf("%d\n",s[1]?-1:ans);

    dfs(1),clr(vis,n),clr(h,n),etop=0;
    for(int i=1;i<=n;i++) if(onc[i]) c++;
    for(x=1;!onc[u[x]]||!onc[v[x]];x++);
    for(int i=1;i<=m;i++) if(x^i) build_edge(u[i],v[i]);

    int rt=u[x];getdpt(1),clr(vis,n),getans(rt),clr(vis,n);
    if(c&1)
    {
        if(gabs(s[rt])&1) return !printf("-1\n");
        val[u[x]]-=s[rt]/2,val[v[x]]-=s[rt]/2,ans=gabs(s[rt])/2;
        return !printf("%d\n",ans+getans(rt));
    }
    if(s[rt]) return !printf("-1\n");
    k[u[x]]=1,k[v[x]]=-1,getans(rt),ans=0,c=0;
    for(int i=1;i<=n;i++)
        if(!k[i]) ans+=gabs(s[i]);
        else a[++c]=(k[i]>0?s[i]:-s[i]);
    a[++c]=0,sort(a+1,a+c+1),x=a[(c+1)/2];
    for(int i=1;i<=c;i++) ans+=gabs(a[i]-x);
    return !printf("%d\n",ans);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值