poj 2828 块状链表 OR 线段树 OR 树状数组

这个题是真正可以体现出块状链表的优点。数组定位快但是插入慢，而链表插入快却定位慢。块状链表正是结合了数组和链表的优点将定位和插入的复杂度变成了sqrt(n)。

块状链表：

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int N = 750;
int b[N][N];
int sz[N];
int next[N];
int ps;
int bl;

void init()
{
    bl = 1;
    ps = 550;
    memset( sz, 0, sizeof(sz) );
    memset( next, -1, sizeof(next) );
}

void spilt( int cur )
{
    int tmp = next[cur];
    int ncur = next[cur] = bl++;
    next[ncur] = tmp;
    for ( int i = sz[cur] / 2; i < sz[cur]; i++ )
    {
        b[ncur][sz[ncur]++] = b[cur][i];
    }
    sz[cur] = sz[cur] / 2;
}

void insert( int pos, int val )
{
    int cur = 0, p = sz[cur];
    while ( p < pos + 1 && next[cur] != -1 )
    {
        cur = next[cur];
        p += sz[cur];
    }
    if ( p < pos + 1 )
    {
        b[cur][sz[cur]++] = val;
    }
    else
    {
        p -= sz[cur];
        pos -= p;
        for ( int j = sz[cur] - 1; j >= pos; j-- )
        {
            b[cur][j + 1] = b[cur][j];
        }
        b[cur][pos] = val;
        sz[cur]++;
    }
    if ( sz[cur] > ps ) spilt(cur);
}

int main ()
{
    int n;
    while ( scanf("%d", &n) != EOF )
    {
        init();
        for ( int i = 0; i < n; i++ )
        {
            int pos, val;
            scanf("%d%d", &pos, &val);
            insert( pos, val );
        }
        int cnt = 0;
        for ( int i = 0; i != -1; i = next[i] )
        {
            for ( int j = 0; j < sz[i]; j++ )
            {
                printf("%d", b[i][j]);
                cnt++;
                if ( cnt != n )
                {
                    putchar(' ');
                }
                else
                {
                    putchar('
');
                }
            }
        }
    }
    return 0;
}

线段树：

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int N = 200001;
int ans[N];
int n;

struct X
{
    int pos, val;
} x[N];

struct Node 
{
    int l, r, sum;
} node[N << 2];

void pushup( int i )
{
    node[i].sum = node[i << 1].sum + node[i << 1 | 1].sum;
}

void build( int i, int l, int r )
{
    node[i].l = l, node[i].r = r, node[i].sum = node[i].r - node[i].l + 1;
    if ( l == r ) return ;
    int mid = ( l + r ) >> 1;    
    build( i << 1, l, mid );
    build( i << 1 | 1, mid + 1, r );
}

void update( int i, int pos, int val )
{
    if ( node[i].l == pos && node[i].r == pos )
    {
        node[i].sum = val;
        return ;
    }
    int mid = ( node[i].l + node[i].r ) >> 1;
    if ( pos <= mid )
    {
        update( i << 1, pos, val );
    }
    else
    {
        update( i << 1 | 1, pos, val );
    }
    pushup(i);
}

int query( int i, int v )
{
    if ( node[i].l == node[i].r ) return node[i].l;
    if ( node[i << 1].sum >= v ) return query( i << 1, v );
    return query( i << 1 | 1, v - node[i << 1].sum );
}

int main ()
{
    while ( scanf("%d", &n) != EOF )
    {
        for ( int i = 1; i <= n; i++ )
        {
            scanf("%d%d", &x[i].pos, &x[i].val);
            x[i].pos++;
        }
        build( 1, 1, n );
        for ( int i = n; i >= 1; i-- )
        {
            int tmp = query( 1, x[i].pos );
            ans[tmp] = x[i].val;
            update( 1, tmp, 0 );
        }    
        for ( int i = 1; i < n; i++ )
        {
            printf("%d ", ans[i]);
        }    
        printf("%d
", ans[n]);
    }
    return 0;
}

树状数组：

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int N = 200001;
int c[N];
int ans[N];
int n;

struct X
{
    int pos, val;
} x[N];

int lb( int i )
{
    return i & -i;
}

void update( int i, int v )
{
    while ( i <= n )
    {
        c[i] += v;
        i += lb(i);
    }
}

int query( int k )
{
    int ans = 0, cnt = 0;
    for ( int i = 17; i >= 0; i-- )
    {
        ans += ( 1 << i );
        if ( ans >= n || cnt + c[ans] >= k )
        {
            ans -= ( 1 << i );
        }
        else
        {
            cnt += c[ans];
        }
    }
    return ans + 1;
}

int main ()
{
    while ( scanf("%d", &n) != EOF )
    {
        memset( c, 0, sizeof(c) );
        for ( int i = 1; i <= n; i++ )
        {
            scanf("%d%d", &x[i].pos, &x[i].val);
            x[i].pos++;
            update( i, 1 );
        }
        for ( int i = n; i >= 1; i-- )
        {
            int tmp = query( x[i].pos );
            ans[tmp] = x[i].val;
            update( tmp, -1 );
        }    
        for ( int i = 1; i < n; i++ )
        {
            printf("%d ", ans[i]);
        }    
        printf("%d
", ans[n]);
    }
    return 0;
}