状态转移有两种,一种是对于一个状态v找出这个状态能够来自哪些状态u,然后更新v,但有时这种转移方法不太好用,那么考虑另一种转移方法,对于一个状态u,我们考虑状态u能够影响的所有状态v,然后不断松弛,像spfa一样...
这题就是第二种转移方法。
另外Inq标记数组确实能节省时间,毕竟减少了节点入队次数... 这题不加inq16ms,加上inq 0ms
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<fstream>
using namespace std;
#define INF 0x7fffffff
int DP[15][1<<10];
bool Inq[15][1<<10];
struct Road
{
int b,c,p,r;
Road()
{}
Road(int b,int c,int p,int r)
{
this->b=b,this->c=c,this->p=p,this->r=r;
}
};
vector<Road> Gra[15];
int work(int n)
{
queue<int> q1,q2;
int i,j,u,v,c,p,r,s;
memset(Inq,0,sizeof(Inq));
for(i=0;i<10;++i)
for(j=0;j<(1<<10);++j)
DP[i][j]=INF;
DP[0][1]=0,Inq[0][1]=true;
q1.push(0),q2.push(1);
while(!q1.empty())
{
u=q1.front();
q1.pop();
s=q2.front();
q2.pop();
Inq[u][s]=false;
for(i=0;i<Gra[u].size();++i)
{
v=Gra[u][i].b,c=Gra[u][i].c,p=Gra[u][i].p,r=Gra[u][i].r;
if(s&(1<<c))
{
if(DP[u][s]+p<DP[v][s|(1<<v)])
{
DP[v][s|(1<<v)]=DP[u][s]+p;
if(!Inq[v][s|(1<<v)])
{
q1.push(v);
q2.push(s|(1<<v));
Inq[v][s|(1<<v)]=true;
}
}
}
else
{
if(DP[u][s]+r<DP[v][s|(1<<v)])
{
DP[v][s|(1<<v)]=DP[u][s]+r;
if(!Inq[v][s|(1<<v)])
{
q1.push(v);
q2.push(s|(1<<v));
Inq[v][s|(1<<v)]=true;
}
}
}
}
}
int ans=INF;
for(i=0;i<(1<<n);++i)
{
if((i&1)&&(i&(1<<(n-1))))
ans=min(ans,DP[n-1][i]);
}
return ans;
}
int main()
{
int n,m,i,j,a,b,c,p,r,ans;
scanf("%d %d",&n,&m);
while(m--)
{
scanf("%d %d %d %d %d",&a,&b,&c,&p,&r);
a--,b--,c--;
Gra[a].push_back(Road(b,c,p,r));
}
ans=work(n);
if(ans==INF)
printf("impossible\n");
else
printf("%d\n",ans);
return 0;
}