ZOJ 3613 Wormhole Transport（DP 斯坦纳树）

题意：有n个星球，其中有些星球上有多个工作，有些星球上有些资源（但是一个星球上的资源只能提供给一个工厂），知道在一些星球边建立边的代价，问使得得到资源的工厂的数目最多的多少，在些前提下代价最小是多少。

还是斯坦纳树，不懂看：http://endlesscount.blog.163.com/blog/static/821197872012525113427573/

#include <iostream>#include <cstdio>#include <cstring>#include <queue>using namespace std;const int N=205;struct Edge{int v,w;Edge *nxt;}memo[N*N],*cur,*adj[N];int n,m,mask;int d[N][1<<10],dp[1<<10];bool vis[N][1<<10];int fac[N],rec[N],st[N];queue<int> que;void addEdge(int u,int v,int w){cur->v=v;cur->w=w;cur->nxt=adj[u];adj[u]=cur++;}void up(int &a,int b){ a=min(a,b); }bool judge(int s){int cnt=0;for(int i=1;i<=n;i++){if(s&st[i]) cnt+=fac[i],cnt-=rec[i];}return cnt>=0;}int getnum(int s){int cnt=0;for(int i=1;i<=n;i++){if(s&st[i]) cnt+=rec[i];}return cnt;}void spfa(){while(que.size()){int u=que.front()/1000,s=que.front()%1000;que.pop();vis[u][s]=0;for(Edge* it=adj[u];it;it=it->nxt){int v=it->v,w=it->w;if(d[v][s|st[v]]>d[u][s]+w){d[v][s|st[v]]=d[u][s]+w;if((s|st[v])==s&&!vis[v][s])que.push(v*1000+s),vis[v][s]=1;}}}}void init(){cur=memo;memset(adj,0,sizeof(adj));memset(vis,0,sizeof(vis));memset(dp,63,sizeof(dp));memset(d,63,sizeof(d));memset(st,0,sizeof(st));}int main(){while(scanf("%d",&n)!=EOF){init();int ans=0,fac_cnt=0,rec_cnt=0;for(int i=1;i<=n;i++){scanf("%d%d",&fac[i],&rec[i]);if(fac[i]&&rec[i]) fac[i]--,rec[i]--,ans++;if(fac[i]) {st[i]=(1<<(fac_cnt++));d[i][st[i]]=0;}if(rec[i]){st[i]=(1<<(4+rec_cnt++));d[i][st[i]]=0;}}scanf("%d",&m);while(m--){int u,v,w;scanf("%d%d%d",&u,&v,&w);addEdge(u,v,w);addEdge(v,u,w);}mask=(1<<(4+rec_cnt))-1;for(int s=1;s<=mask;s++){for(int i=1;i<=n;i++){if(st[i]&&!(st[i]&s)) continue;for(int p=(s-1)&s;p;p=(p-1)&s){int s1=p|st[i],s2=(s-p)|st[i];up(d[i][s],d[i][s1]+d[i][s2]);}if(d[i][s]<1e9&&!vis[i][s])que.push(i*1000+s),vis[i][s]=1;}spfa();}for(int s=1;s<=mask;s++)for(int i=1;i<=n;i++)dp[s]=min(dp[s],d[i][s]);int cnt=0,cost=0;for(int s=1;s<=mask;s++){if(judge(s)&&dp[s]<1e9){for(int p=(s-1)&s;p;p=(p-1)&s){if(judge(p)&&judge(s-p)&&dp[p]<1e9&&dp[s-p]<1e9)up(dp[s],dp[p]+dp[s-p]);}int tmp=getnum(s);if(cnt<tmp||(cnt==tmp&&cost>dp[s]))cnt=tmp,cost=dp[s];}}printf("%d %d\n",cnt+ans,cost);}return 0;}