列表

详情


NC20233. [JSOI2016]飞机调度

描述

作为一个旅行达人以及航空公司的金卡会员,JYY 每一年的飞行里程可以 绕赤道几周了。JYY 发现,航空公司为了提高飞机的使用率,并不是简单的一条 航线使用一架飞机来回飞,而是会让同一架飞机连续不停的飞不同的航线,甚至 有的时候为了能够完成飞机的调度,航空公司还会增开一些临时航线——在飞机 转场的同时顺路捎一些乘客。JYY 研究了一下 JSOI 著名航空公司 JS Airways 的 常规直飞航线,他想计算一下,在最佳调度方案下,JS Airway 最少需要多少架 飞机才能成功执飞这所有的航线。
JSOI 王国里有 N 个机场,编号为 1 到 N。从 i 号机场到 j 号机场需要飞行Ti,j的时间。由于风向,地理位置和航空管制的因素,Ti,j 和 Tj,i并不一定相同。 
此外,由于飞机降落之后需要例行维修和加油。当一架飞机降落 k 号机场时, 需要花费Pk的维护时间才能再次起飞。 
JS Airways 一共运营 M 条航线,其中第 i 条直飞航线需要在Di时刻从Xi机场起飞,不经停,飞往Yi机场。 
为了简化问题,我们假设 JS Airway 可以在 0 时刻在任意机场布置任意多架加油维护完毕的飞机;为了减少飞机的使用数,我们允许 JS Airways 增开任意多条临时航线以满足飞机的调度需求。 
JYY 想知道,理论上 JS Airways 最少需要多少架飞机才能完成所有这 M 个 航班。

输入描述

输入一行包含两个正整数 N 和 M。
接下来一行包含 N 个正整数表示每一个机场的飞机维护时间。
接下来 N 行,每行 N 个非负整数,其中第 i 行第 j 个非负整数为Ti,j,表示 从 i 号机场飞往 j 号机场所需要花费的时间。数据保证Ti,i= 0。
接下来 M 行,每行 3 个正整数,其中第 i 行为Di",表示第 i 条航线的 起飞机场,降落机场,以及起飞时间。数据保证Xi ≠Yi, 。

输出描述

输出文件包含一行一个正整数,表示 JS Airways 理论上最少需要的飞机数。

示例1

输入:

3 3
100 1 1
0 1 1
1 0 5
2 1 0
1 2 1
2 1 1
3 1 9

输出:

2

说明:

在第一个样例中,JS Airways 可以在 0 时刻在 2 号机场安排一架飞机并执飞
第 2 条航线(2→1)。此外还需要在 0 时刻在 1 号机场安排一架飞机,这架飞机
首先执飞第 1 条航线(1→2),然后通过临时新增一条航线从 2 号机场起飞飞往
3 号机场,降落 3 号机场之后执飞第 3 条航线(3→1)。

示例2

输入:

3 3
100 1 1
0 1 1
1 0 5
2 1 0
1 2 1
2 1 1
3 1 8

输出:

3

说明:

在第二个样例中,执行完第 1 条航线的飞机无法赶上第 3 条航线的起飞时
间,因此 JS Airways 必须使用 3 架不同的飞机才能完成所有的航班。

原站题解

上次编辑到这里,代码来自缓存 点击恢复默认模板

C++14(g++5.4) 解法, 执行用时: 214ms, 内存消耗: 7264K, 提交时间: 2019-03-14 17:58:20

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 505;
const int MAXP = 1005;
const int INF = 1e9 + 10;
template <typename T> void chkmax(T &x, T y) {x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {x = min(x, y); } 
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
struct edge {int dest, flow; unsigned home; };
int n, m, v[MAXN];
int s, t, flow, dist[MAXP];
int x[MAXN], y[MAXN], d[MAXN];
int len[MAXN][MAXN], bak[MAXN][MAXN];
unsigned curr[MAXP];
vector <edge> a[MAXP];
int dinic(int pos, int limit) {
	if (pos == t) return limit;
	int used = 0, tmp;
	for (unsigned &i = curr[pos]; i < a[pos].size(); i++) {
		if (a[pos][i].flow && dist[pos] + 1 == dist[a[pos][i].dest] && (tmp = dinic(a[pos][i].dest, min(limit - used, a[pos][i].flow)))) {
			used += tmp;
			a[pos][i].flow -= tmp;
			a[a[pos][i].dest][a[pos][i].home].flow += tmp;
		}
		if (used == limit) return used;
	}
	return used;
}
bool bfs() {
	static int q[MAXP];
	static int l = 0, r = -1;
	for (int i = 0; i <= r; i++)
		dist[q[i]] = 0;
	l = 0, r = -1;
	q[++r] = s; dist[s] = 1;
	while (l <= r) {
		int tmp = q[l++];
		for (unsigned i = 0; i < a[tmp].size(); i++)
			if (a[tmp][i].flow && dist[a[tmp][i].dest] == 0) {
				dist[a[tmp][i].dest] = dist[tmp] + 1;
				q[++r] = a[tmp][i].dest;
			}
	}
	return dist[t] != 0;
}
void addedge(int s, int t, int flow) {
	a[s].push_back((edge) {t, flow, a[t].size()});
	a[t].push_back((edge) {s, 0, a[s].size() - 1});
}
int main() {
	read(n), read(m);
	for (int i = 1; i <= n; i++)
		read(v[i]);
	for (int i = 1; i <= n; i++)
	for (int j = 1; j <= n; j++) {
		read(len[i][j]);
		if (i != j) len[i][j] += v[j];
		bak[i][j] = len[i][j];
	}
	for (int i = 1; i <= n; i++)
	for (int j = 1; j <= n; j++)
	for (int k = 1; k <= n; k++)
		chkmin(len[j][k], len[j][i] + len[i][k]);
	s = 0, t = m * 2 + 1;
	for (int i = 1; i <= m; i++) {
		read(x[i]), read(y[i]), read(d[i]);
		addedge(s, i, 1);
		addedge(i + m, t, 1);
	}
	for (int i = 1; i <= m; i++)
	for (int j = 1; j <= m; j++)
		if (d[j] - d[i] >= bak[x[i]][y[i]] + len[y[i]][x[j]]) addedge(i, j + m, 1);
	while (bfs()) {
		memset(curr, 0, sizeof(curr));
		flow += dinic(s, INF);
	}
	printf("%d\n", m - flow);
	return 0;
}

C++(clang++ 11.0.1) 解法, 执行用时: 161ms, 内存消耗: 3056K, 提交时间: 2022-10-04 22:35:37

#include<bits/stdc++.h>
using namespace std;
#define N 505
struct flight {
	int u,v,t;
}q[N];
vector<int>g[N];
int n,m,ans;
int f[N][N],p[N],d[N][N],match[N];
bool vis[N];
void floyd() {
    for(int k=1;k<=n;k++) 
    	for(int i=1;i<=n;i++)
    		for(int j=1;j<=n;j++)
    			f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
}
int hungary(int u) {
	for(auto v:g[u]) {
		if(!vis[v]) {
			vis[v]=1;
			if(!match[v]||hungary(match[v])) {
				match[v]=u;
				return 1;
			}
		}
	}
	return 0;
}
int main() {
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%d",p+i);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++) {
			scanf("%d",&d[i][j]);
			f[i][j]=d[i][j]+(i==j?0:p[j]);
		}
	floyd();
	for(int i=1;i<=m;i++)
		scanf("%d%d%d",&q[i].u,&q[i].v,&q[i].t);	
	for(int i=1;i<=m;i++)
		for(int j=1;j<=m;j++) {
            int t = q[j].t - q[i].t - d[q[i].u][q[i].v] - p[q[i].v];
			if (t >= f[q[i].v][q[j].u])  g[i].push_back (j);
        }
	for(int i=1;i<=m;i++) {
    	memset(vis,0,sizeof(vis));
    	ans+=hungary(i);
    }
	cout<<m-ans;
}

上一题