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.
我活在这夜里。无论周围多么黑暗,我都要努力发光!我相信着,终有一天,我会在这深邃的夜里,造就一道最美的彩虹。
原文地址:https://www.cnblogs.com/RainbowCrown/p/11148411.html