注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Yavin(某沙茶的代码库)

Star Wars fan and OIer

 
 
 

日志

 
 

ZJOI2013 defend simplex algorithm  

2013-04-07 17:23:31|  分类: BZOJ |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
以下才是正确版本,参考了vani,vfleaking和维基百科的代码
常数很渣,不过轻松过无压力。
几个关键点吧——最大化,于是每次找递减的转轴后值变化最小的

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<fstream>
#include<map>
#include<ctime>
#include<set>
#include<queue>
#include<cmath>
#include<vector>
#include<bitset>
#include<functional>
#define x first
#define y second
#define mp make_pair
#define pb push_back
#define REP(i,l,r) for((i)=(l);(i)<=(r);++(i))
#define REP2(i,l,r) for((i)=(l);(i)!=(r);++(i))
using namespace std;

typedef long long LL;
typedef double ld;

const int MAXM=12000+10;
const int MAXN=1200+10;
const LL MAXS=10000;

int n,m;
LL g[MAXM][MAXN];
pair<LL,LL> q[MAXM];
int top;
int cost[MAXN];

void orz(int A,int B)
{
int i,j;
int tmp=-g[A][B];
g[A][B]=-1;
top=0;
REP(i,0,n)
if(g[A][i])
{
g[A][i]/=tmp;
q[++top]=mp(i,g[A][i]);
}
REP(i,0,m)
if(g[i][B] && i!=A)
{
int tmp=g[i][B];
g[i][B]=0;
REP(j,1,top)
g[i][q[j].x]+=tmp*q[j].y;
}
}

LL solve()
{
while(1)
{
int i,me=0,min=0;
REP(i,1,n)
if(g[0][i]>0)
{
me=i;
break;
}
if(!me)
return g[0][0];
REP(i,1,m)
if(g[i][me]<0 && ( min==0 || g[i][0]*g[min][me]>g[min][0]*g[i][me]))
min=i;
orz(min,me);
}
return 0;
}

int main()
{
#ifndef ONLINE_JUDGE
freopen("input.txt","r",stdin);freopen("output.txt","w",stdout);
#endif
int i,j;
scanf("%d%d",&n,&m);
LL ans=0;
REP(i,1,n)
{
scanf("%d",&cost[i]);
g[0][i]=cost[i];
ans+=cost[i]*MAXS;
}
REP(i,1,m)
{
int l,r,a;
scanf("%d%d%d",&l,&r,&a);
REP(j,l,r)
g[i][j]=-1;
g[i][0]=MAXS*(r-l+1)-a;
}
REP(i,1,n)
{
++m;
g[m][i]=-1;
g[m][0]=MAXS;
}
printf("%lld\n",ans-solve());
return 0;
}



  评论这张
 
阅读(363)| 评论(5)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017