• 2018.07.10【省赛模拟】模拟B组【JSOI2013】吃货JYY


    #Description
    世界上一共有N个JYY愿意去的城市,分别从1编号到N。JYY选出了K个他一定要乘坐的航班。除此之外,还有M个JYY没有特别的偏好,可以乘坐也可以不乘坐的航班。
    一个航班我们用一个三元组(x,y,z)来表示,意义是这趟航班连接城市x和y,并且机票费用是z。每个航班都是往返的,所以JYY花费z的钱,既可以选择从x飞往y,也可以选择从y飞往x。
    南京的编号是1,现在JYY打算从南京出发,乘坐所有K个航班,并且最后回到南京,请你帮他求出最小的花费。

    #Input
    输入数据的第一行包含两个整数N和K;
    接下来K行,每行三个整数x,y,z描述必须乘坐的航班的信息,数据保证在这K个航班中,不会有两个不同的航班在同一对城市之间执飞;
    第K+2行包含一个整数M;
    接下来M行,每行三个整数x,y,z 描述可以乘坐也可以不乘坐的航班信息。

    #Output
    输出一行一个整数,表示最少的花费。数据保证一定存在满足JYY要求的旅行方案。

    #Sample Input
    6 3
    1 2 1000
    2 3 1000
    4 5 500
    2
    1 4 300
    3 5 300

    #Sample Output
    3100

    #Data Constraint
    对于10%的数据满足N≤4;
    对于30%的数据满足N≤ 7;
    对于额外30%的数据满足,JYY可以只通过必须乘坐的K个航班从南京出发到达任意一个城市;
    对于100%的数据满足2≤N≤13,0≤K≤78,2 ≤M ≤ 200,1 ≤x,y ≤N,1 ≤z ≤ 10^4。

    #Hint
    样例说明:一个可行的最佳方案为123541。
    机票所需的费用为1000+1000+300+500+300=3100。

    #题解
    题意就是在一个图中,有k条边必须经过,其他边任意经过,从1号点出发,最终回到1号点的最小费用。
    这题我第一眼——最小生成树?(太天真了)
    第二眼——类似于欧拉回路?(然而我不会)
    第三眼——难道是网络流?(拉倒吧)
    第四眼——dp?但是有后效性,弃疗吧!
    然后,我就不会做了。

    结果正解就是一个dp加上类似于欧拉回路的东东再乱搞即可。
    听起来简单,搞起来很好搞,但是想到真的是神仙才想得到吧。
    我这个凡人算了

    首先,我们先看看画出来的图,我们发现在这个图中会有一个欧拉图,然后,我们就考虑把别的边或点给丢掉,来构成一个欧拉图。
    但是这样似乎很难搞,那么我们逆向想想,不如我们就新建一个图,来把边、点给丢进去,变成一个欧拉图,这样似乎很好搞啦。
    于是呢,我们就需要一个状态来记录。考虑到n极小,那么就是一个状态压缩。
    但是,看到是欧拉图,那么在新的图中的点有两种情况——一是奇点,二是偶点。什么意思?也就是奇点有奇数的边与之相连,偶点就是有偶数的边与之相连。然后还有的点为加入新图,那么就有3种情况——
    0:未加入图
    1:加入图且为奇点
    2:加入图且为偶点
    于是现在,我们就看看如何把一个点加入新图中。
    这里写图片描述
    我们就要考虑吧now给连上边

    第一种情况:
    当有必连边连接,那么就必须连上
    这里写图片描述
    然后,我们同时要改变下这个新图的点的奇偶性,即可

    第二种情况:
    没有必连边时,我们就考虑连一条最近的边。
    这里写图片描述
    然后这个显然是可行的。可以以最小的代价把这个点也给连上。
    但是,问题来了。正确性能否保证呢?因为以下图的走法,我们会发现开始节点的奇偶性好像会改变。
    这里写图片描述
    看看这个走路的方案,好像开始点的奇偶性是奇性,然而,这是个不严谨的欧拉图,这个走路分方案其实是太过于早考虑了,那么我们就先欠下这个账。

    先放下那个账,我们可以发现,我们新图中还剩下一些奇点是没有变成偶点的,根据欧拉图,我们就必须把这些奇点给改成偶点才可以行得通,那么我们就把奇点都给拿出来,然后再搞一遍dp,然后用最短路来把这些奇点以最小的费用给匹配出来。
    然后就可以吧所有的奇点给变成偶点了。
    正确?
    没错,因为每次都是两个点之间再连一条最短路的边,那么就代表这两个点之间返回用最小的费用。由于这是个不严谨的欧拉图,那么这个方法也是正确的。
    感性理解理解。

    欠的旧账——
    我们看看这个图:
    这里写图片描述
    用我们的方法来做,开始点的奇偶性为奇性,然后呢,我们就可以根据我们刚刚讲的方法来把这个开始点与另一个点匹配,那么也就可能走了一条重复的边,由于刚才说这样是正确的,那么这个走法也是正确的。
    也就是说,我们在加入点时,不考虑它走不走重复边,只有在最后把奇点匹配时才从真正意义上来走重复边。
    理解理解。
    然后你就会了。

    代码:

    uses math;
    var
            i,j,k,l,o,n,m,zs,answer,tot,tot1,maxx,zy,gs,add,p:longint;
            x,y,z,x1,y1,z1,tov,next,last,a,b:array[1..10000] of longint;
            map:array[1..13,1..13] of longint;
            f:array[0..1594323] of longint;
            pp:array[0..13] of longint;
            mi,mi2,apple:array[0..14] of longint;
            bz:array[1..14] of boolean;
            dp:array[0..8191] of longint;
            flag,flag1:boolean;
    procedure insert(x,y:longint);
    begin
            inc(tot);
            tov[tot]:=y;
            next[tot]:=last[x];
            last[x]:=tot;
    end;
    begin
            fillchar(map,sizeof(map),127);
            maxx:=map[1,1];
            readln(zs,n);
            for i:=1 to n do
            begin
                    readln(x[i],y[i],z[i]);
                    insert(x[i],y[i]);
                    insert(y[i],x[i]);
                    bz[x[i]]:=true;
                    bz[y[i]]:=true;
                    x[n+i]:=y[i];
                    y[n+i]:=x[i];
                    z[n+i]:=z[i];
                    map[x[i],y[i]]:=min(z[i],map[x[i],y[i]]);
                    map[y[i],x[i]]:=map[x[i],y[i]];
                    answer:=answer+z[i];
            end;
            readln(m);
            for i:=1 to m do
            begin
                    readln(x1[i],y1[i],z1[i]);
                    map[x1[i],y1[i]]:=min(map[x1[i],y1[i]],z1[i]);
                    map[y1[i],x1[i]]:=map[x1[i],y1[i]];
            end;
            for k:=1 to zs do
            begin
                    for i:=1 to zs do
                    begin
                            if i<>k then
                            for j:=1 to zs do
                            begin
                                    if (i<>j) and (k<>j) then
                                    begin
                                            if (map[i,k]<maxx) and (map[k,j]<maxx) then
                                            map[i,j]:=min(map[i,j],map[i,k]+map[k,j]);
                                    end;
                            end;
                    end;
            end;
            mi[1]:=1;
            mi2[1]:=1;
            for i:=2 to zs+1 do
            begin
                    mi[i]:=mi[i-1]*3;
                    mi2[i]:=mi2[i-1]*2;
            end;
            fillchar(f,sizeof(f),127);
            f[2]:=0;
            
            for i:=1 to mi[zs+1]-1 do
            begin
                    zy:=i;
                    flag1:=true;
                    for j:=1 to zs do
                    begin
                            apple[j]:=i mod mi[j+1] div mi[j];
                    end;
                    for j:=1 to zs do
                    begin
                            zy:=i;
                            if ((apple[j])=0) and (bz[j]) then
                            begin
                                    flag1:=false;
                                    flag:=true;
                                    gs:=0;
                                    k:=last[j];
                                    while k>0 do
                                    begin
                                            l:=apple[tov[k]];
                                            if l>0 then
                                            begin
                                                    flag:=false;
                                                    if l=1 then zy:=zy+mi[tov[k]]
                                                    else
                                                    if l=2 then zy:=zy-mi[tov[k]];
                                                    inc(gs);
                                            end;
                                            k:=next[k];
                                    end;
                                    if flag then
                                    begin
                                            p:=j;
                                            for k:=1 to zs do
                                            begin
                                                    if apple[k]>0 then
                                                    begin
                                                            if map[j,k]<map[j,p] then
                                                            begin
                                                                    p:=k;
                                                            end;
                                                    end;
                                            end;
                                            zy:=zy+mi[j];
                                            l:=i mod mi[p+1] div mi[p];
                                            if l=1 then zy:=zy+mi[p]
                                            else
                                            if l=2 then zy:=zy-mi[p];
                                            if map[j,p]<maxx then
                                            f[zy]:=min(f[zy],f[i]+map[j,p]);
                                    end
                                    else
                                    begin
                                            if gs mod 2=0 then gs:=2
                                            else gs:=1;
                                            zy:=zy+gs*mi[j];
                                            f[zy]:=min(f[zy],f[i]);
                                    end;
                            end;
                    end;
            end;
    
    
            for i:=1 to mi2[zs+1]-1 do
            begin
                    j:=0;
                    for k:=1 to zs do
                    begin
                            if i and mi2[k]>0 then
                            begin
                                    j:=k;
                                    break;
                            end;
                    end;
                    dp[i]:=maxx;
                    for k:=1 to zs do
                    begin
                            if (i and mi2[k]>0) and (k<>j) then
                            begin
                                    dp[i]:=min(dp[i],dp[i-mi2[k]-mi2[j]]+map[j,k]);
                            end;
                    end;
            end;
            tot1:=mi[zs+1]-1;
            for i:=1 to zs do
            begin
                    if not bz[i] then
                    begin
                            tot1:=tot1-mi[i]*2;
                    end;
            end;
            for i:=2 to mi[zs+1]-1 do
            begin
                    if f[i]=maxx then continue;
                    for j:=1 to zs do
                    begin
                            apple[j]:=i mod mi[j+1] div mi[j];
                    end;
                    zy:=i;
                    flag1:=true;
                    for j:=1 to zs do
                    begin
                            zy:=i;
                            if ((apple[j])=0) and (bz[j]) then
                            begin
                                    flag1:=false;
                                    break;
                            end;
                    end;
                    if flag1 then
                    begin
                            pp[0]:=0;
                            add:=0;
                            zy:=0;
                            for j:=1 to zs do
                            begin
                                    if (apple[j])=1 then
                                    begin
                                            inc(pp[0]);
                                            pp[pp[0]]:=j;
                                            zy:=zy+mi2[j];
                                    end;
                            end;
                            if pp[0] mod 2=1 then continue;
                            add:=dp[zy];
                            if add<maxx then
                            begin
                                    f[tot1]:=min(f[tot1],f[i]+add);
                            end;
                    end;
            end;
            writeln(f[tot1]+answer);
    end.
    
    我活在这夜里。无论周围多么黑暗,我都要努力发光!我相信着,终有一天,我会在这深邃的夜里,造就一道最美的彩虹。
  • 相关阅读:
    FLV视频转换的利器 ffmpeg.exe
    ffmpeg参数设定解说
    SQL里加减日期
    SQL Server 获得影响行数
    CheckBoxList RadioButtonList 不生成table 表示
    SQL语句 从一个表读取数据,写入到另一个表的相同字段中
    ffmpeg和Mencoder使用实例小全
    执行存储过程出现:"不是有效的标识符。"
    SQL 格式化超长的字段
    js遍历选中的dom元素
  • 原文地址:https://www.cnblogs.com/RainbowCrown/p/11148411.html
Copyright © 2020-2023  润新知