題意:一個n個點m條邊的無向無環圖,在盡量少的節點上放燈,使得所有邊都被照亮,每盞燈將照亮以它為一個端點的所有邊,在總燈數最小的前提下,被兩盞燈同時照亮的邊數應盡量大。
思路:無向無環圖就是“森林”,常用樹形dp,本題要優化的目標有兩個,放置的燈數a應盡量少,被兩盞燈同時照亮的邊數b應盡量大,為了統一,我們把b替換成”恰好被一盞燈照亮的邊數c盡量小“。然後設x=Ma+c為最終的優化目標,M是一個很大的正整數。當x取最小值的時候,x/M就是a的最小值,x%M就是c最小值。
定義dp(i,j),其中i表示節點i,j表示節點i的父節點是否放置了街燈,0代表沒放,1代表放了,則dp(i,j)代表在上述下x的最小值。
實際上,對於每個節點而言,只有兩種決策:在i處放或者不放街燈。
決策一:節點i處不放街燈,那麼i是根或者父親節點放了街燈。所以dp(i,j)=sum{ dp(v,0) | v取遍i的所有兒子節點 },如果i不是根節點,那麼結果+1,因為i和父親連接的這條邊只被一盞燈照亮。
決策二:節點i處放街燈,dp(i,j)=sum{ dp(v,1)| v取遍i的所有兒子節點 } + M,如果i不是根節點而且j=0,那麼結果+1。
總結:以後遇到需要同時優化兩個量v1,v2的問題,要求首先滿足v1最小,在這個前提下v2最小的問題,可以考慮優化x=M*v1+v2,其中M是比"比v2的最大理論值和v2的最小理論值之差"還要大的數。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<vector>
#include<queue>
#include<stack>
#define INF 0x3f3f3f3f
#define NODENUM 1005
#define EDGENUM 1005
#define MAXN 1005
using namespace std;
int root;
const int m=2000;
struct EdgeNode{int to,next;} E[2*EDGENUM];
int edgenum,head[NODENUM],N,T,M;
bool vis[NODENUM];
int ans,dp[NODENUM][2];
void init()
{
edgenum=0;
memset(head,-1,sizeof(head));
memset(vis,0,sizeof(vis));
ans=0;
}
void add(int x,int y)
{
edgenum++;
E[edgenum].next=head[x];
head[x]=edgenum;
E[edgenum].to=y;
}
void dfs(int s)
{
vis[s]=1;
int sum0=0,sum1=0;
for(int p=head[s];p!=-1;p=E[p].next)
{
int v=E[p].to;
if(!vis[v])
{
dfs(v);
sum0+=dp[v][0];
sum1+=dp[v][1];
}
}
if(s==root) dp[s][0]=min(sum1+m,sum0),ans+=dp[s][0];
else dp[s][1]=min(sum0+1,sum1+m),dp[s][0]=sum1+m+1;
}
void build()
{
scanf("%d%d",&N,&M);
for(int i=0;i<M;++i)
{
int a,b;
scanf("%d%d",&a,&b);
add(a,b);
add(b,a);
}
}
int main()
{
scanf("%d",&T);
while(T--)
{
init();
build();
for(int i=0;i<N;++i) if(!vis[i]) dfs(root=i);
printf("%d %d %d\n",ans/m,M-ans%m,ans%m);
}
return 0;
}
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<vector>
#include<queue>
#include<stack>
#define INF 0x3f3f3f3f
#define NODENUM 1005
#define EDGENUM 1005
#define MAXN 1005
using namespace std;
int root;
const int m=2000;
struct EdgeNode{int to,next;} E[2*EDGENUM];
int edgenum,head[NODENUM],N,T,M;
bool vis[NODENUM];
int ans,dp[NODENUM][2];
void init()
{
edgenum=0;
memset(head,-1,sizeof(head));
memset(vis,0,sizeof(vis));
ans=0;
}
void add(int x,int y)
{
edgenum++;
E[edgenum].next=head[x];
head[x]=edgenum;
E[edgenum].to=y;
}
void dfs(int s)
{
vis[s]=1;
int sum0=0,sum1=0;
for(int p=head[s];p!=-1;p=E[p].next)
{
int v=E[p].to;
if(!vis[v])
{
dfs(v);
sum0+=dp[v][0];
sum1+=dp[v][1];
}
}
if(s==root) dp[s][0]=min(sum1+m,sum0),ans+=dp[s][0];
else dp[s][1]=min(sum0+1,sum1+m),dp[s][0]=sum1+m+1;
}
void build()
{
scanf("%d%d",&N,&M);
for(int i=0;i<M;++i)
{
int a,b;
scanf("%d%d",&a,&b);
add(a,b);
add(b,a);
}
}
int main()
{
scanf("%d",&T);
while(T--)
{
init();
build();
for(int i=0;i<N;++i) if(!vis[i]) dfs(root=i);
printf("%d %d %d\n",ans/m,M-ans%m,ans%m);
}
return 0;
}