首先 化简一下式子;
log2[(1+a1)/1]+log2[(1+a1+a2)/(1+a1)]+log2[(1+a1+a2+a3)/(1+a1+a2)]+....+log2[(1+a1+a2+...+an)/(1+a1+a2+...+an-1)]=log2(1+a1+a2+a3+a4+a5+...an);
直接维护 1+a1+a2+...+an的 最小值 ,相当于 最短路;
其次 log2()手写
AC代码
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#include <math.h>
using namespace std;
typedef unsigned long long ll;
const int maxn=1e5+10;
const ll inf=1e19;
struct point{
int to;
int next;
int val;
int val1;
}pt[3*maxn];
struct node{
int id;
ll vv;
ll level;
bool operator < (const struct node &a ) const
{
return vv > a.vv;
}
};
int q,head[3*maxn],vis[maxn];
ll dis[maxn];
void init()
{
q=0;
memset(head,-1,sizeof(head));
}
void add(int u,int v,int val,int val1)
{
pt[q].next=head[u];
pt[q].to=v;
pt[q].val=val;
pt[q].val1=val1;
head[u]=q++;
}
void dij(int st)
{
priority_queue<node>q;
for (int i=0;i<maxn;i++)
dis[i]=inf;
memset(vis,0,sizeof(vis));
while (!q.empty())
q.pop();
node no;
dis[st]=1;
no.id=st;
no.level=1;
no.vv=dis[st];
q.push(no);
while (!q.empty())
{
node nn=q.top();
q.pop();
if (vis[nn.id])
continue;
vis[nn.id]=1;
node noo;
for (int i=head[nn.id];i!=-1;i=pt[i].next)
{
int v=pt[i].to;
if (!vis[v] && ((ll)pt[i].val)/nn.level >= (1ll<<pt[i].val1)-1 && dis[nn.id]+pt[i].val<dis[v])
{
dis[v]=dis[nn.id]+pt[i].val;
noo.id=v;
noo.level=pt[i].val+nn.level;
noo.vv=dis[nn.id]+pt[i].val;
q.push(noo);
}
}
}
}
ll log2(ll n)
{
for(int i=1;;i++)
if((1ll<<i)>n)
return i-1;
}
int main()
{
int n,m,t,u,v,a,b;
scanf("%d",&t);
while (t--)
{
init();
scanf("%d%d",&n,&m);
for (int i=0;i<m;i++)
{
scanf("%d%d%d%d",&u,&v,&a,&b);
add(u,v,a,b);
}
dij(1);
if (dis[n]==inf)
printf("-1\n");
else
{
//cout<<dis[n]<<endl;
ll vvv=log2(dis[n]);
printf("%lld\n",vvv);
}
}
return 0;
}