BZOJ 3265 志愿者招募增强版 单
2015-12-09 12:29
405 查看
标题效果:同1061 只是间隔为每种类型的志愿工作是多级
这是卡网络流量?未知
所有在所有的1061您将可以更改为在稍微改变- -
这是卡网络流量?未知
所有在所有的1061您将可以更改为在稍微改变- -
#include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define INF 1e10 #define EPS 1e-7 using namespace std; int m,n; namespace Linear_Programming{ double A[10100][1010],b[10100],c[1010],v; void Pivot(int l,int e) { int i,j; b[l]/=A[l][e]; for(i=1;i<=n;i++) if(i!=e) A[l][i]/=A[l][e]; A[l][e]=1/A[l][e]; for(i=1;i<=m;i++) if(i!=l&&fabs(A[i][e])>EPS) { 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]*A[l][e]; } v+=c[e]*b[l]; for(i=1;i<=n;i++) if(i!=e) c[i]-=c[e]*A[l][i]; c[e]=-c[e]*A[l][e]; } double Simplex() { int i,l,e; while(1) { for(i=1;i<=n;i++) if(c[i]>EPS) break; if((e=i)==n+1) return v; double temp=INF; for(i=1;i<=m;i++) if( A[i][e]>EPS && b[i]/A[i][e]<temp ) temp=b[i]/A[i][e],l=i; if(temp==INF) return INF; Pivot(l,e); } } } int main() { using namespace Linear_Programming; int i,j,k,x,y; cin>>n>>m; for(i=1;i<=n;i++) scanf("%lf",&c[i]); for(i=1;i<=m;i++) { scanf("%d",&k); for(j=1;j<=k;j++) { scanf("%d%d",&x,&y); for(;x<=y;x++) A[i][x]+=1.0; } scanf("%lf",&b[i]); } double ans=Simplex(); printf("%d\n",int(ans+0.5)); return 0; }
相关文章推荐
- 打包android img文件
- Android自定义Button的“款式”
- Android自定义Button的“款式”
- JS实现字符串转日期并比较大小实例分析
- 天声人語 20151209 漱石没後99年
- fedora21 环境下编译 Google V8
- 激活Navicat?如何注册Navicat?
- 最详尽OC 函数库方法翻译
- Android网络编程之UDP通信----简单客户端与服务端通信
- 不该被遗忘的nodeName、nodeValue和nodeType!
- DOM中元素节点、属性节点、文本节点的理解
- 基于QTP的自动化测试留下一份向数据库自动导入测试数据的VBS脚本
- 关于TI CC2540 BLE协议栈中数据传输的几个函数实现(草稿)
- oracle db_link的使用说明
- [Machine Learning] 转载:Kd Tree算法原理和开源实现代码
- springMVC中使用List和Map绑定对象作为参数的使用方法
- CentOS ls 颜色
- PHP php面向对象的入门
- pAdPt_-6 LearnBdMapBasic[(百度)地图(Android SDK)基础]
- javascript中parentNode,childNodes,children的应用详解