2019年7月训练(壹)

2019-07-25

luogu P3627 [APIO2009]抢掠计划

卡了三个小时,看了题解才作出来的(菜)

前驱知识:

壹~邻接表存储/遍历

贰~SPFA跑最长路(<改>就行了)

叄~Tarjan缩点



壹.邻接表储存

两个,add存无边权,未缩点;build有边权,已缩点。

void add(int u,int v)
 {
    cnt++;
    e[cnt].to=v;
    e[cnt].next=head[u];
    head[u]=cnt;
 }

  

 void build(int u,int v,int w)
 {
    cnt++;
    e[cnt].to=v;
    e[cnt].val=w;
    e[cnt].next=head[u];
    head[u]=cnt;
 }


贰.输入

懒得开两个head数组了,所以memset了。

cnt记得重置为零。

   cnt=0;
    memset(e,0,sizeof(e));
    memset(head,0,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&u[i],&v[i]);
        add(u[i],v[i]);
    }
    for(int i=1;i<=n;i++) scanf("%d",&w[i]);
    scanf("%d%d",&s,&p);
    for(int i=1;i<=p;i++) scanf("%d",&bar[i]);


叄.tarjan 缩点

 void Tarjan(int x)
 {
    dfn[x]=low[x]=++total;
    stk[++top]=x;vis[x]=true;
    for(int i=head[x];i;i=e[i].next)
    {
        int t=e[i].to;
        if(!dfn[t])
        {
            Tarjan(t);
            low[x]=min(low[x],low[t]);
        }
        else if(vis[t])
         low[x]=min(low[x],dfn[t]);
    }
    if(dfn[x]==low[x])
    {
        tot++;
        do{
            int tp=stk[top];
            sum[tot]+=w[tp];
            vis[tp]=false;
            be[tp]=tot;
        }while(stk[top--]!=x);
    }//数组模拟栈
 }

 

这里我是用的数组模拟栈,首先int tp=stk[top]取出栈顶

sum表示缩完点后这个点的点权
不是很懂Tarjan的好好理解下Tarjan弹栈的部分再看这里就懂了
懂Tarjan的模拟一下应该就懂了,每次弹栈时,所有被弹出的点都是缩完点后的一个点
即sum[tot]+=w[tp],缩完点后的点权+=原点权
很好理解吧
不懂的说明你对Tarjan还是理解不到位……这篇题解不是讲Tarjan的,楼下大佬应该有详细讲解。

然后vis[tp]=false表示tp已经出栈

g表示缩完点后每个点在哪个点中
即g[tp]=tot,tp这个点在缩完点后的第tot个点里
然后用栈顶和x比较,标准Tarjan操作

stk[top--]就相当于pop弹栈了

肆.Spfa跑最长路
 void Spfa(int s)
 {
    for(int i=1;i<=tot;i++) dis[i]=0;
    int bes=be[s];//不想认真想变量名了
    q.push(bes);
    vis[bes]=true;
    dis[bes]=sum[bes];
    while(!q.empty())
    {
        int h=q.front();q.pop();
        vis[h]=false;
        for(int i=head[h];i;i=e[i].next)
        {
            int t=e[i].to;
            if(dis[t]<dis[h]+e[i].val)
            {
                dis[t]=dis[h]+e[i].val;
                if(!vis[t])
                {
                    q.push(t);
                    vis[t]=true;
                }
            }
        }
    }
 }

起点边权解决方法:

起点从s变为be[s]了,dis[s]不是0了

因为缩点

我们刚刚build建的图是一张缩完点后的缩点图,所以我们push起点的时候当然是push缩完点后s所在的点啊

为什么呢?万一s本身就在一个环里,push(s)的话问题就大了

所以我们把起点一律改成g[s]来操作

并且,我们把起点的dis值直接加上点权,这样就不会漏掉起点点权了

而且放心,一开始更改dis值不会对后续松弛操作造成无法更改的影响,毕竟这是最长路



伍.main函数部分

 for(int i=1;i<=m;i++)
     if(be[u[i]]!=be[v[i]])
      build(be[u[i]],be[v[i]],sum[be[v[i]]]);

这时候前面的u,v数组便起作用了

我们通过判断u[i],v[i]是否在一个新点里,如果不在的话就build一条新边,这样就能建一个缩完点后的图

!!!build:起点是u[i]所在的缩点,终点是v[i]所在的缩点(如果u[i]和v[i]在一个缩点里就不会执行build了),边权是v[i]所在缩点的点权(sum)!!!

这样,这张新图就表明:从U走到V可以抢劫W的钱,钱数当然是缩点点权啊

但是-------------------------------------------------------------------------这样每次把边权设置为后头那个缩点的点权,前头那个点(起点)被忽略了,点权不就没用了吗?

对√的确有这个问题

但是在上一部分我们已经简单粗暴の解决了



以下为AC码:

#include<cstdio>
#include<cstring>
#include<stack>
#include<algorithm>
#include<queue>
#define maxn 500010
using namespace std;
struct edbee
{
    int to,val,next;
} e[maxn];
int m,n,p,s,cnt,be[maxn],u[maxn],v[maxn],w[maxn],head[maxn],bar[maxn],dis[maxn],dfn[maxn],low[maxn],stk[maxn],sum[maxn];
bool vis[maxn];
queue <int> q;
int ans,top,tot,total;

 void add(int u,int v)
 {
    cnt++;
    e[cnt].to=v;
    e[cnt].next=head[u];
    head[u]=cnt;
 }

 void build(int u,int v,int w)
 {
    cnt++;
    e[cnt].to=v;
    e[cnt].val=w;
    e[cnt].next=head[u];
    head[u]=cnt;
 }

 void Tarjan(int x)
 {
    dfn[x]=low[x]=++total;
    stk[++top]=x;vis[x]=true;
    for(int i=head[x];i;i=e[i].next)
    {
        int t=e[i].to;
        if(!dfn[t])
        {
            Tarjan(t);
            low[x]=min(low[x],low[t]);
        }
        else if(vis[t])
         low[x]=min(low[x],dfn[t]);
    }
    if(dfn[x]==low[x])
    {
        tot++;
        do{
            int tp=stk[top];
            sum[tot]+=w[tp];
            vis[tp]=false;
            be[tp]=tot;
        }while(stk[top--]!=x);
    }
 }

 void Spfa(int s)
 {
    for(int i=1;i<=tot;i++) dis[i]=0;
    int bes=be[s];
    q.push(bes);
    vis[bes]=true;
    dis[bes]=sum[bes];
    while(!q.empty())
    {
        int h=q.front();q.pop();
        vis[h]=false;
        for(int i=head[h];i;i=e[i].next)
        {
            int t=e[i].to;
            if(dis[t]<dis[h]+e[i].val)
            {
                dis[t]=dis[h]+e[i].val;
                if(!vis[t])
                {
                    q.push(t);
                    vis[t]=true;
                }
            }
        }
    }
 }

int main()
{
    cnt=0;
    memset(e,0,sizeof(e));
    memset(head,0,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&u[i],&v[i]);
        add(u[i],v[i]);
    }
    for(int i=1;i<=n;i++) scanf("%d",&w[i]);
    scanf("%d%d",&s,&p);
    for(int i=1;i<=p;i++) scanf("%d",&bar[i]);
    for(int i=1;i<=n;i++)
     if(!dfn[i]) Tarjan(i);
    cnt=0;
    memset(e,0,sizeof(e));
    memset(head,0,sizeof(head));
    for(int i=1;i<=m;i++)
     if(be[u[i]]!=be[v[i]])
      build(be[u[i]],be[v[i]],sum[be[v[i]]]);
    Spfa(s);
    for(int i=1;i<=p;i++)
     ans=max(ans,dis[be[bar[i]]]);
    printf("%d",ans);
    return 0;
}

谢谢你能看到这里。

2019-07-25 22:35:43

以下为私货,可以不看

REOL要出新专啦!!!

唱歌超好听,人长得超好看的!!!

国内Apple music:https://music.apple.com/cn/album/1472096049?app=music

再次谢谢你能看到这里。

2019-07-2522:36:16

原文地址:https://www.cnblogs.com/plzplz/p/11247078.html

时间: 2024-11-09 09:50:43

2019年7月训练(壹)的相关文章

2019年9月训练(壹)数位DP (HDU 2089)

开学之后完全没时间写博客.... HDU 2089 不要62(vjudge) 数位DP 思路: 题目给出区间[n,m] ,找出不含4或62的数的个数 用一个简单的差分:先求0~m+1的个数,再减去0~n的个数. 但问题依旧不简单,再次简化为求0~i位数中不含4或62的数的个数. i=1 //0~9中 i=2 //0~99中 i=3 //0~999中 ...... dp[i][0] //0~i位数中的吉利数 dp[i][1] //0~i位数中以2打头的吉利数 dp[i][2] //0~i位数中的非

2019年8月训练(壹)二分,三分

二分查找 P1024 一元三次方程求解 题目给出范围[-100,100],同时两根绝对值之差<=1,保证了每一个大小为1的区间里至多有1个解,也就是说当区间的两个端点的函数值异号时区间内一定有一个解,同号时一定没有解. 也就可以二分去做查找. AC码: #include<cstdio> #include<cstring> #include<algorithm> using namespace std; double a,b,c,d; double f(double

2019年1月训练记录(更新ing)

前言 时光飞逝,转眼间,便到了\(2019\)年. 这一年里,还是要继续努力吧,好好学习一些新的算法. 当然,还是要好好准备一下期末考试啦... ... \(Jan\ 1st\) 原文地址:https://www.cnblogs.com/chenxiaoran666/p/2019Jan.html

2019年7月训练(陆)

模板:luogo P3379 [模板]最近公共祖先(LCA) 今天讲的时候有点跑神,现在卑微地来补习(菜) LCA指的是最近公共祖先(Least Common Ancestors). 最简单的算法无疑是从两个点一个个往上走,出现的第一个两个点都走过的点即为两点的LCA. 但是时间很长. 所以起用倍增,倍增的作用就是将两点上升所需的复杂度减低. 大致流程为:将deep不同的两个点跳到同一层,再跳到deep[lca-1]的那层,再向上跳一层就是lca了. 加速跳的方法就是每次向上跳的层数为2的i次方

周记 - 2019年11月03日

2019年11月05日 2019年徐州区域赛结束了.封榜前3题铜牌前部,封榜后最后27分钟罚了4次通过E题.虽然4题罚时爆炸,不过万幸得了银牌后部.目前看应该还会再参加一年的,这个博客会不断更新记录最后一年参加比赛的学习进度(以及最后两年本科的其他事情),今年的目标是做一个真正的全能选手,首先希望在寒假结束前把Codeforces打到橙色(2100+),证明自己思维还行吧. 既然想做全能选手,每个专题都要学到省选级别吧,学习的分类就参照OI-Wiki的分类,题单的话,模板题从洛谷找,思维题从Co

蔡康永的说话之道——2019年12月15日

.bodyContainer { font-family: Arial, Helvetica, sans-serif; text-align: center; padding-left: 32px; padding-right: 32px; } .notebookFor { font-size: 18px; font-weight: 700; text-align: center; color: rgb(119, 119, 119); margin: 24px 0px 0px; padding:

西安活动 | 2019年1月13号 &quot;拥抱开源, 又见.NET&quot; 线下交流活动报名进行中

随着.NET Core的发布和开源,.NET又重新回到人们的视野..NET Core的下个3.0即将release,加入非常多的新功能,越来越拥抱变化,DevOps和Microservice的最佳实践已经在.NET Core落地,比如 Ocelot网关.Grpc+Consul 服务注册发现.Apworks CQRS实现.Xigadee 微服务工具库.脚手架. 西安.NET社区组织发起了此次“拥抱变化, 又见.NET”线下交流活动,邀请了三位资深.NET开发者作为分享讲师,他们将从架构.原理.语言

2019年2月26日【整理物品,下载收集考研资料,明天正式开始复习】

2019年1月26日星期六 一:一句话木马重学习 1.网站安全狗网马查杀 http://download.safedog.cn/download/software/safedogwzApache.exe 2.D盾 Web 查杀 http://www.d99net.net/down/WebShellKill_V2.0.9.zip 3 深信服WebShellKillerTool http://edr.sangfor.com.cn/tool/WebShellKillerTool.zip 4 BugSc

java 手机号正则表达式 2019年1月

  import java.util.regex.Matcher; import java.util.regex.Pattern; import org.apache.commons.lang3.StringUtils; /**  * @author wucai  *三大运营商号码均可验证(不含卫星通信1349)  */ public class mobile {     /*  <br> 2019年1月16日已知     中国电信号段         133,149,153,173,174,