Codeforces #81E Pairs

2 篇文章 0 订阅

dp1[x][0]表示x的子树所能组合的最大对数,dp1[x][1]表示x的子树所能组合的最大对数的情况下的最多 男女 数

dp2[x][0]表示x和它的子树所能组合的最大对数,dp2[x][1]表示x和它的子树所能组合的最大对数的情况下的最多 男女 数

得:

dp1[x][0]+=dp2[y][0];    //y为x的子结点

dp1[x][1]+=dp2[y][1];

dp2[x][0]=dp1[x][0]+增量;

dp2[x][1]=dp1[x][1]+增量;

每个圈里的数字匹配方式有两种,故需要做两次DP。

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<map>
#include<vector>
using namespace std;

const int N=100010;

int n,dp1[N][2],dp2[N][2],ans[N][2],len;
int tmp[N][2],w;
int head[N],adj[N],v[N],s[N],g[N];
int root,ansa,ansb;
bool vis[N];

void DFS(int x)
{
    vis[x]=true;
    g[x]=0;
    dp1[x][0]=dp1[x][1]=0;
    dp2[x][0]=dp2[x][1]=0;
    int a,b;
    for(int y=head[x];y!=-1;y=adj[y]) {
        if(y!=root) {
            DFS(y);
            dp1[x][0]=dp1[x][0]+dp2[y][0];
            dp1[x][1]=dp1[x][1]+dp2[y][1];
            a=dp1[y][0]+1-dp2[y][0];
            b=dp1[y][1]+(s[x]^s[y])-dp2[y][1];
            if(a>dp2[x][0]||(a==dp2[x][0]&&b>dp2[x][1])) {
                dp2[x][0] = a;
                dp2[x][1] = b;
                g[x]=y;
            }
        }
    }
    dp2[x][0]+=dp1[x][0];
    dp2[x][1]+=dp1[x][1];
}

void calculate(int x,int flag)
{
    for(int y=head[x];y!=-1;y=adj[y]) {
        if(y!=root) {
            if(flag&&g[x]==y) {
                tmp[w][0]=x;
                tmp[w][1]=y;
                ++w;
                calculate(y,0);
            }
            else calculate(y,1);
        }
    }
}

void work(int x)
{
    while(!vis[x]) {
        vis[x]=true;
        x=v[x];
    }
    int cnt = 2;
    int a=0,b=0;
    while(cnt--) {
        root = x;
        DFS(x);
        if(dp2[x][0]>a||(dp2[x][0]==a&&dp2[x][1]>b)) {
            a=dp2[x][0];
            b=dp2[x][1];
            w=0;
            calculate(x,1);
        }
        x=v[x];
    }
    ansa+=a;
    ansb+=b;
    for(int i=0;i<w;++i) {
        ans[len][0]=tmp[i][0];
        ans[len][1]=tmp[i][1];
        ++len;
    }
}

int main()
{
    memset(head,-1,sizeof head);
    cin>>n;
    for(int i=1;i<=n;++i) {
        scanf("%d%d",v+i,s+i);
        --s[i];
        adj[i]=head[v[i]];
        head[v[i]]=i;
    }
    for(int i=1;i<=n;++i)
        if(!vis[i]) work(i);
    printf("%d %d\n",ansa,ansb);
    for(int i=0;i<len;++i)
        printf("%d %d\n",ans[i][0],ans[i][1]);
}


 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值