CDQ分治2

CDQ套CDQ：四维偏序问题

题目来源：COGS 2479 偏序

#define LEFT 0
#define RIGHT 1struct Node{int a,b,c,d,bg;};
Node q[_],tmp1[_],tmp2[_];  int aa,bb,cc,dd,n; long long Ans;void cdq2(RG int L,RG int R){if(L == R)return;RG int mid = (L+R)>>1; cdq2(L,mid); cdq2(mid+1,R);RG int oo = L , l = L , r = mid + 1;while(l <= mid && r <= R){if(tmp1[l].c < tmp1[r].c){if(tmp1[l].bg == LEFT)BIT:: Update(tmp1[l].d);tmp2[oo++] = tmp1[l++];}else {if(tmp1[r].bg == RIGHT)Ans += BIT::Query(tmp1[r].d);tmp2[oo++] = tmp1[r++];}}while(l <= mid)tmp2[oo++] = tmp1[l++];while(r <= R){if(tmp1[r].bg == RIGHT)Ans += BIT::Query(tmp1[r].d) ;tmp2[oo++] = tmp1[r++];}for(RG int i = L; i <= R; i ++){if(tmp2[i].bg == LEFT)BIT::Clear(tmp2[i].d);tmp1[i] = tmp2[i];}return;
}void cdq1(int L,int R){if(L == R)return;RG int mid = (L+R)>>1; cdq1(L,mid);  cdq1(mid+1,R);RG int oo = L , l = L ,r = mid + 1;while(l <= mid && r <= R){if(q[l].b < q[r].b)q[l].bg = LEFT ,tmp1[oo++] = q[l++] ;else q[r].bg = RIGHT , tmp1[oo++] = q[r++] ;}while(l <= mid)q[l].bg = LEFT , tmp1[oo++] = q[l++];while(r <= R)q[r].bg = RIGHT , tmp1[oo++] = q[r++];for(RG int i = L; i <= R; i ++)q[i] = tmp1[i];cdq2(L,R);
}int main(){freopen("partial_order.in","r",stdin);freopen("partial_order.out","w",stdout);n = gi();for(RG int i = 1; i <= n; i ++)q[i].a = i; for(RG int i = 1; i <= n; i ++)q[i].b = gi(); for(RG int i = 1; i <= n; i ++)q[i].c = gi(); for(RG int i = 1; i <= n; i ++)q[i].d = gi(); Ans = 0; cdq1(1,n);cout<<Ans;  return 0;
}

CDQ套CDQ套CDQ：五维偏序问题

题目来源：COGS 2580 偏序 \(II\)

#define LEFT 0
#define RIGHT 1struct Node{int d1,d2,d3,d4,d5,bg[2];};
Node q[_] , tmp1[_] , tmp2[_] , tmp3[_];  int n; long long ans;IL bool Check(RG int id,RG int opt)
{return tmp2[id].bg[0]==opt && tmp2[id].bg[1]==opt;}void cdq3(RG int L,RG int R){      //cdq分治解决三维偏序问题 if(L == R)return;RG int mid = (L+R)>>1; cdq3(L,mid); cdq3(mid+1,R);RG int l = L , r = mid+1 , oo = L-1;while(l <= mid && r <= R){if(tmp2[l].d4 < tmp2[r].d4){if(Check(l,LEFT))BIT::Update(tmp2[l].d5);tmp3[++oo] = tmp2[l++];}else{if(Check(r,RIGHT))ans += BIT::Query(tmp2[r].d5);tmp3[++oo] = tmp2[r++];}}while(l <= mid)tmp3[++oo] = tmp2[l++];while(r <= R){if(Check(r,RIGHT))ans += BIT::Query(tmp2[r].d5);tmp3[++oo] = tmp2[r++];}for(RG int i = L; i <= R; i ++){if(Check(i,LEFT))BIT::Clear(tmp2[i].d5);tmp2[i] = tmp3[i];}return;}void cdq2(RG int L,RG int R){       //消去第三维的影响 if(L == R)return;RG int mid = (L+R)>>1; cdq2(L,mid); cdq2(mid+1,R);RG int l = L , r = mid+1 , oo = L-1;while(l <= mid && r <= R){if(tmp1[l].d3 < tmp1[r].d3)tmp1[l].bg[1] = LEFT , tmp2[++oo] = tmp1[l++];elsetmp1[r].bg[1] = RIGHT , tmp2[++oo] = tmp1[r++];}while(l <= mid)tmp1[l].bg[1] = LEFT , tmp2[++oo] = tmp1[l++];while(r <= R)tmp1[r].bg[1] = RIGHT , tmp2[++oo] = tmp1[r++];for(RG int i = L; i <= R; i ++)tmp1[i] = tmp2[i];cdq3(L,R); }void cdq1(RG int L,RG int R){         //消去第二维的影响 if(L == R)return;RG int mid = (L+R)>>1; cdq1(L,mid); cdq1(mid+1,R);RG int l = L , r = mid+1 , oo = L-1;while(l <= mid && r <= R){if(q[l].d2 < q[r].d2)q[l].bg[0] = LEFT , tmp1[++oo] = q[l++];else q[r].bg[0] = RIGHT , tmp1[++oo] = q[r++];} while(l<=mid)q[l].bg[0] = LEFT , tmp1[++oo] = q[l++];while(r<=R)q[r].bg[0] = RIGHT , tmp1[++oo] = q[r++];for(RG int i = L; i <= R; i ++)q[i] = tmp1[i];      cdq2(L,R);
}int main(){freopen("partial_order_two.in","r",stdin);freopen("partial_order_two.out","w",stdout);n = gi();for(RG int i = 1; i <= n; i ++)q[i].d1 = i;for(RG int i = 1; i <= n; i ++)q[i].d2 = gi();for(RG int i = 1; i <= n; i ++)q[i].d3 = gi();for(RG int i = 1; i <= n; i ++)q[i].d4 = gi();for(RG int i = 1; i <= n; i ++)q[i].d5 = gi();//sort_a  消去第一维的影响 ans = 0; cdq1(1,n);cout << ans;  return 0;
}