树上差分水题:

一、Network POJ - 3417:求多少条树边被非树边覆盖了零次或者一次:注意,若一刀就能使树不连通,那么也需要第二刀去切附加边。1号节点的 cf 不算#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#include<cmath&...
摘要由CSDN通过智能技术生成

一、Network POJ - 3417:
求多少条树边被非树边覆盖了零次或者一次:
注意,若一刀就能使树不连通,那么也需要第二刀去切附加边。
1号节点的 cf 不算

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
#include<queue>
#include<set>
#define ll long long
#define llu unsigned ll
using namespace std;
const int mod=1000000007;
const int maxn=200100;
const int inf=0x3f3f3f3f;
int head[maxn],edge[maxn],ver[maxn],nt[maxn];
int f[maxn][20],d[maxn];
int cf[maxn];
int tot=1;
int n,m;
int t;
ll cnt=0;

void init(void)
{
    memset(head,0,sizeof(head));
    memset(cf,0,sizeof(cf));
    memset(f,0,sizeof(f));
    memset(d,0,sizeof(d));
    tot=1;cnt=0;
}

void add(int x,int y)
{
    ver[++tot]=y,nt[tot]=head[x],head[x]=tot;
}

void bfs(void)
{
    queue<int>q;
    q.push(1);
    d[1]=1;
    while(q.size())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i;i=nt[i])
        {
            int y=ver[i];
            if(d[y]) continue;
            d[y]=d[x]+1;
            f[y][0]=x;
            for(int j=1;j<=t;j++)
                f[y][j]=f[f[y][j-1]][j-1];
            q.push(y);
        }
    }
    return ;
}

int lca(int x,int y)
{
    if(d[x]>d[y]) swap(x,y);
    for(int i=t;i>=0;i--)
        if(d[f[y][i]]>=d[x]) y=f[y][i];

    if(x==y) return x;

    for(int i=t;i>=0;i--)
        if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}

void dfs(int x,int fa)
{

    for(int i=head[x];i;i=nt[i])
    {
        int y=ver[i];
        if(y==fa) continue;
        dfs(y,x);
        if(cf[y]==0) cnt=cnt+m;
        else if(cf[y]==1) cnt++;
        //cout<<"y:  "<<y<<"   cf[y]:  "<<cf[y]<<endl;
        cf[x]+=cf[y];
    }
    return ;
}


int main(void)
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        t=log2(n)+1;
        int x,y;
        for(int i=1;i<n;i++)
        {
            scanf("%d%d",&x,&y);
            add(x,y);
            add(y,x);
        }
        bfs();

        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&x,&y);
            cf[x]++,cf[y]++,cf[lca(x,y)]-=2;
        }

        dfs(1,0);


        printf("%lld\n",cnt);
    }
    return 0;
}

二、P3128 [USACO15DEC]最大流Max Flow
随便搞搞就可以了:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
#include<queue>
#include<set>
#define ll long long
#define llu unsigned ll
using namespace std;
const int mod=1000000007;
const int maxn=200100;
const int inf=0x3f3f3f3f;
int head[maxn],edge[maxn],ver[maxn],nt[maxn];
int f[maxn][20],d[maxn];
int cf[maxn];
int tot=1;
int n,m;
int t;
ll cnt=0;

void init(void)
{
    memset(head,0,sizeof(head));
    memset(cf,0,sizeof(cf));
    memset(f,0,sizeof(f));
    memset(d,0,sizeof(d));
    tot=1;cnt=0;
}

void add(int x,int y)
{
    ver[++tot]=y,nt[tot]=head[x],head[x]=tot;
}

void bfs(void)
{
    queue<int>q;
    q.push(1);
    d[1]=1;
    while(q.size())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i;i=nt[i])
        {
            int y=ver[i];
            if(d[y]) continue;
            d[y]=d[x]+1;
            f[y][0]=x;
            for(int j=1;j<=t;j++)
                f[y][j]=f[f[y][j-1]][j-1];
            q.push(y);
        }
    }
    return ;
}

int lca(int x,int y)
{
    if(d[x]>d[y]) swap(x,y);
    for(int i=t;i>=0;i--)
        if(d[f[y][i]]>=d[x]) y=f[y][i];

    if(x==y) return x;

    for(int i=t;i>=0;i--)
        if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}

void dfs(int x,int fa)
{

    for(int i=head[x];i;i=nt[i])
    {
        int y=ver[i];
        if(y==fa) continue;
        dfs(y,x);
        cf[x]+=cf[y];
    }
    cnt=max(cnt,(ll)cf[x]);
    return ;
}


int main(void)
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        t=log2(n)+1;
        int x,y;
        for(int i=1;i<n;i++)
        {
            scanf("%d%d",&x,&y);
            add(x,y);
            add(y,x);
        }
        bfs();

        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&x,&y);
            cf[x]++,cf[y]++;
            int pp=lca(x,y);
            cf[pp]--,cf[f[pp][0]]--;
        }

        dfs(1,0);


        printf("%lld\n",cnt);
    }
    return 0;
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值