bzoj1061、3265 志愿者招募【单纯形法】

解题思路：

如果不清楚单纯形法的可以看这里：
http://blog.sina.com.cn/s/blog_4b1046f80101m9ou.html
题目很像标准单纯形形式但求的是最大值，且所有限制条件都是大于等于。
根据对偶原理，我们只需要把单纯形法整个ABC矩阵转置再求单纯形即是答案。

这是1061的代码，3265的改一下输入即可。

#include<bits/stdc++.h>
using namespace std;  

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=1005,M=10005;
const double eps=1e-8;
int n,m;
double ans,a[M][N],b[M],c[N];

void pivot(int l,int e)
{
    b[l]/=a[l][e];
    for(int i=1;i<=n;i++)
        if(i!=e)a[l][i]/=a[l][e];
    a[l][e]=1.0/a[l][e];

    for(int i=1;i<=m;i++)
        if(i!=l&&abs(a[i][e])>eps)
        {
            b[i]-=a[i][e]*b[l];
            for(int j=1;j<=n;j++)
                if(j!=e)a[i][j]-=a[i][e]*a[l][j];
            a[i][e]=-a[i][e]*a[l][e];   
        }

    ans+=c[e]*b[l];
    for(int i=1;i<=n;i++)
        if(i!=e)c[i]-=c[e]*a[l][i];
    c[e]=-c[e]*a[l][e];
}

void simplex()
{
    int i,e,l;
    while(true)
    {
        for(i=1;i<=n;i++)
            if(c[i]>eps){e=i;break;}
        if(i==n+1)return;
        double t=1e20;
        for(i=1;i<=m;i++)
            if(a[i][e]>eps&&b[i]/a[i][e]<t)
                l=i,t=b[i]/a[i][e];
        if(t==1e20){ans=t;return;}
        pivot(l,e);
    }
}

int main()
{
    //freopen("lx.in","r",stdin);
    n=getint(),m=getint();
    for(int i=1;i<=n;i++)c[i]=getint();
    for(int i=1;i<=m;i++)
    {
        int l=getint(),r=getint();b[i]=getint();
        for(int j=l;j<=r;j++)a[i][j]++;
    }
    simplex();
    printf("%d",(int)(ans+0.5));
    return 0;
}