BZOJ 3265: 志愿者招募加强版 [单纯形法]

传送门


 

一个人多段区间,一样....

不过国家队论文上说这道题好像不能保证整数解....

复制代码
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=1005,M=1e4+5;
const double INF=1e15,eps=1e-8;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

int n,m;
double a[M][N];
int q[N];
void Pivot(int l,int e){
    double t=a[l][e];a[l][e]=1;
    for(int j=0;j<=n;j++) a[l][j]/=t;
    int p=0;
    for(int j=0;j<=n;j++) if(abs(a[l][j])>eps) q[++p]=j;
    for(int i=0;i<=m;i++) if(i!=l && abs(a[i][e])>eps){
        double t=a[i][e];a[i][e]=0;
        for(int j=1;j<=p;j++) a[i][q[j]]-=t*a[l][q[j]];
    }
}
void simplex(){
    while(true){
        int l=0,e=0; double mn=INF;
        for(int j=1;j<=n;j++) if(a[0][j]>eps) {e=j;break;}
        if(!e) return;
        for(int i=1;i<=m;i++) 
            if(a[i][e]>eps && a[i][0]/a[i][e]<mn) {mn=a[i][0]/a[i][e];l=i;}
        if(!l) return;//unbounded
        Pivot(l,e);
    }
}
int main(){
    freopen("in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;i++)a[0][i]=read();
    for(int i=1;i<=m;i++){
        int k=read();
        while(k--){
            int l=read(),r=read();
            for(int j=l;j<=r;j++) a[i][j]=1;
        }
        a[i][0]=read();
    }
    simplex();
    printf("%d",int(-a[0][0]+0.5));
}
复制代码

 

posted @   Candy?  阅读(334)  评论(0编辑  收藏  举报
编辑推荐:
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
阅读排行:
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
· 程序员转型AI:行业分析
· 重磅发布!DeepSeek 微调秘籍揭秘,一键解锁升级版全家桶,AI 玩家必备神器!
点击右上角即可分享
微信分享提示