poj 3621(参数搜索+二分+spfa负环求解)

原创 2012年07月02日 15:44:35
【题目大意】
          给出一个有向图,问求一个回路,使得回路上的点权之和/边权之和 最大。
【解题思路】
          转:此题是对01分数规划的应用,那么首先明白01分数规划的思想.
01整数规划问题就是求解方程(a1*x1+a2*x2+..+an*xn)/(b1*x1+b2*x2+..+bn*xn)的最小值/最大值问题。其中xi = 0或1(i=1,2...n)对于此类问题我们可以通过二分来求解很接近答案的近似值。我们可以先令:
(a1*x1+a2*x2+..+an*xn)/(b1*x1+b2*x2+..+bn*xn)=L,则我们可以将此式转换为:x1*(a1-b1*L)+x2*(a2-b2*L)+...xn*(an-bn*L)=0,我们先定义一个估计值val,如果这个值使得上面的式子小于0我们就可以知道val>L,如果上式等于0,则val = L;如果大于0,则val<L,显然我们可以采用二分的思想求解次问题。

         对于此题,设happy[u]为点u的欢乐值,w[u][v]为u-->v的边权值。要求的是happy[1]+happy[2]+...+happy[n] / w[1][2]+...+w[n-1][n] = ans,设ans就是所求的最大值。则移项,ans*w[u][v] - happy[v]  = 0 .

          建图:我们重新构造一幅图,使得边权为happy[v] - ans*w[u][v]。用SPFA算法,二分枚举ans,判断是否存在负权回路,若存在,说明ans偏小了,则增大ans,若不存在,则减小ans。

#include <iostream>
#include <string>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const double eps=1e-3;
const int inf=100000000000;
struct nod
{
    int u;
    int len;
    int next;
}edge [5008];
int head[1014];
int cnt;
int n,m;
double weigth[1012];
void init()
{
    cnt=0;
    for(int i=1;i<=n;i++)
    head[i]=-1;
}
void insert(int u,int v,int len)
{
    edge[cnt].u=v;
    edge[cnt].len=len;
    edge[cnt].next=head[u];
    head[u]=cnt;
    cnt++;
}


bool spfa(double mid)
{
     double dist[1012];
     int visited[1012];
     int used[1016];
     int que[1002*1002];
     int i, v;
     double newdist;
    for (i = 1; i <= n; i++)
    {
     dist[i] = inf;
     visited[i] = false;
     used[i] = 0;
    }
 dist[1] = 0;
 int hea = 0, tail = 0;
 que[tail++] = 1;
 visited[1] = true;
 used[1]++;
 while (hea < tail)
 {
  int u = que[hea];
  visited[u] = false;
  hea++;
  for(i=head[u];i!=-1;i=edge[i].next)
  {
    v = edge[i].u;
   newdist = mid * edge[i].len - weigth[v];  //新的边权值
   if (dist[u] + newdist < dist[v])
   {
    dist[v] = dist[u] + newdist;
    if (!visited[v])
    {
     que[tail++] = v;
     visited[v] = true;
     used[v]++;
     if (used[v] >= n)           //有负权环路
      return false;
    }
   }
  }
 }
 return true;     //无负权环路
}


int main()
{
  while(cin >>n>>m)
  {
    init();
   for(int i=1;i<=n;i++)
   cin>>weigth[i];
   for(int i=0;i<m;i++)
   {
       int uu,vv,cc;
       cin>>uu>>vv>>cc;
       insert(uu,vv,cc);
   }
   double mid;
   double L=0;
   double R=10066;
   double ans=0;
   while(R-L>=0.001)//二分
   {
      mid=(L+R)/2;
      if(spfa(mid))
      {
         R=mid;
      }
      else
      {
         ans=mid;
         L=mid;
      }
   }
   printf("%.2lf\n",ans);
 }
    return 0;
}

POJ3259----Wormholes(SPFA判断负环)

判断给定的有向图中是否存在负环。       利用 spfa 算法判断负环有两种方法:       1) spfa 的 dfs 形式,判断条件是存在一点在一条路径上出现多次。    ...
  • hy1405430407
  • hy1405430407
  • 2016年04月07日 20:58
  • 1307

POJ 3259 Wormholes(SPFA算法判断是否存在负环)

Wormholes Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 36755   Accepted: 13457...
  • zwj1452267376
  • zwj1452267376
  • 2015年08月19日 21:03
  • 775

用SPFA判断负环的方法及其优化

判负环的优化
  • wyh0410
  • wyh0410
  • 2016年11月03日 15:40
  • 2829

poj 3621 二分+spfa判负环

http://poj.org/problem?id=3621
  • u012774187
  • u012774187
  • 2014年11月09日 15:36
  • 399

poj 3621(参数搜索 + spfa)

01整数规划问题就是求解方程(a1*x1+a2*x2+..+an*xn)/(b1*x1+b2*x2+..+bn*xn)的最小值/最大值问题。其中xi = 0或1(i=1,2...n) 对于此类问...
  • water_glass
  • water_glass
  • 2011年10月05日 16:27
  • 411

poj 3621 二分+spfa

题意:给出一个有向图,问求一个回路,使得回路上的点权之和/边权之和最大。 这题主要是分析出如何确定ans值。我们将(a1*x1+a2*x2+..+an*xn)/(b1*x1+b2*x2+..+b...
  • amourjun
  • amourjun
  • 2013年07月29日 17:23
  • 306

SPFA判负环|BFS|DFS

判断给定的有向图中是否存在负环。       利用 spfa 算法判断负环有两种方法:       1) spfa 的 dfs 形式,判断条件是存在一点在一条路径上出现多次。       2)...
  • ws_yzy
  • ws_yzy
  • 2015年11月05日 07:33
  • 3450

easy sssp(spfa判断负环)

easy sssp描述: 输入数据给出一个有N(2
  • cax1165
  • cax1165
  • 2016年08月05日 17:57
  • 583

01参数搜索+二分+spfa判负环

【题目大意】 给出一个有向图,问求一个回路,使得回路上的点权之和/边权之和 最大。 【解题思路】 转:此题是对01分数规划的应用,那么首先明白01分数规划的思想. 01整数规划问题就是求解方程...
  • Weiguang_123
  • Weiguang_123
  • 2012年09月29日 15:39
  • 1024

Bellman_Ford和Spfa两种方法判断负环

Bellman_Ford是通过N-1次循环求出每个点到原点的最短路的,每次循环遍历所有的边,如果能更新就更新。前面说过为什么N-1次就行(若一个点到起点的最短路需要经过N个点,包括他自己,那么第N次循...
  • corncsd
  • corncsd
  • 2014年01月17日 17:09
  • 581
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:poj 3621(参数搜索+二分+spfa负环求解)
举报原因:
原因补充:

(最多只允许输入30个字)