单调栈应用[通俗易懂]

单调栈应用[通俗易懂]题意大概让算t秒每一秒的最短路和,每条边每秒dis++;dp算距离dis[i][j]:1到i点走j条边的最短路O(n(n+m));单调栈维护个上凸;等差数列n^2求解;代码:#include#include#include#include#includeusingnamespacest

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

题意大概让算t秒每一秒的最短路和,每条边每秒dis++;

dp算距离dis[i][j]:1到i点走j条边的最短路O(n(n+ m));
单调栈维护个上凸;
等差数列n^2求解;
代码

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<vector>
using namespace std;

long long  mod=1e9+7;
int head[5505],tot,to[15005],val[15005],nex[15005];
int dis[5505][5505];
int n,m,t;
int inf=99999999; 
long long ans;
void add(int x,int y,int v)
{
    int tmp=head[x];
    head[x]=++tot;
    to[tot]=y;
    val[tot]=v;
    nex[tot]=tmp;   
}
void cal()
{
    for(int i=1;i<=n;i++)
    for(int j=0;j<=n;j++)
    dis[i][j]=inf;
    dis[1][0]=0;
    for(int i=0;i<n-1;i++)
    for(int j=1;j<=n;j++)
    if(dis[j][i]<inf)
    {
        for(int now=head[j];now!=-1;now=nex[now])
        {
                int next=to[now];

                if(dis[next][i+1]==inf||dis[next][i+1]>dis[j][i]+val[now])
                { 
                    dis[next][i+1]=dis[j][i]+val[now];  
                }
        }

    }
}
int q[5005],from[5005];
long long  C(int k0,int a0,int k1,int a1,int k2,int a2)   //locate x cordinate
{
    return (long long)(a0-a1)*(k2-k0)-(long long)(a0-a2)*(k1-k0);
}
void solve(int x)
{
    tot=0;
    for(int i=n-1;i>=0;i--)
    {

        if(dis[x][i]==inf) continue;
        while(tot>=2&&C(q[tot-1],dis[x][q[tot-1]],q[tot],dis[x][q[tot]],i,dis[x][i])>=0)    tot--;
        q[++tot]=i; 

    }
    from[1]=0;
    for(int i=2;i<=tot;i++)
    from[i]=(dis[x][q[i]]-dis[x][q[i-1]]+q[i-1]-q[i]-1)/(q[i-1]-q[i]);  //trans to int(up)a
    int to=t+1; 
    for(int i=1;i<=tot;i++) from[i]=max(0,from[i]);
    for(int i=tot;i>=1;i--)
    {
        if(from[i]>=to) continue;
        ans=(ans+(long long)(dis[x][q[i]])*(to-from[i])%mod)%mod;
        ans=(ans+(long long )(to+from[i]-1)* ((to-from[i])/2)%mod*q[i])%mod;
        to=from[i]; 
    }
}
int main()
{
    int aa,bb,cc;
    scanf("%d%d%d",&n,&m,&t);
    memset(head,-1,sizeof(head));
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&aa,&bb,&cc);
        add(aa,bb,cc);
        add(bb,aa,cc);
    }
    cal();
    for(int i=1;i<=n;i++)
    solve(i);
    cout<<(ans+mod)%mod<<endl;

}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/190778.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号