bzoj2259：[Oibh]新型计算机（最短路）

2259：[Oibh]新型计算机

时间限制: 1000 ms 内存限制: 131072 KB

题目描述

Tim正在摆弄着他设计的“计算机”，他认为这台计算机原理很独特，因此利用它可以解决许多难题。

但是，有一个难题他却解决不了，是这台计算机的输入问题。新型计算机的输入也很独特，假设输入序列中有一些数字（都是自然数——自然数包括0），计算机先读取第一个数字S1，然后顺序向后读入S1个数字。接着再读一个数字S2，顺序向后读入S2个数字……依此类推。不过只有计算机正好将输入序列中的数字读完，它才能正确处理数据，否则计算机就会进行自毁性操作！

Tim现在有一串输入序列。但可能不是合法的，也就是可能会对计算机造成破坏。于是他想对序列中的每一个数字做一些更改，加上一个数或者减去一个数，当然，仍然保持其为自然数。使得更改后的序列为一个新型计算机可以接受的合法序列。

不过Tim还希望更改的总代价最小，所谓总代价，就是对序列中每一个数操作的参数的绝对值之和。

写一个程序：
从文件中读入原始的输入序列；
计算将输入序列改变为合法序列需要的最小代价；
向输出文件打印结果。

输入

输入文件包含两行，第一行一个正整数N，N<1 000 001。
输入文件第二行包含N个自然数，表示输入序列。

输出

仅一个整数，表示把输入序列改变为合法序列需要的最小代价，保证最小代价小于$10^9$。

输入样例

4
2 2 2 2

输出样例

1

解：

这道题乍一眼看是一道dp题。
当然也可以使用dp，然而我并不是用dp做的。
我们发现一个奇妙的性质。我们从i号点可以跳到i+a[i]+1，并不付出任何代价，我们走到旁边的点就需要多付出1的代价。然而只有在我们跳了一下之后才能走1权边。所以不是所有相邻点都可以连边。
序列的某些贪心dp问题都可以用最短路做，那么二维的可不可以呢？我觉得可以思考。
code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
struct lxy{
    int next,to,len;
}eg[3000005];
struct lxy1{
    int id,num;
    bool operator < (const lxy1 &QAQ)const{
      return num>QAQ.num;
    }
}yyy;

int head[1000005],dis[1000005],n,l,r,x,cnt;
bool vis[1000005];
priority_queue <lxy1> d;

int readit()
{
    int a=0;char s=getchar();
    while(s<'0'||s>'9') s=getchar();
    while(s>='0'&&s<='9') a=a*10+s-'0',s=getchar();
    return a;
}

void add(int op,int ed,int len){
    eg[++cnt].next=head[op];
    eg[cnt].to=ed;
    eg[cnt].len=len;
    head[op]=cnt;
    //cout<<op<<" "<<ed<<" "<<len<<endl;
}

void spfa()
{
    memset(dis,0x7f,sizeof(dis));
    dis[1]=0;
    yyy.id=1,yyy.num=0,d.push(yyy);
    while(!d.empty())
    {
        lxy1 now=d.top();
        for(int i=head[now.id];i!=-1;i=eg[i].next)
          if(dis[now.id]+eg[i].len<dis[eg[i].to])
          {
            dis[eg[i].to]=dis[now.id]+eg[i].len;
            yyy.id=eg[i].to;yyy.num=dis[eg[i].to];
            d.push(yyy);
          }
        vis[now.id]=1;
        while(!d.empty()&&vis[d.top().id]==1)d.pop();
    }
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d",&n);
    r=n,l=n;
    for(int i=1;i<=n;i++)
    {
        x=readit();
        if(x+i+1<=n) add(i,x+i+1,0),r=min(r,x+i+1),l=min(l,i);
        else add(i,n+1,i+x-n);
    }
    for(int i=l;i<=n;i++)
      add(i+1,i,1);
    for(int i=r;i<=n;i++)
      add(i,i+1,1);
    spfa();
    printf("%d",dis[n+1]);
}