HihoCoder - 1513 bitset處理五維偏序
阿新 • • 發佈:2018-08-06
時空復雜度 偏序 題目 get 時空 維護 hoc mat cin 個關鍵字的第\(i\)名前面的集合狀態
題意:給出\(n<3e4\)個有序組\((a,b,c,d,e)\),求對第\(i\)個有序組有多少個\(j\)滿足\((a_j<a_i,b_j<b_i,c_j<c_i,d_j<d_i,e_j<e_i)\)
五維偏序問題按套路來可以排序+樹套樹套樹套樹(打死
然而這是顯然連\(O(n^2)\)暴力都不如的
可是題目給4s,\(O(n^2)\)是不可能的,但在神奇的bitset加持下\(O(5*n^2/32)\)的時空復雜度是可以卡過去的!
用bitset表示集合,\(bit[i]:\)如果\(i\)在集合中就設為1,否則0
維護\(bit[i][j]\),表示排在第\(j\)
如果\(rank[k]<rank[i]\),則置\(bit[i][j][k]=1\)
這裏如果是有序的就可以直接利用bitset的\(O(n^2/32)\)構造進行暴力傳遞
然後對5個關鍵字的集合求交就可以得出各個答案
#include<bits/stdc++.h> #include<unordered_set> #include<unordered_map> #define rep(i,j,k) for(register int i=j;i<=k;i++) #define rrep(i,j,k) for(register int i=j;i>=k;i--) #define erep(i,u) for(register int i=head[u];~i;i=nxt[i]) #define iter(i,j) for(int i=0;i<(j).size();i++) #define print(a) printf("%lld",(ll)a) #define println(a) printf("%lld\n",(ll)a) #define printbk(a) printf("%lld ",(ll)a) #define IOS ios::sync_with_stdio(0) using namespace std; const int MAXN = 3e4+11; const int oo = 0x3f3f3f3f; typedef long long ll; ll read(){ ll x=0,f=1;register char ch=getchar(); while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } int a[MAXN][6],ra[MAXN][6]; bitset<MAXN> bs[MAXN][6],t; int main(){ int n; while(cin>>n){ rep(i,0,n) rep(j,1,5) bs[i][j].reset(); rep(i,1,n) rep(j,1,5) ra[a[i][j]=read()][j]=i; rep(i,2,n) rep(j,1,5){ bs[i][j]=bs[i-1][j]; bs[i][j][ra[i-1][j]]=1; } rep(i,1,n){ t=bs[a[i][1]][1]; rep(j,2,5) t&=bs[a[i][j]][j]; println(t.count()); } } return 0; }
HihoCoder - 1513 bitset處理五維偏序