BZOJ3671: [Noi2014]随机数生成器(贪心)

Time Limit: 50 Sec  Memory Limit: 256 MB
Submit: 2098  Solved: 946
[Submit][Status][Discuss]

Description

 

Input

第1行包含5个整数，依次为 x_0,a,b,c,d ，描述小H采用的随机数生成算法所需的随机种子。第2行包含三个整数 N,M,Q ，表示小H希望生成一个1到 N×M 的排列来填入她 N 行 M 列的棋盘，并且小H在初始的 N×M 次交换操作后，又进行了 Q 次额外的交换操作。接下来 Q 行，第 i 行包含两个整数 u_i,v_i，表示第 i 次额外交换操作将交换 T_(u_i )和 T_(v_i ) 的值。

Output

输出一行，包含 N+M-1 个由空格隔开的正整数，表示可以得到的字典序最小的路径序列。

Sample Input

1 3 5 1 71
3 4 3
1 7
9 9
4 9

Sample Output

1 2 6 8 9 12

HINT

本题的空间限制是 256 MB，请务必保证提交的代码运行时所使用的总内存空间不超过此限制。

一个32位整数（例如C/C++中的int和Pascal中的Longint）为4字节，因而如果在程序中声明一个长度为 1024×1024 的32位整型变量的数组，将会占用 4 MB 的内存空间。

2≤N,M≤5000

0≤Q≤50000

0≤a≤300

0≤b,c≤108

0≤x0<d≤1081≤ui,vi≤N×

Source

这题究竟是几个意思？？？

为了科普随机化算法？？(雾)

字典序最小？那不就是贪心走小的就行了么？

顺便维护一下每一个$x$，对应那些$y$是能选的

选了一个数之后它左下和右上的矩阵就都不能选了

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<ctime>
#include<cstdlib>
#define LL long long
using namespace std;
const int MAXN = 1e5 + , INF = 1e8 + , mod = 1e9 + , mod2 = 1e9 + ;
inline int read() {
    char c = getchar(); int x = , f = ;
    while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x * f;
}
int N, M, Q, tot = ;
LL Xi, X0, a, b, c, d;
int A[ * ], Map[][], l[MAXN], r[MAXN];
int xx[] = {, +, }, yy[] = {, , +};
main() {
#ifdef WIN32
    ///freopen("a.in", "r", stdin);
#endif
    X0 = read(); a = read(), b = read(), c = read(), d = read();
    N = read(); M = read(); Q = read();
    for(int i = ; i <= N * M; i++) A[i] = i;
    for(int i = ; i <= N * M; i++) swap(A[i], A[(Xi = (a * X0 * X0 + b * X0 + c) % d) % i + ]), X0 = Xi;

    //for(int i = 1; i <= N * M; i++)    printf("%d\n", A[i]);
    while(Q--) {int x = read(), y = read(); swap(A[x], A[y]);}
    for(int i = ; i <= N; i++) for(int j = ; j <= M; j++) Map[i][j] = A[++tot];
    for(int i = ; i <= N; i++)
        for(int j = ; j <= M; j++)
            A[Map[i][j]] =(i - ) * M + j;
    memset(Map, , sizeof(Map));

    for(int i = ; i <= N; i++) l[i] = , r[i] = M + ;
    for(int i = ; i <= N * M; i++) {
        int y = A[i] % M, x = A[i] / M + ;
        if(y == ) x--, y = M;
        if(y <= l[x] || y >= r[x]) continue;
        printf("%d ", i);
        for(int j = x - ; j >= ; j--) r[j] = min(y + , r[j]);
        for(int j = x + ; j <= N; j++) l[j] = max(y - , l[j]);
    }
    return ;
}