一、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;
}