题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3118

题意:给出一个图以及图中指定的n-1条边组成的生成树。每条边权值加1或者减去1都有相应的代价。求一个最小代价使得给出的边是最小生成树。

思路:对于每条非树边,必与某些树边形成环。设该非树边的权值为w2,某树边的权值为 w1。最后非树边增加x2,树边减少x1,那么w1-x1<=w2+x2。这样我们可以得到一些式子。代价也知道,这样就转化成线性规划问题。题目求的是最小值,我们可以将目标方程的系数取反求最大值。

单纯形的步骤:

(1)求出一个初始解;

(2)迭代。

(1)这个题的系数矩阵A是全么模:1、元素都是0,-1,1;2、任意子方阵的行列式为0,-1,1。

(2)据说A是全么模时解是整数解,因此此题可直接用单纯形。

const int COL=1005;
const int ROW=30005;

int n,m,B[ROW],N[COL];
double A[ROW][COL],b[ROW],c[COL],v;
double ans[COL];

int sgn(double x)
{
    if(x>1e-8) return 1;
    if(x<-1e-8) return -1;
    return 0;
}

//B中第l个替换N中第e个

void pivot(int l,int e)
{
    int i,j;
    double temp=A[l][e];
    b[l]/=temp; A[l][e]=1/temp;
    for(i=1;i<=n;i++) if(i!=e) A[l][i]/=temp;
    for(i=1;i<=m;i++) if(i!=l)
    {
        b[i]-=A[i][e]*b[l];
        for(j=1;j<=n;j++) if(j!=e) A[i][j]-=A[i][e]*A[l][j];
        A[i][e]=-A[i][e]/temp;
    }
    v+=b[l]*c[e];
    for(i=1;i<=n;i++) if(i!=e) c[i]-=c[e]*A[l][i];
    c[e]*=-A[l][e];
    swap(B[l],N[e]);
}

void simplex()
{
    int i,j,k,x;
    int l,s;
    double temp,temp1,temp2,temp3;

    while(1)
    {
        temp2=-dinf; s=-1;
        for(i=1;i<=n;i++) if(sgn(c[i])>0)
        {
            temp=dinf;
            for(k=1;k<=m;k++) if(sgn(A[k][i])>0)
            {
                temp3=b[k]/A[k][i];
                if(temp3<temp) temp=temp3,x=k;
            }
            if(temp2<temp*c[i])
            {
                s=i,l=x,temp2=temp*c[i];
            }
        }
        if(s==-1) break;
        pivot(l,s);
    }
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++) if(B[j]==i) break;
        if(j<=m) ans[i]=b[j];
        else ans[i]=0;
    }
}

void print()
{
    int i,j;
    printf("v: %.3lf\n",v);
    printf("c:\n");
    for(i=1;i<=n;i++) printf("%.3lf ",c[i]); puts("");
    printf("A:\n");
    for(i=1;i<=m;i++)
    {
        for(j=1;j<=n;j++) printf("%.3lf ",A[i][j]);
        puts("");
    }
    printf("b:\n");
    for(i=1;i<=m;i++) printf("%.3lf ",b[i]); puts("");
    printf("B:\n");
    for(i=1;i<=m;i++) printf("%d ",B[i]); puts("");

    printf("N:\n");
    for(i=1;i<=n;i++) printf("%d ",N[i]); puts("");

}

int init()
{
    int i,j;
    int k=1;
    for(i=1;i<=m;i++) if(b[i]<b[k]) k=i;
    if(sgn(b[k])>=0)
    {
        for(i=1;i<=n;i++) N[i]=i;
        for(i=1;i<=m;i++) B[i]=n+i;
        v=0;
        simplex();
        return 1;
    }

    static double tmpC[COL];
    for(i=1;i<=n;i++) tmpC[i]=c[i];
    tmpC[n+1]=0;

    n++;
    for(i=1;i<=m;i++) A[i][n]=-1;
    for(i=1;i<=n;i++) N[i]=i;
    for(i=1;i<=m;i++) B[i]=n+i;
    v=0;
    for(i=1;i<=n;i++)
    {
        if(i<n) c[i]=0;
        else c[i]=-1;
    }
    pivot(k,n);
    simplex();

    if(sgn(ans[n])!=0) return 0;

    static int belongB[COL];
    clr(belongB,0);
    for(i=1;i<=m;i++)
    {
        if(B[i]>n) continue;
        belongB[B[i]]=i;
    }
    map<int,int> mp;
    for(i=1;i<=n;i++) mp[N[i]]=i;
    clr(c,0);
    v=0;
    for(i=1;i<=n;i++)
    {
        if(!belongB[i])
        {
            c[mp[i]]+=tmpC[i];
        }
        else
        {
            v+=tmpC[i]*b[belongB[i]];
            int j;
            for(j=1;j<=n;j++)
            {
                c[j]+=tmpC[i]*(-A[belongB[i]][j]);
            }
        }
    }
    c[mp[n]]=0;
    for(i=1;i<=m;i++) A[i][mp[n]]=0;

    simplex();
    n--;
    return 1;
}

struct node
{
    int u,v,id,w,next;
};

node edges[ROW];
int e;
int head[COL];

void add(int u,int v,int w,int id)
{
    e++;
    edges[e].u=u;
    edges[e].v=v;
    edges[e].w=w;
    edges[e].id=id;
    edges[e].next=head[u];
    head[u]=e;
}

int h[COL],up[COL],down[COL];
int eNum;
int inq[COL],KK;
int pre[COL];

void build(int s,int t,int p)
{
    int i;
    for(i=pre[t];i!=-1;)
    {
        int w1=edges[i].w;
        int w2=edges[p*2].w;

        m++;
        A[m][edges[i].id]=-1;
        A[m][p]=-1;
        b[m]=w2-w1;

        int u=edges[i].u;
        i=pre[u];
    }
}

void bfs(int s,int t,int p)
{
    queue<int> Q;
    Q.push(s);
    KK++;
    inq[s]=KK;
    pre[s]=-1;
    while(!Q.empty())
    {
        int u=Q.front();
        Q.pop();

        int i;
        for(i=head[u];i!=-1;i=edges[i].next)
        {
            int v=edges[i].v;
            int id=edges[i].id;
            if(!h[id]||KK==inq[v]) continue;
            pre[v]=i;
            if(v==t)
            {
                build(s,t,p);
                return;
            }
            Q.push(v);
            inq[v]=KK;
        }
    }
}

int main()
{
    n=myInt();
    eNum=myInt();
    clr(head,-1);
    int i;
    for(i=1;i<=eNum;i++)
    {
        int u,v,w;
        scanf("%d%d%d%d%d%d",&u,&v,&w,&h[i],&up[i],&down[i]);
        add(u,v,w,i);
        add(v,u,w,i);
    }
    for(i=1;i<=eNum;i++) if(!h[i])
    {
        int t=i*2;
        bfs(edges[t].u,edges[t].v,i);
    }
    for(i=1;i<=eNum;i++) c[i]=h[i]?-down[i]:-up[i];

    n=eNum;

    if(!init())
    {
        puts("no solution");
    }

    double res=-v;
    if(res<1e-10) res=0;
    printf("%.0lf\n",-v);
}

bzoj 3118: Orz the MST(单纯形)的更多相关文章

  1. BZOJ 3118 Orz the MST

    权限题qwq 如果我们要使得某棵生成树为最小生成树,那么上面的边都不能被替代,具体的,对于一个非树边,它的权值要\(\ge\)它两端点在树上的路径上的所以边的权值,所以对于每个非树边就可以对一些树边列 ...

  2. BZOJ3118 Orz the MST 【单纯形 + 生成树】

    题目链接 BZOJ3118 题解 少有的单纯形好题啊 我们先抽离出生成树 生成树中的边只可能减,其它边只可能加 对于不在生成树的边,其权值一定要比生成树中其端点之间的路径上所有的边都大 然后就是一个最 ...

  3. BZOJ 2654 & 玄学二分+MST

    题意: 给一张图,边带权且带颜色黑白,求出一棵至少包含k条白边的MST SOL: 正常人都想优先加黑边或者是白边,我也是这么想的...你看先用白边搞一棵k条边的MST...然后维护比较黑边跟白边像堆一 ...

  4. BZOJ 3112 Zjoi2013 防守战线 单纯形

    题目大意: 单纯形*2.. . #include <cmath> #include <cstdio> #include <cstring> #include < ...

  5. BZOJ 3265 志愿者招募加强版(单纯形)

    3265: 志愿者招募加强版 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 848  Solved: 436[Submit][Status][Disc ...

  6. BZOJ3118 : Orz the MST

    对于树边显然只需要减少权值,对于非树边显然只需要增加权值 设i不为树边,j为树边 X[i]:i增加量 X[j]:j减少量 C[i]:修改1单位i的代价 对于每条非树边i(u,v),在树上u到v路径上的 ...

  7. BZOJ 2654: tree( 二分 + MST )

    我们给白色的边增加权值 , 则选到的白色边就会变多 , 因此可以二分一下. 不过这道题有点小坑... ------------------------------------------------- ...

  8. bzoj [Noi2008] 1061 志愿者招募 单纯形

    [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 5437  Solved: 3267[Submit][Status][Di ...

  9. bzoj3118: Orz the MST(线性规划+单纯形法)

    传送门 不难发现,对于每一条树边肯定要减小它的权值,对于每一条非树边要增加它的权值 对于每一条非树边\(j\),他肯定与某些树边构成了一个环,那么它的边权必须大于等于这个环上的所有边 设其中一条边为\ ...

随机推荐

  1. 夺命雷公狗—angularjs—9—ng-class的自定义函数的用法

    angularjs里面其实给我们留下了一个很不错的地方,他就是可以直接调用函数从而对该位置进行处理, 被点击后展示效果如下所示: 开始走代码吧.... <!doctype html> &l ...

  2. 夺命雷公狗---node.js---10之POST的接收

    首先我们在项目下创建一个表单,代码如下所示: <!DOCTYPE html> <html lang="en"> <head> <meta ...

  3. ThinkPHP讲解(五)——数据库配置及Model数据模型层、查询

    数据库配置 在TP框架中要进行连接数据库操作,要进行配置 要在convertion.php中找到“数据库配置”,并复制到项目配置文件config.php中 Model模型层制作 model:数据库中每 ...

  4. [ubuntu] Can not run OpenProj on Ubuntu

    Download openproj from http://sourceforge.net/projects/openproj/ Yes, you can only download the rpm ...

  5. Sql server analysis service 通过IIS连接时的最大连接数问题

    做过SSAS项目的大部分人都应该知道SSAS是可以通过在IIS上建立代理站点的方式来建立远程连接的,这样可以绕过连接到SSAS时需要在同一个域环境下通过域用户来验证的问题,这样即使连接到SSAS的客户 ...

  6. 用javascript在客户端删除某一个cookie键值对

    下面这个方法展示如何在客户端浏览器上用javascript删除某一个cookie键值对. //用javascript删除某一个cookie的方法,该方法传入要删除cookie的名称 function ...

  7. android自定义view属性

    第一种 /MainActivity/res/values/attrs.xml <?xml version="1.0" encoding="utf-8"?& ...

  8. 鸟哥的linux私房菜之档案与文件系统的压缩与打包

    00000001 节约空间 其实简单的说压缩就是把没有用到的0给去掉,解压的时候在加上 在linux中,压缩文件档案的扩展名大多是.tar,.tar.gz,tgz,gz,.Z,.bz2 compres ...

  9. java中length,length(),size()的区别

    1. java中的length属性是针对数组说的,比如说你声明了一个数组,想知道这个数组的长度则用到了length这个属性.2. java中的length()方法是针对字符串String说的,如果想看 ...

  10. Java 基本数据类型 sizeof 功能【转】

    转自:http://blog.csdn.net/sunboy_2050/article/details/7310008 版权声明:本文为博主原创文章,未经博主允许不得转载. Java基本数据类型int ...