Problem【★★★】
边的权值为l/v+c,其中已知l,求那些点不可能成为1 to n最短路上的点。
v>0,c>=0,v与c是确定的数。
Model Building【★★★★】
spfa+单调栈+dfs。
Solution【★★★★】
dis[i][j]表示经过j条边到达i的只考虑l下的最短路,经过k条边到达n的最短路为dis[n][k]/v+k*c,以c的大小为x轴,1 to n的最短路为y轴建立平面直角坐标系(v不影响直线与y轴交点的相对位置,因此也不影响直线间交点的相对位置):
每条直线代表一个路径方案,只有当该直线在x∈[0,∞]的某个位置为所有直线的最小值时该直线所代表的最短路方案才是合法的。因此只要spfa处理出dis,按k从大到小,单调栈维护凸壳,最后dfs出所有可能的点即可。注意,dfs时需记忆化。
总时间复杂度O(nm)。
Code【★★★☆】
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
using namespace std;
const int maxn=2005,maxm=20005;
const long long inf=1e15;
int n,m,cnt,last[maxn],cnt1,sol[maxn];
long long dis[maxn+5][maxn+5];
bool inq[maxn+5][maxn+5],ans[maxn];
struct edge{
int v,next;
long long w;
}e[maxm];
struct node{
int u,k;
};
inline void add(int u,int v,long long w)
{
e[++cnt].v=v;
e[cnt].w=w;
e[cnt].next=last[u];
last[u]=cnt;
}
void spfa()
{
queue<node>q;
for(int i=0;i<=n+5;i++)
for(int j=0;j<=n+5;j++)
dis[i][j]=inf;
dis[1][0]=0;
q.push({1,0});
while(!q.empty())
{
node nw=q.front();
q.pop();
int u=nw.u,k=nw.k;
inq[u][k]=0;
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v;
long long w=e[i].w;
if(dis[v][k+1]>dis[u][k]+w)
{
dis[v][k+1]=dis[u][k]+w;
if(!inq[v][k+1])
{
inq[v][k+1]=1;
q.push({v,k+1});
}
}
}
}
}
void dfs(int u,int k)
{
ans[u]=1;
if(inq[u][k]||k<=0)
return ;
inq[u][k]=1;
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v;
long long w=e[i].w;
if(dis[v][k-1]+w==dis[u][k])
dfs(v,k-1);
}
}
double getpos(int k1,int k2)
{
return (double)(dis[n][k2]-dis[n][k1])/(k1-k2);
}
void solve()
{
spfa();
long long maxm=inf,minm=inf;
for(int i=n;i>=1;i--)
{
if(dis[n][i]==inf) continue;
while(cnt1)
{
if(cnt1 == 1 && getpos(sol[cnt1],i) < 0) cnt1--;
else if(cnt1 > 1 && getpos(sol[cnt1],sol[cnt1-1]) > getpos(sol[cnt1],i)) cnt1--;
else break;
}
sol[++cnt1]=i;
}
memset(inq,0,sizeof(inq));
for(int i=1;i<=cnt1;i++)
dfs(n,sol[i]);
ans[1]=1;
int cnt2=0;
for(int i=1;i<=n;i++) if(!ans[i]) cnt2++;
cout<<cnt2<<endl;
for(int i=1;i<=n;i++)
if(!ans[i])
cout<<i<<' ';
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1,u,v,w;i<=m;i++)
{
scanf("%d%d%d",&u,&v,&w);
add(u,v,w);
add(v,u,w);
}
solve();
return 0;
}