HDU 3339 In Action 最短路+01背包

http://acm.hdu.edu.cn/showproblem.php?pid=3339

 

Since 1945, when the first nuclear bomb was exploded by the Manhattan Project team in the US, the number of nuclear weapons have soared across the globe.
Nowadays,the crazy boy in FZU named AekdyCoin possesses some nuclear weapons and wanna destroy our world. Fortunately, our mysterious spy-net has gotten his plan. Now, we need to stop it.
But the arduous task is obviously not easy. First of all, we know that the operating system of the nuclear weapon consists of some connected electric stations, which forms a huge and complex electric network. Every electric station has its power value. To start the nuclear weapon, it must cost half of the electric network's power. So first of all, we need to make more than half of the power diasbled. Our tanks are ready for our action in the base(ID is 0), and we must drive them on the road. As for a electric station, we control them if and only if our tanks stop there. 1 unit distance costs 1 unit oil. And we have enough tanks to use.
Now our commander wants to know the minimal oil cost in this action.

Input

The first line of the input contains a single integer T, specifying the number of testcase in the file.
For each case, first line is the integer n(1<= n<= 100), m(1<= m<= 10000), specifying the number of the stations(the IDs are 1,2,3...n), and the number of the roads between the station(bi-direction).
Then m lines follow, each line is interger st(0<= st<= n), ed(0<= ed<= n), dis(0<= dis<= 100), specifying the start point, end point, and the distance between.
Then n lines follow, each line is a interger pow(1<= pow<= 100), specifying the electric station's power by ID order.

Output

The minimal oil cost in this action.
If not exist print "impossible"(without quotes).

Sample Input

2
2 3
0 2 9
2 1 3
1 0 2
1
3
2 1
2 1 3
1
3

Sample Output

5
impossible

题目大意:给n个点和m条边,每个点有一个能量站,能量站有ai的能量,你有无数辆坦克可以从点0出发,你需要摧毁一些能量站使得这些能量站的总和大于等于总能量的一半,但是你的坦克每移动1m就要消耗1L升油,让你求出可以完成要求的最少油量,若不可能则输出impossible。

思路:最短路+01背包。这道题的数据有重边!!!让我找bug找了好久。(心态大崩)背包应该是有两种方法的:(1)把距离当做背包容量,每个点的能量当做价值,做01背包,同时记录价值大于等于能量一半的情况下所花费的最小容量。(2)把能量当做背包容量,到每个点的距离当做价值,做01背包,同时记录消耗背包体积大于等于总体积一半情况下所获得的最小价值。提供四种方法吧。

floyd+第一种背包方法:

#include<iostream>
#include<cstdio>
#include<stack>
#include<cmath>
#include<cstring>
#include<queue>
#include<set>
#include<algorithm>
#include<iterator>
#define INF 0x3f3f3f3f
typedef long long ll;
typedef unsigned long long ull;
using namespace std;

int s[105][105];
int a[105];
int dp[100005];
int n,m;

void floyd()
{
    for(int k=0;k<=n;k++)
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++)
                s[i][j]=min(s[i][j],s[i][k]+s[k][j]);
}
int main()
{
    int t;
    scanf("%d",&t);
	while(t--)
	{
		int b,e,v;
		memset(s,INF,sizeof(s));
		for(int i=0;i<=n;i++)
            s[i][i]=0;
		scanf("%d %d",&n,&m);
		for(int i=0;i<m;i++)
		{
			scanf("%d %d %d",&b,&e,&v);
            if(s[b][e]>v)
                s[b][e]=s[e][b]=v;
		}
		floyd();
		int sum=0,MAX=0;
		for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            sum+=a[i];
            if(s[0][i]!=INF)
                MAX+=s[0][i];
        }
        memset(dp,0,sizeof(dp));
        sum/=2;
        int MIN=INF;
        for(int i=1;i<=n;i++)
        {
            for(int j=MAX;j>=s[0][i];j--)
            {
                dp[j]=max(dp[j],dp[j-s[0][i]]+a[i]);
                if(dp[j]>sum)
                    MIN=min(MIN,j);
            }
        }
        if(MIN==INF)
            printf("impossible\n");
        else
            printf("%d\n",MIN);
	}
	return 0;
}

floyd+ 第二种背包方法:

#include<iostream>
#include<cstdio>
#include<stack>
#include<cmath>
#include<cstring>
#include<queue>
#include<set>
#include<algorithm>
#include<iterator>
#define INF 0x3f3f3f3f
typedef long long ll;
typedef unsigned long long ull;
using namespace std;

int s[105][105];
int a[105];
int dp[100005];
int n,m;

void floyd()
{
    for(int k=0;k<=n;k++)
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++)
                s[i][j]=min(s[i][j],s[i][k]+s[k][j]);
}
int main()
{
    int t;
    scanf("%d",&t);
	while(t--)
	{
		int b,e,v;
		memset(s,INF,sizeof(s));
		for(int i=0;i<=n;i++)
            s[i][i]=0;
		scanf("%d %d",&n,&m);
		for(int i=0;i<m;i++)
		{
			scanf("%d %d %d",&b,&e,&v);
            if(s[b][e]>v)
                s[b][e]=s[e][b]=v;
		}
		floyd();
		int sum=0;
		for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            sum+=a[i];
        }
        memset(dp,INF,sizeof(dp));
        dp[0]=0;
        int MIN=INF;
        for(int i=1;i<=n;i++)
        {
            for(int j=sum;j>=a[i];j--)
            {
                dp[j]=min(dp[j],dp[j-a[i]]+s[0][i]);
                if(j>sum/2)
                    MIN=min(MIN,dp[j]);
            }
        }
        if(MIN==INF)
            printf("impossible\n");
        else
            printf("%d\n",MIN);
	}
	return 0;
}

dijkstra+ 第一种背包方法:

#include<iostream>
#include<cstdio>
#include<stack>
#include<cmath>
#include<cstring>
#include<queue>
#include<set>
#include<algorithm>
#include<iterator>
#define INF 0x3f3f3f3f
typedef long long ll;
typedef unsigned long long ull;
using namespace std;

int s[105][105];
int d[105];
int a[105];
int vis[105];
int dp[100005];
int n,m;

void dijkstra()
{
    memset(vis,0,sizeof(vis));
    memset(d,INF,sizeof(d));
    d[0]=0;
    while(1)
    {
        int u=-1,MIN=INF;
        for(int i=0;i<=n;i++)
        {
            if(!vis[i]&&d[i]<MIN)
            {
                MIN=d[i];
                u=i;
            }
        }
        if(u==-1)
            break;
        vis[u]=1;
        for(int i=0;i<=n;i++)
            if(!vis[i]&&d[i]>d[u]+s[u][i])
                d[i]=d[u]+s[u][i];
    }
}

int main()
{
    int t;
    scanf("%d",&t);
	while(t--)
	{
		int b,e,v;
		memset(s,INF,sizeof(s));
		scanf("%d %d",&n,&m);
		edge temp;
		for(int i=0;i<m;i++)
		{
			scanf("%d %d %d",&b,&e,&v);
            if(s[b][e]>v)
                s[b][e]=s[e][b]=v;
        }
		dijkstra();
		int sum=0,MAX=0;
		for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            sum+=a[i];
            if(d[i]!=INF)
                MAX+=d[i];
        }
        memset(dp,0,sizeof(dp));
        sum/=2;
        int MIN=INF;
        for(int i=1;i<=n;i++)
        {
            for(int j=MAX;j>=d[i];j--)
            {
                dp[j]=max(dp[j],dp[j-d[i]]+a[i]);
                if(dp[j]>sum)
                    MIN=min(MIN,j);
            }
        }
        if(MIN==INF)
            printf("impossible\n");
        else
            printf("%d\n",MIN);
	}
	return 0;
}

堆优化dijkstra+ 第一种背包方法:

#include<iostream>
#include<cstdio>
#include<stack>
#include<cmath>
#include<cstring>
#include<queue>
#include<set>
#include<algorithm>
#include<iterator>
#define INF 0x3f3f3f3f
typedef long long ll;
typedef unsigned long long ull;
using namespace std;

struct edge
{
    int to,v;
};

struct node
{
    int dis;
    int id;
    bool operator <(const node&a)const
    {
        return dis>a.dis;
    }
};

vector<edge> vec[105];
int d[105];
int a[105];
int vis[105];
int dp[100005];
int n,m;

void dijkstra()
{
    memset(vis,0,sizeof(vis));
    memset(d,INF,sizeof(d));
    d[0]=0;
    node t;
    t.id=0,t.dis=0;
    priority_queue<node> q;
    q.push(t);
    while(!q.empty())
    {
        node h=q.top();
        q.pop();
        int id=h.id;
        if(d[id]<h.dis)
            continue;
        for(int i=0;i<vec[id].size();i++)
        {
            int to=vec[id][i].to;
            if(d[to]>d[id]+vec[id][i].v)
            {
                d[to]=d[id]+vec[id][i].v;
                node temp;
                temp.id=to;
                temp.dis=d[to];
                q.push(temp);
            }
        }
    }
}

int main()
{
    int t;
    scanf("%d",&t);
	while(t--)
	{
		int b,e,v;
		scanf("%d %d",&n,&m);
		for(int i=0;i<=n;i++)
            while(vec[i].size())
                vec[i].pop_back();
		edge temp;
		for(int i=0;i<m;i++)
		{
			scanf("%d %d %d",&b,&e,&v);
            temp.v=v;
            temp.to=e;
            vec[b].push_back(temp);
            temp.to=b;
            vec[e].push_back(temp);
        }
		dijkstra();
		int sum=0,MAX=0;
		for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            sum+=a[i];
            if(d[i]!=INF)
                MAX+=d[i];
        }
        memset(dp,0,sizeof(dp));
        sum/=2;
        int MIN=INF;
        for(int i=1;i<=n;i++)
        {
            for(int j=MAX;j>=d[i];j--)
            {
                dp[j]=max(dp[j],dp[j-d[i]]+a[i]);
                if(dp[j]>sum)
                    MIN=min(MIN,j);
            }
        }
        if(MIN==INF)
            printf("impossible\n");
        else
            printf("%d\n",MIN);
	}
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值