hdu 4435 charge-station

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

有题意可以知道 对应第k个城市 即使前(k-1)个城市都建立加油站 也没有在第k个城市建立加油站花的钱多

思路 ：

如果全部都建立加油站 还是不能满足条件的话 则说明没有解决方案

否则 从后面逐个向前尝试看是否可以不建立

对于 是否能达到目的 我用的方法比较麻烦

两次 spfa 分别求到某个城市时 还未在此站加油（也可能没有加油站）最多可以有多少油 和 在某个城市想回来的话最少需要的油量

每个点必须满足 有的油量可以回来

代码：

#include<iostream>

#include<cstdio>

#include<cstring>

#include<cmath>

#include<algorithm>

#include<vector>

#include<set>

#include<queue>

#include<stack>

#include<map>

#include<string>

#include<iomanip>

using namespace std;



#define LL long long

const int INF=0x5fffffff;

const double FINF=1e9;

const int N=205;

int d[N][N];

int x[N];

int y[N];

int dist[N];

int need[N];

int had[N];

bool in[N];

void spfa1(int s,int n,int k)

{

    queue<int>qt;

    memset(in,false,sizeof(in));

    for(int i=1;i<=n;++i)

    dist[i]=-1;

    dist[s]=0;

    in[s]=true;

    qt.push(s);

    while(!qt.empty())

    {

        int x=qt.front();

        qt.pop();

        in[x]=false;

        for(int i=1;i<=n;++i)

        {

            if(i==s||i==x||d[x][i]>min(dist[x]+had[x],k))

            continue;

            int tmp=min(dist[x]+had[x],k)-d[x][i];

            if(dist[i]<tmp)

            {

                dist[i]=tmp;

                if(!in[i])

                {

                    in[i]=true;

                    qt.push(i);

                }

            }

        }

    }

}

void spfa2(int s,int n,int k)

{

    queue<int>qt;

    memset(in,false,sizeof(in));

    for(int i=1;i<=n;++i)

    need[i]=INF;

    need[s]=0;

    in[s]=true;

    qt.push(s);

    while(!qt.empty())

    {

        int x=qt.front();

        qt.pop();

        in[x]=false;

        for(int i=1;i<=n;++i)

        {

            if(d[x][i]+need[x]>k)

            continue;

            int tmp=need[x]+d[x][i];

            tmp=tmp-had[i];

            if(tmp<0)

            tmp=0;

            if(need[i]>tmp)

            {

                need[i]=tmp;

                if(!in[i])

                {

                    in[i]=true;

                    qt.push(i);

                }

            }

        }

    }

}

bool Can(int n,int k)

{

    spfa1(1,n,k);

    spfa2(1,n,k);

    for(int i=1;i<=n;++i)

    {

        if(need[i]>dist[i])

        return false;

    }

    return true;

}

int main()

{

    //freopen("data.txt","r",stdin);

    int n,k;

    while(scanf("%d %d",&n,&k)!=EOF)

    {

        for(int i=1;i<=n;++i)

        {

            scanf("%d %d",&x[i],&y[i]);

            for(int j=1;j<i;++j)

            {

                d[i][j]=d[j][i]=ceil(sqrt(1.0*(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])));

            }

        }

        for(int i=1;i<=n;++i)

        had[i]=k;

        if(!Can(n,k))

        {printf("-1\n");continue;}

        int ans[N];

        memset(ans,0,sizeof(ans));

        ans[1]=1;

        for(int i=n;i>1;--i)

        {

            had[i]=0;

            if(!Can(n,k))

            {had[i]=k;ans[i]=1;}

        }

        int i=n;

        while(!ans[i])

        --i;

        for(int l=i;l>=1;--l)

        printf("%d",ans[l]);

        printf("\n");

    }

    return 0;

}