【BZOJ】【1391】【CEOI2008】order

网络流/最小割


  暴力建图就好了……S->i 容量为收益,i->j+n 容量为租金,j+n->T容量为购买所花的钱。

  如果亏钱的话那么割掉的就是收益,表示不赚钱。

  如果租金大于购买所花的钱就会割掉购买的钱(因为流量限制住了……)

  表示T到死啊……copy了Hzwer……orz

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4 #define inf 0x7fffffff
 5 using namespace std;
 6 int T,n,m,cnt=1,ans,cur[2501],q[2505],head[2505],h[2505];
 7 struct data{int to,next,v;}e[3000001];
 8 void ins(int u,int v,int w)
 9 {cnt++;e[cnt].to=v;e[cnt].v=w;e[cnt].next=head[u];head[u]=cnt;}
10 void insert(int u,int v,int w)
11 {ins(u,v,w);ins(v,u,0);}
12 bool bfs()
13 {
14      int t=0,w=1,i,now;
15      for(int i=1;i<=T;i++)h[i]=-1;
16      q[0]=h[0]=0;
17      while(t!=w)
18      {
19             now=q[t];t++;if(t==2501)t=0;
20             for(i=head[now];i;i=e[i].next)
21             {
22                   if(e[i].v&&h[e[i].to]<0)
23                         {h[e[i].to]=h[now]+1;q[w++]=e[i].to;if(w==2501)w=0;}
24              }
25      }
26      if(h[T]==-1)return 0;return 1;
27      }
28 int dfs(int x,int f)
29 {
30     if(x==T)return f;
31     int w,used=0;
32     for(int i=cur[x];i;i=e[i].next)
33     {
34             if(e[i].v&&h[e[i].to]==h[x]+1)
35             {
36                 w=f-used;
37                 w=dfs(e[i].to,min(w,e[i].v));   
38                 e[i].v-=w;if(e[i].v>0)cur[x]=i;e[i^1].v+=w;
39                 used+=w;if(used==f)return f;                      
40                 }
41             }
42     if(!used)h[x]=-1;
43     return used;
44     }
45 void dinic(){while(bfs()){for(int i=0;i<=T;i++)cur[i]=head[i];ans-=dfs(0,inf);}}
46 int main()
47 {
48     scanf("%d%d",&n,&m);
49     T=n+m+1;
50     int a,b,c,d;
51     for(int i=1;i<=n;i++)
52     {
53         scanf("%d%d",&a,&b);
54         insert(0,i,a);ans+=a;
55         for(int j=1;j<=b;j++)
56         {
57             scanf("%d%d",&c,&d);
58             insert(i,n+c,d);
59         }
60     }
61     for(int i=1;i<=m;i++)
62     {
63         scanf("%d",&a);
64         insert(n+i,T,a);
65     }
66     dinic();
67     printf("%d",ans);
68     return 0;
69 }
View Code(Hzwer)
 1 //BZOJ 2768
 2 #include<vector>
 3 #include<cstdio>
 4 #include<cstring>
 5 #include<cstdlib>
 6 #include<iostream>
 7 #include<algorithm>
 8 #define rep(i,n) for(int i=0;i<n;++i)
 9 #define F(i,j,n) for(int i=j;i<=n;++i)
10 #define D(i,j,n) for(int i=j;i>=n;--i)
11 #define pb push_back
12 using namespace std;
13 inline int getint(){
14     int v=0,sign=1; char ch=getchar();
15     while(ch<'0'||ch>'9'){ if (ch=='-') sign=-1; ch=getchar();}
16     while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();}
17     return v*sign;
18 }
19 const int N=2510,M=3000000,INF=~0u>>2;
20 typedef long long LL;
21 /******************tamplate*********************/
22 int n,m,ans;
23 struct edge{
24     int to,v;
25 };
26 edge E[M];
27 int head[N],next[M],cnt;
28 void add(int x,int y,int v){
29     E[++cnt]=(edge){y,v};
30     next[cnt]=head[x]; head[x]=cnt;
31     E[++cnt]=(edge){x,0};
32     next[cnt]=head[y]; head[y]=cnt;
33 }
34 int s,t,cur[N],d[N],Q[N];
35 void init(){
36     n=getint();m=getint();
37     cnt=1;ans=0;
38     s=0; t=n+m+1;
39     int x,y,a,b;
40     F(i,1,n){
41         x=getint(); y=getint();
42         ans+=x;
43         add(s,i,x);
44         F(j,1,y){
45             a=getint(); b=getint();
46             add(i,n+a,b);
47         }
48     }
49     F(i,1,m){
50         x=getint();
51         add(n+i,t,x);
52     }
53 }
54 bool mklevel(){
55     F(i,s,t) d[i]=-1;
56     d[s]=0;
57     int l=0,r=-1;
58     Q[++r]=s;
59     while(l<=r){
60         int x=Q[l++];
61         for(int i=head[x];i;i=next[i])
62             if (d[E[i].to]==-1 && E[i].v){
63                 d[E[i].to]=d[x]+1;
64                 Q[++r]=E[i].to;
65             }
66     }
67     return d[t]!=-1;
68 }
69 int dfs(int x,int a){
70     if (x==t||a==0) return a;
71     int flow=0;
72     for(int &i=cur[x];i && flow<a;i=next[i])
73         if (d[E[i].to]==d[x]+1 && E[i].v){
74             int f=dfs(E[i].to,min(a-flow,E[i].v));
75             E[i].v-=f;
76             E[i^1].v+=f;
77             flow+=f;
78         }
79     if (!flow) d[x]=-1;
80     return flow;
81 }
82 void Dinic(){
83     while(mklevel()){
84         F(i,s,t) cur[i]=head[i];
85         ans-=dfs(s,INF);
86     }
87 }
88 int main(){
89 #ifndef ONLINE_JUDGE
90     freopen("2768.in","r",stdin);
91     freopen("2768.out","w",stdout);
92 #endif
93     init();Dinic();
94     printf("%d\n",ans);
95     return 0;
96 }
View Code(My)

 

1391: [Ceoi2008]order

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 746  Solved: 230
[Submit][Status][Discuss]

Description

有N个工作,M种机器,每种机器你可以租或者买过来. 每个工作包括若干道工序,每道工序需要某种机器来完成,你可以通过购买或租用机器来完成。 现在给出这些参数,求最大利润

Input

第一行给出 N,M(1<=N<=1200,1<=M<=1200) 下面将有N块数据,每块数据第一行给出完成这个任务能赚到的钱(其在[1,5000])及有多少道工序 接下来若干行每行两个数,分别描述完成工序所需要的机器编号及租用它的费用(其在[1,20000]) 最后M行,每行给出购买机器的费用(其在[1,20000])

Output

最大利润

Sample Input

2 3
100 2
1 30
2 20
100 2
1 40
3 80
50
80
110

Sample Output

50

HINT

Source

[ Submit][ Status][ Discuss]

转载于:https://www.cnblogs.com/Tunix/p/4338469.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 有一个 $n$ 个点的棋盘,每个点上有一个数字 $a_i$,你需要从 $(1,1)$ 走到 $(n,n)$,每次只能往右或往下走,每个格子只能经过一次,路径上的数字和为 $S$。定义一个点 $(x,y)$ 的权值为 $a_x+a_y$,求所有满足条件的路径中,所有点的权值和的最小值。 输入格式 第一行一个整数 $n$。 接下来 $n$ 行,每行 $n$ 个整数,表示棋盘上每个点的数字。 输出格式 输出一个整数,表示所有满足条件的路径中,所有点的权值和的最小值。 数据范围 $1\leq n\leq 300$ 输入样例 3 1 2 3 4 5 6 7 8 9 输出样例 25 算法1 (树形dp) $O(n^3)$ 我们可以先将所有点的权值求出来,然后将其看作是一个有权值的图,问题就转化为了在这个图中求从 $(1,1)$ 到 $(n,n)$ 的所有路径中,所有点的权值和的最小值。 我们可以使用树形dp来解决这个问题,具体来说,我们可以将这个图看作是一棵树,每个点的父节点是它的前驱或者后继,然后我们从根节点开始,依次向下遍历,对于每个节点,我们可以考虑它的两个儿子,如果它的两个儿子都被遍历过了,那么我们就可以计算出从它的左儿子到它的右儿子的路径中,所有点的权值和的最小值,然后再将这个值加上当前节点的权值,就可以得到从根节点到当前节点的路径中,所有点的权值和的最小值。 时间复杂度 树形dp的时间复杂度是 $O(n^3)$。 C++ 代码 算法2 (动态规划) $O(n^3)$ 我们可以使用动态规划来解决这个问题,具体来说,我们可以定义 $f(i,j,s)$ 表示从 $(1,1)$ 到 $(i,j)$ 的所有路径中,所有点的权值和为 $s$ 的最小值,那么我们就可以得到如下的状态转移方程: $$ f(i,j,s)=\min\{f(i-1,j,s-a_{i,j}),f(i,j-1,s-a_{i,j})\} $$ 其中 $a_{i,j}$ 表示点 $(i,j)$ 的权值。 时间复杂度 动态规划的时间复杂度是 $O(n^3)$。 C++ 代码

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值