bzoj2253

cdq分治+dp

看见三维偏序是cdq，互相包含是最长上升子序列

这个代码是错的

交了两份代码，发现手动出数据是不一样的。。。

不调了

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = ;
int n, lim, ans;
struct BIT {
    int tree[N];
    int lowbit(int i)
    {
        return i & (-i);
    }
    void update(int pos, int delta)
    {
        for(int i = pos; i <= lim; i += lowbit(i)) tree[i] = max(tree[i], delta);
    }
    int query(int pos)
    {
        int ret = ;
        for(int i = pos; i; i -= lowbit(i)) ret = max(ret, tree[i]);
        return ret;
    }
    void clr(int pos)
    {
        for(int i = pos; i <= lim; i += lowbit(i)) tree[i] = ;
    }
} tr;
struct data {
    ll x, y, z;
    int type, id;
} f[N], g[N];
ll a, p;
map<ll, int> mp;
vector<ll> vt;
ll t[], c[];
int dp[N];
bool cpx(data A, data B)
{
    if(A.x == B.x)
    {
        if(A.y == B.y) return A.z < B.z;
        else return A.y < B.y;
    }
    return A.x < B.x;
}
bool cpy(data A, data B)
{
    return A.y == B.y ? A.z < B.z : A.y < B.y;
}
void cdq(int l, int r)
{
    if(l >= r) return;
    int mid = (l + r) >> ;
    while(f[mid].x == f[mid - ].x) --mid;
    if(l > mid) return;
    cdq(l, mid);
    for(int i = l; i <= r; ++i) g[i] = f[i];
    sort(f + l + , f + mid + , cpy);
    sort(f + mid + , f + r + , cpy);
    int i = l, j = mid + ;
    while(j <= r)
    {
        while(i <= mid && f[i].y < f[j].y)
        {
            tr.update(f[i].z, dp[f[i].id]);
            ++i;
        }
        dp[f[j].id] = max(dp[f[j].id], tr.query(f[j].z - ) + );
        ans = max(ans, dp[f[j].id]);
        ++j;
    }
    for(int i = l; i <= mid; ++i)
        tr.clr(f[i].z);
    for(int i = l; i <= r; ++i) f[i] = g[i];
    cdq(mid + , r);
}
int main()
{
    scanf("%lld%lld%d", &a, &p, &n);
    ll p1 = a % p * a % p * a % p;
    f[].id = ;
    t[] = a % p;
    t[] = a % p * a % p;
    t[] = a % p * a % p * a % p;
    for(int i = ; i <= ; ++i) c[i] = t[i];
    sort(c + , c + );
    f[].x = c[] * p1 % p;
    f[].y = c[] * p1 % p;
    f[].z = c[] * p1 % p;
    vt.push_back(f[].z);
    for(int i = ; i <= n; ++i)
    {
        t[] = t[] * p1 % p;
        t[] = t[] * p1 % p;
        t[] = t[] * p1 % p;
        for(int j = ; j <= ; ++j) c[j] = t[j];
        sort(c + , c + );
        f[i].x = c[];
        f[i].y = c[];
        f[i].z = c[];
        f[i].id = i;
        vt.push_back(f[i].z);
    }
    sort(vt.begin(), vt.end());
    vt.erase(unique(vt.begin(), vt.end()), vt.end());
    for(int i = ; i < vt.size(); ++i) mp[vt[i]] = i + ;
    for(int i = ; i <= n; ++i) f[i].z = mp[f[i].z];
    lim = vt.size();
    sort(f + , f + n + , cpx);
    cdq(, n);
    printf("%d\n", ans);
     return ;
}