随笔 - 531  文章 - 0  评论 - 3  阅读 - 10215 

神犇航空有K架飞机,为了简化问题,我们认为每架飞机都是相同的。神犇航空的世界中有N个机场,以0⋯N-1编号,

其中0号为基地机场,每天0时刻起飞机才可以从该机场起飞,并不晚于TT时刻回到该机场。

一天,神犇航空接到了MM个包机请求,每个请求为在ss时刻从aa机场起飞,在恰好tt时刻到达bb机场,可以净获利cc。

设计一种方案,使得总收益最大。

 

 

 任务当作点,点权为1,然后拆点,连边  (1, Val [i]  )

 S-> begin( 1, 0) ,end-> T (1, 0) 

 

任务之间放边(  1,   - cost[ i ][ j ] ) 

最后再建一个 S‘' --> S ,( K, 0)

 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
#include <iostream>
#include <queue>
#include <cstring>
using namespace std;
const int N =500, M=1e5+100;
 
 #define inf 1e9
 int all=1,hd[N],go[M],w[M],nxt[M],cost[M];
  
 int a[N][N],tm[N][N];
 int qa[N],qb[N],qs[N],qe[N],qc[N];
 int TIME,K;
 int S,T,n,m;
 int pre[N],mn[N],dis[N],vis[N],ans=0;
  
 void add_(int x,int y,int z,int c){
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=z; cost[all]=c;
      
     swap(x,y);
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=0; cost[all]= -c;
 }
  
  bool spfa__max(){
    int i;
    queue<int> q;
    q.push(S);
    for(i=0;i<N;i++) dis[i]= -inf;
    for(i=0;i<N;i++) vis[i]=0;
    mn[S]=inf, dis[S]=0, vis[S]=1;
      
    while(!q.empty()){
        int x=q.front();
        vis[x]=0;
        q.pop();
        for(i=hd[x];i;i=nxt[i]){
            int y=go[i],z=w[i],c=cost[i];
            if(w[i]==0) continue;
              
            if(dis[x]+c>dis[y]){
                dis[y]=dis[x]+c;
                mn[y]=min(mn[x],z);
                pre[y]=i;
                if(vis[y]==0) vis[y]=1,q.push(y);
            }
        }
    }
    return dis[T]== -inf ?0:1;
 }
 
 void solve(){
    int maxf=0, minc=0;
    while(spfa__max()){
        int x= T;
        maxf+=  mn[T] ,minc+= mn[T]*dis[T];
         
        while(x!=S){
            int i=pre[x];
            w[i]-= mn[T];
            w[i^1]+=mn[T];
            x=go[i^1];
        }
    }
    cout<< minc<<endl;
 }
  
 signed main(){
    int i,j,x;
    cin>>n>>m>>K>>TIME;
    S=m*2+1,T=S+2;
     
    for(i=1;i<=n;++i)
     for(j=1;j<=n;++j){
        cin>>tm[i][j];
     }
     for(i=1;i<=n;++i)
     for(j=1;j<=n;++j)
        cin>>a[i][j];
      
     add_(S,S+1,K,0);
     for(int i=1;i<=m;i++)
     cin>>qa[i]>>qb[i]>>qs[i]>>qe[i]>>qc[i],
     qa[i]++,qb[i]++;
      
     for(int i=1;i<=m;i++){
        add_(i,i+m,1, qc[i]);
         
        if(tm[1][qa[i]]<=qs[i])
             add_(S+1,i,1, -a[1][qa[i]]);
        if(qe[i]+tm[qb[i]][1]<=TIME)
             add_(i+m,T,1, -a[qb[i]][1]);
              
        for(int j=1;j<=m;j++){
            if(i!=j && qe[i]+tm[qb[i]][qa[j]]<=qs[j])
            add_(i+m,j,1,-a[qb[i]][qa[j]]);
        }
    }
     solve();
  }

 

 

 

 

    

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
#include <iostream>
#include <queue>
#include <cstring>
using namespace std;
const int N =500, M=1e5+100;
 
 #define inf 1e9
 int all=1,hd[N],go[M],w[M],nxt[M],cost[M];
  
 int a[N][N],tm[N][N];
 int qa[N],qb[N],qs[N],qe[N],qc[N];
 int TIME,K;
 int S,T,n,m;
 int pre[N],mn[N],dis[N],vis[N],ans=0;
  
 void add_(int x,int y,int z,int c){
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=z; cost[all]=c;
      
     swap(x,y);
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=0; cost[all]= -c;
 }
  
  bool spfa__max(){
    int i;
    queue<int> q;
    q.push(S);
    for(i=0;i<N;i++) dis[i]= -inf;
    for(i=0;i<N;i++) vis[i]=0;
    mn[S]=inf, dis[S]=0, vis[S]=1;
      
    while(!q.empty()){
        int x=q.front();
        vis[x]=0;
        q.pop();
        for(i=hd[x];i;i=nxt[i]){
            int y=go[i],z=w[i],c=cost[i];
            if(w[i]==0) continue;
              
            if(dis[x]+c>dis[y]){
                dis[y]=dis[x]+c;
                mn[y]=min(mn[x],z);
                pre[y]=i;
                if(vis[y]==0) vis[y]=1,q.push(y);
            }
        }
    }
    return dis[T]== -inf ?0:1;
 }
 
 void solve(){
    int maxf=0, minc=0;
    while(spfa__max()){
        int x= T;
        maxf+=  mn[T] ,minc+= mn[T]*dis[T];
         
        while(x!=S){
            int i=pre[x];
            w[i]-= mn[T];
            w[i^1]+=mn[T];
            x=go[i^1];
        }
    }
    cout<< minc<<endl;
 }
  
 signed main(){
    int i,j,x;
    cin>>n>>m>>K>>TIME;
    S=m*2+1,T=S+2;
     
    for(i=1;i<=n;++i)
     for(j=1;j<=n;++j){
        cin>>tm[i][j];
     }
     for(i=1;i<=n;++i)
     for(j=1;j<=n;++j){
        cin>>a[i][j];
     }
     add_(S,S+1,K,0);
     for(int i=1;i<=m;i++)
     cin>>qa[i]>>qb[i]>>qs[i]>>qe[i]>>qc[i],
     qa[i]++,qb[i]++;
      
     for(int i=1;i<=m;i++){
        add_(i,i+m,1, qc[i]);
         
        if(tm[1][qa[i]]<=qs[i])
             add_(S+1,i,1, - a[1][qa[i]]);
        if(qe[i]+tm[qb[i]][1]<=TIME)
             add_(i+m,T,1,  -a[qb[i]][1]);
              
        for(int j=1;j<=m;j++){
            if(i!=j && qe[i]+tm[qb[i]][qa[j]]<=qs[j])
            add_(i+m,j,1,-a[qb[i]][qa[j]]);
        }
    }
     solve();
  }

 

posted on   towboat  阅读(35)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!
点击右上角即可分享
微信分享提示