BZOJ 1001 [BeiJing2006] 狼抓兔子(平面图最大流)
2015-02-03 21:12
337 查看
边数巨多, 好一点的最大流模板可能会过,但这个题的正解是平面最大流,具体讲解百度就好了,周冬的课件讲的非常详细
然后 + 一个队列优化的DIJ就好了, 不优化的话会TLE, 还要注意的就是n = 1 和 m = 1的情况
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <vector>
#include <queue>
using namespace std;
#define MAXM 2111111
#define MAXN 2111111
#define INF 0x1f1f1f1f
struct Edge
{
int v, val;
int next;
} E[MAXM << 2];
bool vis[MAXN];
int dis[MAXN];
int head[MAXN], ecnt ;
deque<int> Q;
void InitGra( )
{
ecnt = 0;
memset(head, -1, sizeof(head));
}
int sta, end;
void AddEdge( int u, int v, int val)
{
E[ecnt].v = v;
E[ecnt].val = val;
E[ecnt].next = head[u];
head[u] = ecnt++;
}
int Dijkstra( int n )
{
memset(vis, 0, sizeof(vis));
memset(dis, INF, sizeof(dis));
while(!Q.empty()) Q.pop_front();
vis[sta] = 1;
dis[sta] = 0;
Q.push_back(sta);
while(!Q.empty())
{
int u=Q.front();Q.pop_front();
for(int k=head[u];k!=-1;k=E[k].next)
{
int v=E[k].v;
if(dis[v]>dis[u]+E[k].val)
{
dis[v]=dis[u]+E[k].val;
if(!vis[v])
{
vis[v]=true;
if(!Q.empty())
dis[v]>dis[Q.front()]?Q.push_back(v):Q.push_front(v);
else Q.push_back(v);
}
}
}
vis[u]=false;
}
printf("%d\n",dis[end]);
return 0;
}
int main( )
{
int u, v, val;
int n, m;
while(scanf("%d %d",&n, &m) != EOF)
{
InitGra();
int num = n*m+2;
int u, v, val;
end = num ;
sta = num - 1;
int ans = INF;
if( n == 1 || m == 1)
{
for(int i=1; i<=n; i++)
for(int j=1; j<m; j++)
{
scanf("%d",&val);
ans=min(ans,val);
}
for(int i=1; i<n; i++)
for(int j=1; j<=m; j++)
{
scanf("%d",&val);
ans=min(ans,val);
}
printf("%d\n",ans);
continue;
}
for(int i = 1; i <= n ; i++)
{
for( int j = 1; j < m ; j++)
{
scanf("%d",&val);
if( i == 1) u = sta, v = j;
else if(i == n) u = (i*2 - 3)*(m-1) + j, v = end;
else u = (i*2 - 3)*(m-1) + j, v = (i-1)*2*(m-1) + j;
AddEdge(u ,v, val);
}
}
for( int i = 1; i < n; i++)
{
for( int j = 1; j <= m; j++)
{
scanf("%d",&val);
if(j == 1) u =(i*2 - 1)*(m-1) + 1, v = end;
else if(j == m) u = sta, v = (i*2 - 1)*(m-1);
else u = (i-1)*2*(m-1) + j - 1 , v = (i*2 - 1)*(m-1)+ j;
AddEdge(u ,v, val);
}
}
for( int i = 1; i < n; i++)
{
for( int j = 1; j < m; j++)
{
scanf("%d",&val);
u = (i-1)*2*(m-1) + j, v = (i*2 - 1)*(m-1) +j;
AddEdge(u ,v, val);
}
}
Dijkstra(end);
}
return 0;
}
然后 + 一个队列优化的DIJ就好了, 不优化的话会TLE, 还要注意的就是n = 1 和 m = 1的情况
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <vector>
#include <queue>
using namespace std;
#define MAXM 2111111
#define MAXN 2111111
#define INF 0x1f1f1f1f
struct Edge
{
int v, val;
int next;
} E[MAXM << 2];
bool vis[MAXN];
int dis[MAXN];
int head[MAXN], ecnt ;
deque<int> Q;
void InitGra( )
{
ecnt = 0;
memset(head, -1, sizeof(head));
}
int sta, end;
void AddEdge( int u, int v, int val)
{
E[ecnt].v = v;
E[ecnt].val = val;
E[ecnt].next = head[u];
head[u] = ecnt++;
}
int Dijkstra( int n )
{
memset(vis, 0, sizeof(vis));
memset(dis, INF, sizeof(dis));
while(!Q.empty()) Q.pop_front();
vis[sta] = 1;
dis[sta] = 0;
Q.push_back(sta);
while(!Q.empty())
{
int u=Q.front();Q.pop_front();
for(int k=head[u];k!=-1;k=E[k].next)
{
int v=E[k].v;
if(dis[v]>dis[u]+E[k].val)
{
dis[v]=dis[u]+E[k].val;
if(!vis[v])
{
vis[v]=true;
if(!Q.empty())
dis[v]>dis[Q.front()]?Q.push_back(v):Q.push_front(v);
else Q.push_back(v);
}
}
}
vis[u]=false;
}
printf("%d\n",dis[end]);
return 0;
}
int main( )
{
int u, v, val;
int n, m;
while(scanf("%d %d",&n, &m) != EOF)
{
InitGra();
int num = n*m+2;
int u, v, val;
end = num ;
sta = num - 1;
int ans = INF;
if( n == 1 || m == 1)
{
for(int i=1; i<=n; i++)
for(int j=1; j<m; j++)
{
scanf("%d",&val);
ans=min(ans,val);
}
for(int i=1; i<n; i++)
for(int j=1; j<=m; j++)
{
scanf("%d",&val);
ans=min(ans,val);
}
printf("%d\n",ans);
continue;
}
for(int i = 1; i <= n ; i++)
{
for( int j = 1; j < m ; j++)
{
scanf("%d",&val);
if( i == 1) u = sta, v = j;
else if(i == n) u = (i*2 - 3)*(m-1) + j, v = end;
else u = (i*2 - 3)*(m-1) + j, v = (i-1)*2*(m-1) + j;
AddEdge(u ,v, val);
}
}
for( int i = 1; i < n; i++)
{
for( int j = 1; j <= m; j++)
{
scanf("%d",&val);
if(j == 1) u =(i*2 - 1)*(m-1) + 1, v = end;
else if(j == m) u = sta, v = (i*2 - 1)*(m-1);
else u = (i-1)*2*(m-1) + j - 1 , v = (i*2 - 1)*(m-1)+ j;
AddEdge(u ,v, val);
}
}
for( int i = 1; i < n; i++)
{
for( int j = 1; j < m; j++)
{
scanf("%d",&val);
u = (i-1)*2*(m-1) + j, v = (i*2 - 1)*(m-1) +j;
AddEdge(u ,v, val);
}
}
Dijkstra(end);
}
return 0;
}
相关文章推荐
- bzoj 1001: [BeiJing2006]狼抓兔子(s-t平面图求最大流)
- bzoj1001 [BeiJing2006]狼抓兔子(最小割/平面图最大流转对偶图最短路)
- BZOJ 1001 [BeiJing2006] 狼抓兔子(平面图最大流)
- _bzoj1001 [BeiJing2006]狼抓兔子【平面图】
- 【平面图转对偶图】【最短路】【Beijing 2006】【bzoj 1001】狼抓兔子
- 【BZOJ】 1001 [BeiJing2006]狼抓兔子 网络最大流
- BZOJ 1001: [BeiJing2006]狼抓兔子(平面图最小割转对偶图最短路)
- bzoj1001: [BeiJing2006]狼抓兔子(最大流)
- bzoj 1001: [BeiJing2006]狼抓兔子 平面图最小割
- BZOJ 1001: [BeiJing2006]狼抓兔子【最大流/SPFA+最小割,多解】
- BZOJ1001 [BeiJing2006]狼抓兔子(平面图最小割转最短路)
- [bzoj1001][BJOI2006]狼抓兔子——最大流转最短路,平面图
- 【平面图】【最小割】【最短路】【Heap-Dijkstra】bzoj1001 [BeiJing2006]狼抓兔子
- 【BZOJ1001】【Beijing2006】狼抓兔子(平面图转对偶图:最小割+最短路)
- BZOJ 1001([BeiJing2006]狼抓兔子-最大流转对偶图最短路)
- 【平面图最小割】BZOJ1001- [BeiJing2006]狼抓兔子
- BZOJ 1001: [BeiJing2006]狼抓兔子(s-t平面图+最短路求最小割)
- BZOJ 1001([BeiJing2006]狼抓兔子-最大流转对偶图最短路)
- 【bzoj1001】【BeiJing2006】【狼抓兔子】【网络流】【平面图转对偶图】
- [平面图转最短路]BZOJ1001:[BeiJing2006]狼抓兔子