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

Yavin(某沙茶的代码库)

Star Wars fan and OIer

 
 
 

日志

 
 

USACO 4.2 Ditch  

2013-04-07 15:36:06|  分类: 默认分类 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
嗯,回到了一年前刷usaco的时光
感触颇多啊,这题我至少提交了60+次。
每次都有新的感受。。。代码如下,可以看看注释

/*
TASK:ditch
LANG:C++
*/
#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 MAXN=600+10;
const int MAXM=200+10;

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

void zhuanzhou(int A,int B)
{
int i,j;
top=0;
//xA=g[A][0]+sigma( g[A][i]*x[i] )
//-g[A][B]*x[B]=g[A][0] ...-1*xA...
//x[B]=g[A][0]+...+
int tmp=-g[A][B];
g[A][B]=-1;
REP(i,0,n)
if(g[A][i])
q[++top]=mp(i,g[A][i]/tmp);
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;
}
}

int danchunxing()
{
//g[0][i]记录目标函数
//g[i][0]+g[i][1]*x[1]+...>=0
int i;
while(1)
{
int me=0;
REP(i,1,n)
if(g[0][i]>0)
{
me=i;
break;
}
if(!me)
return g[0][0];
int min=0;
//选择g[k][0]/g[k][me] 最大的一个
REP(i,1,m)
if(min==0 || ( g[i][me]<0 && g[i][0]*g[min][me]>=g[min][0]*g[i][me] ))
min=i;
zhuanzhou(min,me);
}
return g[0][0];
}

int main()
{
freopen("ditch.in","r",stdin);
freopen("ditch.out","w",stdout);
int i,j,a,b,c,e,v;
scanf("%d%d",&e,&v);
REP(i,1,e)
{
scanf("%d%d%d",&a,&b,&c);
//出入流量相等
g[a-1][i]=1;//Ci,j=g[a][i]*x[i]
g[b-1][i]=-1;//同上
//xi<=c
g[(v-2)*2+i][0]=c;//常数项
g[(v-2)*2+i][i]=-1;
}
REP(i,1,v-1-1)
REP(j,1,e)
g[i+v-2][j]=-g[i][j];
n=e;
m=2*(v-2) + e;
printf("%d\n",danchunxing());
return 0;
}





  评论这张
 
阅读(271)| 评论(1)
推荐 转载

历史上的今天

评论

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

页脚

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