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)); }
Copyright:http://www.cnblogs.com/candy99/
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
· 程序员转型AI:行业分析
· 重磅发布!DeepSeek 微调秘籍揭秘,一键解锁升级版全家桶,AI 玩家必备神器!