Cod sursa(job #133943)

Utilizator dominoMircea Pasoi domino Data 10 februarie 2008 00:33:24
Problema Eprubeta Scor 5
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.01 kb
#include <stdio.h>
#include <stdlib.h>
#include <algorithm>

using namespace std;

#define MAX_N 2048
#define MAX_LIST 4000005
#define FIN "eprubeta.in"
#define FOUT "eprubeta.out"
#define SQRT 10
#define NIL 0

int val0[MAX_LIST], val1[MAX_LIST], next0[MAX_LIST], next1[MAX_LIST], list_idx;
int N, M, A, B, Z, L[MAX_N], R[MAX_N], V[MAX_N], C[MAX_N];

inline int find(int L, int idx)
{
    int p, t;

    if (idx <= 0) return L;
    t = idx/SQRT; idx %= SQRT;
    for (p = L; t; --t, p = next1[p]);
    for (; idx; --idx, p = next0[p]);
    return p;
}

inline void update(int &L, int end)
{
    int p, cnt = SQRT, t = 0;

    for (p = L; p && cnt; --cnt, p = next0[p])
        t += val0[p];
    val1[L] = t;
    next1[L] = p;
    for (p = L, t = next0[L]; p != end; p = next0[p])
    {
        next1[t] = next0[next1[p]];
        val1[t] = val1[p]-val0[p]+val0[next1[p]];
    }
}

void make_list(int a, int b, int &L)
{
    int i, p;

    if (a >= b) return;
    for (i = b-1; i >= a; --i)
    {
        p = ++list_idx;
        val0[p] = C[i];
        val1[p] = 0;
        next0[p] = L;
        next1[p] = NIL;
        L = p;
    }
    update(L, NIL);
}


inline void rotate(int &L1, int &L2, int size)
{
    int p1, q1, p2, q2;

    if (size <= 0) return;
    p1 = find(L1, size-1);
    q1 = find(L1, size-SQRT);
    p2 = find(L2, size-1);
    q2 = find(L2, size-SQRT);
    swap(next0[p1], next0[p2]);
    swap(L1, L2);
    update(q1, p1);
    update(q2, p2);
}

inline int get_sum(int L, int idx)
{
    int p, sum = 0, t = idx/SQRT;

    idx %= SQRT;
    for (p = L; t; --t, p = next1[p])
        sum += val1[p];
    for (; idx; --idx, p = next0[p])
        sum += val0[p];
    return sum;
}

int main(void)
{
    int i, j, k, x;

    freopen(FIN, "r", stdin);
    freopen(FOUT, "w", stdout);

    scanf("%d %d %d %d %d", &N, &M, &Z, &A, &B);

    for (i = 0; i < N; ++i)
    {
        for (j = 0; j < N; ++j)
        {
            C[j] = ((i+A)*(j+B)/4) % Z;
            if (j == i) V[i] = C[j];
        }
        reverse(C, C+i);
        make_list(0, i, L[i]);
        make_list(i+1, N, R[i]);
    }

    for (; M; --M)
    {
        scanf("%d %d %d", &x, &i, &j);
        if (x == 1)
        {
            for (x = k = j-i; i <= j; ++i, --j, --k)
            {
                rotate(L[j], R[i], k);
                if (i != j) rotate(L[i], R[j], x-k);
                swap(V[i], V[j]);
            }
        }
        else 
        {
            int sum = 0;
            for (x = k = j-i; i <= j; ++i, --k)
            {
                sum += get_sum(L[i], x-k); 
                sum += get_sum(R[i], k);
                sum += V[i];
            }
            printf("%d\n", sum);
        }
    }

    unsigned total_sum = 0;
    for (i = 0; i < N; ++i)
    {
        x = V[i]+get_sum(L[i], i)+get_sum(R[i], N-1-i);
        total_sum += (unsigned)x*(unsigned)x*(unsigned)(i+1);
    }
    printf("%u\n", total_sum);

    return 0;
}