题目链接: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. BZOJ3118 Orz the MST 【单纯形 + 生成树】

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

  2. BZOJ 2654 &amp; 玄学二分+MST

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

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

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

  4. BZOJ3118 : Orz the MST

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

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

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

  6. 【BZOJ 2753】 2753: [SCOI2012]滑雪与时间胶囊 (分层最小树形图,MST)

    2753: [SCOI2012]滑雪与时间胶囊 Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 2457  Solved: 859 Descriptio ...

  7. BZOJ 1083: [SCOI2005]繁忙的都市(MST)

    裸的最小生成树..直接跑就行了 ---------------------------------------------------------------------- #include<c ...

  8. BZOJ 3479: [Usaco2014 Mar]Watering the Fields( MST )

    MST...一开始没注意-1结果就WA了... ---------------------------------------------------------------------------- ...

  9. BZOJ 1601: [Usaco2008 Oct]灌水( MST )

    MST , kruskal 直接跑 ---------------------------------------------------------------------- #include< ...

随机推荐

  1. word2013删除下载的模板

    word2013删除下载的模板 删除步骤: 1)关闭相关的word文档. 2)按照以下的路径找到相应的位置:"%系统根目录%\Users\Administrator\AppData\Roam ...

  2. ZooKeeper 配置文件(zoo.cfg)详解

    参数名 说明 clientPort 客户端连接server的端口,即对外服务端口,一般设置为2181吧. dataDir 存储快照文件snapshot的目录.默认情况下,事务日志也会存储在这里.建议同 ...

  3. bzoj2592: [Usaco2012 Feb]Symmetry

    Description After taking a modern art class, Farmer John has become interested in finding geometric ...

  4. MVC4 图片上传

    新增 new { enctype = "multipart/form-data" } 这个必须要有 @using (Html.BeginForm(Html.BeginForm(&q ...

  5. CI Weekly #19 | 关于软件开发模型的思考,以及最新 CI/CD 实践分享

    五月一来,夏天便悄然而至.flow.ci 也带来了几个新的变化,帮你进一步优化开发工作流.一起来看看这几个重点功能: 支持 iOS 项目 Xcode8.3 构建 iOSer 们重点来了,flow.ci ...

  6. runtime实现对象存储型数据库——LHDB

    前言 最近在GitHub上看了一份关于基于runtime封装的对象存储型数据库的开源代码,觉得非常值得分享记录一下,在IOS中对数据库的操作一般通过CoreData和SQLite,CoreData 虽 ...

  7. python视频教程全集

    由于Python语言的简洁.易读以及可扩展性,在国外用Python做科学计算的研究机构日益增多,一些知名大学已经用它来教授程序设计课程.如果你想找python视频教程全集,来这就对了!这是链接: py ...

  8. kernel解析dtb为节点

    title: 解析dtb为节点 date: 2019/4/26 14:02:18 toc: true --- kernel解析dtb为节点 head.s入口传递 回顾 看以前的笔记 kernel(二) ...

  9. Ribbon使用Hystrix

    1.导入依赖spring-cloud-starter-hystrix <dependency> <groupId>org.springframework.cloud</g ...

  10. ASP.NET 工作流:支持长时间运行操作的 Web 应用程序

    ASP.NET 工作流 支持长时间运行操作的 Web 应用程序 Michael Kennedy   代码下载位置:MSDN 代码库 在线浏览代码 本文将介绍以下内容: 独立于进程的工作流 同步和异步活 ...