CCF计算机职业资格认证考试题解系列文章为meelo原创,请务必以链接形式注明本文地址

CCF CSP 201703-5 引水入城

问题描述

  MF城建立在一片高原上。由于城市唯一的水源是位于河谷地带的湖中,人们在坡地上修筑了一片网格状的抽水水管,以将湖水抽入城市。如下图所示:

  这片管网由 n 行 m 列节点(红色,图中 n = 5,m = 6),横向管道(紫色)和纵向管道(橙色)构成。
  行和列分别用 1 到 n 的整数和 1 到 m 的整数表示。第 1 行的任何一个节点均可以抽取湖水,湖水到达第 n 行的任何一个节点即算作引入了城市。
  除第一行和最后一行外,横向相邻或纵向相邻的两个节点之间一定有一段管道,每一段管道都有各自的最大的抽水速率,并需要根据情况选择抽水还是放水。对于纵向的管道(橙色),允许从上方向下方抽水或从下方向上方放水;如果从图中的上方向下方抽水,那么单位时间内能通过的水量不能超过管道的最大速率;如果从下方向上方放水,因为下方海拔较高,因此可以允许有任意大的水量。对于横向的管道(紫色),允许从左向右或从右向左抽水,不允许放水,两种情况下单位时间流过的水量都不能超过管道的最大速率。
  现在MF城市的水务负责人想知道,在已知每个管道单位时间容量的情况下,MF城每单位时间最多可以引入多少的湖水。

输入格式

  由于输入规模较大,我们采用伪随机生成的方式生成数据。
  每组数据仅一行包含 6 个非负整数 nmABQX0。其中,n 和 m 如前文所述,表示管网的大小,保证 2 ≤ nm ≤ 5000;保证 1 ≤ ABQX0 ≤ 109
  ABQX0 是数据生成的参数,我们用如下的方式定义一个数列 { Xi }:
  Xi+1 = ( AXi + Bmod Q, (i ≥ 0)
  我们将数列的第 1 项到第 (n-1)m 项作为纵向管道的单位时间容量,其中 X(s-1)m+t 表示第 s 行第 t 列的节点到第 s+1 行第 t 列管道单位时间的容量;将数列的第 (n-1)m+1 项到第 (n-1)m+(n-2)(m-1) 项(即接下来的 (n-2)(m-1) 项)作为横向管道的单位时间容量,其中 X(n-1)m+(s-2)(m-1)+t 表示第 s 行第 t 列的节点到第 s 行第 t+1 列管道单位时间的容量。

输出格式

  输出一行一个整数,表示MF城每单位时间可以引入的水量。
  注意计算过程中有些参数可能超过32位整型表示的最大值,请注意使用64位整型存储相应数据。

样例输入

3 3 10 3 19 7

样例输出

38

样例说明

  使用参数得到数列 { Xi }={ 7, 16, 11, 18, 12, 9, 17, 2, 4, … },按照输入格式可以得到每个管道的最大抽水量如下图所示:

  在标准答案中,单位时间可以引水 38 单位。所有纵向管道均向下抽水即可,不需要横向管道抽水,也不需要向上放水。

样例输入

2 5 595829232 749238243 603779819 532737791

样例输出

1029036148

样例输入

5 2 634932890 335818535 550589587 977780683

样例输出

192923706

样例输入

5 5 695192542 779962396 647834146 157661239

样例输出

1449991168

评测用例规模与约定

  共有10组评测数据,每组数据的参数规模如下所示:

测试点编号 n m
1 =2 =1000
2 =1000 =2
3 =1000 =2
4 =5 =5
5 =10 =10
6 =100 =100
7 =500 =500
8 =1000 =1000
9 =2000 =2000
10 =5000 =5000

解析

这是一个最大流的问题,湖是源,城市是汇。

下面实现了ford-fulkerson算法,只能通过50%的数据。

有更好的方法求告知!

代码

C++

#include <iostream>
#include <vector>
#include <queue>
#include <utility>
#include <climits>
using namespace std; long long A, B, Q, X;
int numVertex; int nextRandom() {
X = (A * X + B) % Q;
return X;
} struct Edge {
int v; // vertex
int w; // weight
Edge(int v_, int w_) : v(v_), w(w_) {}
}; bool bfs(vector<vector<Edge> > &rgraph, int s, int t, vector<pair<int,int> > &parents) {
queue<int> q;
vector<bool> visited(numVertex);
q.push(s);
while(!q.empty()) {
int u = q.front();
q.pop();
for(int i=; i<rgraph[u].size(); i++) {
int v = rgraph[u][i].v;
if(!visited[v] && rgraph[u][i].w>) {
visited[v] = true;
parents[v] = make_pair(u,i);
q.push(v);
if(v == t) return true;
} }
}
return false;
} long long fordFulkerson(vector<vector<Edge> > &rgraph, int s, int t) {
long long maxFlow = ;
vector<pair<int,int> > parents(numVertex);
while(bfs(rgraph, , , parents)) {
int pathFlow = INT_MAX;
for(int v=t; v!=s; ) {
int u=parents[v].first; int ui = parents[v].second;
pathFlow = min(pathFlow, rgraph[u][ui].w);
v = u;
}
maxFlow += pathFlow;
// cout << pathFlow << " " << maxFlow << endl;
for(int v=t; v!=s; ) {
int u = parents[v].first;
int ui = parents[v].second;
if(rgraph[u][ui].w!=INT_MAX) rgraph[u][ui].w -= pathFlow;
int vi = -;
for(int i=; i<rgraph[v].size(); i++) {
if(rgraph[v][i].v == u) {
vi = i;
}
}
if(vi!=- && rgraph[v][vi].w!=INT_MAX) rgraph[v][vi].w += pathFlow;
v = u;
}
}
return maxFlow;
} int main() {
int N, M;
cin >> N >> M >> A >> B >> Q >> X;
numVertex = N * M + ;
// 0:source, 1:sink,
vector<vector<Edge> > graph(numVertex, vector<Edge>()); int offset = ;
// construct graph
for(int n=; n<N-; n++) {
for(int m=; m<M; m++) {
int from = n*M+m+offset;
int to = from+M;
nextRandom();
graph[from].push_back(Edge(to, X));
graph[to].push_back(Edge(from, INT_MAX));
}
} for(int m=; m<M; m++) {
int from = ;
int to = m+offset;
graph[from].push_back(Edge(to, INT_MAX));
} for(int m=; m<M; m++) {
int from = (N-)*M+m+offset;
int to = ;
graph[from].push_back(Edge(to, INT_MAX));
} long long maxFlow = ; for(int n=; n<N-; n++) {
for(int m=; m<M-; m++) {
int from = n*M+m+offset;
int to = from+;
nextRandom();
graph[from].push_back(Edge(to, X));
graph[to].push_back(Edge(from, X));
}
} maxFlow += fordFulkerson(graph, , ); cout << maxFlow;
}

CCF CSP 201703-5 引水入城(50分)的更多相关文章

  1. CCF(引水入城:60分):最大流+ISAP算法

    引水入城 201703-5 这从题目分析来看很像最大流的问题,只需要增加一个超级源点和一个超级汇点就可以按照题意连边再跑最大流算法. 因为数据量太大了,肯定会超时.但是没有想到可行的解决方法. #in ...

  2. Codevs 1066 引水入城 2010年NOIP全国联赛提高组

    1066 引水入城 2010年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 传送门 题目描述 Description 在一个遥远的国度 ...

  3. CODEVS 1066/洛谷 P1514引水入城

    1066 引水入城 2010年NOIP全国联赛提高组  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 钻石 Diamond   题目描述 Description 在一个遥远的国 ...

  4. vijos p1777 引水入城(bfs+贪心)

    引水入城   描述 在一个遥远的国度,一侧是风景秀美的湖泊,另一侧则是漫无边际的沙漠.该国的行政区划十分特殊,刚好构成一个N行M列的矩形,其中每个格子都代表一座城市,每座城市都有一个海拔高度. 为了使 ...

  5. 521. [NOIP2010] 引水入城 cogs

    521. [NOIP2010] 引水入城 ★★★   输入文件:flow.in   输出文件:flow.out   简单对比时间限制:1 s   内存限制:128 MB 在一个遥远的国度,一侧是风景秀 ...

  6. 引水入城 2010年NOIP全国联赛提高组(bfs+贪心)

    1066 引水入城 2010年NOIP全国联赛提高组  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 钻石 Diamond     题目描述 Description 在一个遥远 ...

  7. NOIP2010 引水入城

    4引水入城 题目描述 在一个遥远的国度,一侧是风景秀美的湖泊,另一侧则是漫无边际的沙漠.该国的行政区划十分特殊,刚好构成一个N 行M 列的矩形,如上图所示,其中每个格子都代表一座城市,每座城市都有一个 ...

  8. Luogu 1514 引水入城 (搜索,动态规划)

    Luogu 1514 引水入城 (搜索,动态规划) Description 在一个遥远的国度,一侧是风景秀美的湖泊,另一侧则是漫无边际的沙漠.该国的行政区划十分特殊,刚好构成一个N行M列的矩形,如上图 ...

  9. 洛谷P1514 引水入城

    洛谷P1514 引水入城 原题链接 一道好题...细节真多 第一次提交90分,然后就GG了,不知从何改起 其实比较简单吧... 首先,一个点的水流向最后一排,一定可以形成一个区间. 不行的话肯定GG ...

随机推荐

  1. django中的认证与登录

    认证登录 django.contrib.auth中提供了许多方法,这里主要介绍其中的三个: 1  authenticate(**credentials)    提供了用户认证,即验证用户名以及密码是否 ...

  2. C语言复习---判断素数

    一般 int main01() { ; scanf("%d", &a); n_sqrt = sqrt(a); ; i <= n_sqrt; i++) ) { flag ...

  3. JAVA-JSP隐式对象

    JSP隐式对象 在本章中,我们将讨论和学习JSP中的隐式对象.这些对象是JSP容器为每个页面中的开发人员提供的Java对象,开发人员可以直接调用它们而不用显式地声明它们再调用. JSP隐式对象也称为预 ...

  4. tips 前端 bootstrap 嵌套行 嵌套列 溢出 宽度不正确 栅格化系统计算

    bootstrap 当嵌套列时 有时会出现很奇异的row 的width不对问题出现的情况时 <div class="row" > <!--row a--> ...

  5. 2015/11/7用Python写游戏,pygame入门(7):碰撞检测

    我们已经完成了飞机大战的大部分东西,但是游戏还是没有办法正式开玩,因为子弹并不能打掉飞机.只有完成了这一个工作,游戏才算基本成型. 今天的内容就非常简单了,就是做到这个碰撞检测,以及控制好子弹和飞机的 ...

  6. [Luogu 3128] USACO15DEC Max Flow

    [Luogu 3128] USACO15DEC Max Flow 最近跟 LCA 干上了- 树剖好啊,我再也不想写倍增了. 以及似乎成功转成了空格选手 qwq. 对于每两个点 S and T,求一下 ...

  7. Linux_用户管理&权限管理

    2017年1月11日, 星期三 Linux_用户管理&权限管理 1.  Linux用户管理&权限管理 终端的概念: tty  查看登录的终端 类型  user group    oth ...

  8. Eclipse升级后导入插件的方法

    昨天把eclipse 从 indigo(3.7) 升级到了 juno(4.2). 启动 eclipse 发现从前装的 android 和 maven 插件都不见了. $ ls .eclipse/ or ...

  9. xgboost与sklearn的接口

    xgb使用sklearn接口(推荐) XGBClassifier from xgboost.sklearn import XGBClassifier clf = XGBClassifier( sile ...

  10. G - Pandaland HDU - 6005 (找最小环)

    题目链接:https://cn.vjudge.net/contest/275153#problem/G 具体思路: 我们可以按照暴力的方法进行做 , 我们可以枚举每一条边,将这条边的权值设置为inf, ...