[CF1641D]Two Arrays

OneInDark

已于 2022-02-24 11:40:01 修改

阅读量520

点赞数

分类专栏： C++ 数据结构 # 容斥原理文章标签：我不会做 DYM 得分

于 2022-02-24 11:20:00 首次发布

本文链接：https://blog.csdn.net/qq_42101694/article/details/123107717

版权

C++ 同时被 3 个专栏收录

547 篇文章 3 订阅

订阅专栏

数据结构

137 篇文章 1 订阅

订阅专栏

容斥原理

23 篇文章 0 订阅

订阅专栏

题目

传送门 to CF

思路

做过类似的题，有试图容斥。即 $[S=\varnothing]=\sum_{T\subseteqq S}(-1)^{|T|}$ ，二项式反演，用这个式子翻译 “交集为空”。

本来考虑 $\text{min-max}$ 反演，而后发现不对。

最后还是要 回归原问题的性质：求最小 $w_i+w_j$ 。所以我们只需找到最小的 $w_j$ 。二分后变为判定性问题，只需数个数。于是线段树上二分， $\tt set$ 维护桶，复杂度 $\mathcal O(n2^m\log^2n)$ ，空间复杂度 $\mathcal O(n2^m\log n)$ 。先是 $\rm MLE$ ，后来又 $\rm TLE$ ，确实卡不过去。

很可惜，我对原问题的性质的发掘还不足——找的是全局最小值，而非每个 $w_i$ 都必须求解最小 $w_j$ 。所以可以 双指针。时间复杂度 $\mathcal O(n2^m\log n)$ ，空间复杂度 $\mathcal O(n2^m)$ 。

然而这仍然不能通过！该死的 $\textrm{CF}$ 这都要卡是吧！必须把 $\log n$ 给搞掉。最简单的方式是 $\texttt{unordered map}$ ，只要你的哈希值足够微妙就可以通过（如果是 $\textrm{CF}$ 就要冒着被 $\rm hack$ 的风险，~~除非你使用随机哈希~~）。复杂但正确的方式是，用 $\tt trie$ 维护。

直觉上，时空复杂度皆为 $\mathcal O(n2^mm)$ 。可仔细研究后发现，将 $m$ 个元素的每个子集都插入 $\tt trie$ 是 $\mathcal O(2^m)$ 的（查询复杂度同理），即每次找到 $\tt lowbit$ ，去掉 $\tt lowbit$ 对应的子集已经在 $\tt trie$ 中，直接在其后接上 $\tt lowbit$ 就行。时间、空间都是 $\mathcal O(2^m)$ 。（可是儿子数组还是要用 $\texttt{unordered map}$ 维护……）

代码

~~为了避免代码过长，这里使用了哈希版~~。

#include <cstdio> // JZM yydJUNK!!!
#include <iostream> // XJX yyds!!!
#include <algorithm> // XYX yydLONELY!!!
#include <cstring> // (the STRONG long for LONELINESS)
#include <cctype> // ZXY yydSISTER!!!
#include <unordered_map> // I hate this motherf**ker
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
typedef long long llong;
inline int readint(){
	int a = 0, c = getchar(), f = 1;
	for(; !isdigit(c); c=getchar())
		if(c == '-') f = -f;
	for(; isdigit(c); c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}

typedef unsigned long long ullong;
const unsigned RAD = 19260817;
ullong getHash(int a[], int len){
	ullong v = 0;
	for(int i=0; i!=len; ++i)
		v = RAD*v+unsigned(a[i]);
	return v;
}

const int MAXN = 100005, MAXM = 5;
int n, m, a[MAXN][MAXM+1];

ullong s[1<<MAXM]; int len, cnt[1<<MAXM];
void got(int id){
	static int sta[MAXM], top;
	for(int S=1; S!=(1<<m); ++S){
		for(int j=top=0; j!=m; ++j) if(S>>j&1)
			sta[top ++] = a[id][j];
		s[S-1] = getHash(sta,top);
	}
}

std::unordered_map<ullong,int> mp;
void insert(){ for(int j=0; j!=len; ++j) ++ mp[s[j]]; }
void erase(){ for(int j=0; j!=len; ++j) -- mp[s[j]]; }
inline int query(){
	int now = 0;
	for(int i=0; i!=len; ++i){
		if(!mp.count(s[i])) continue;
		if(cnt[i]&1) now += mp[s[i]];
		else now -= mp[s[i]];
	}
	return now;
}

bool cmp(const int &x, const int &y){
	return a[x][m] < a[y][m];
}
int id[MAXN];
int main(){
	n = readint(), m = readint();
	rep(i,1,n) rep(j,0,m) a[i][j] = readint();
	len = (1<<m)-1; // fixed value
	for(int S=1; S!=(1<<m); ++S)
		cnt[S-1] = __builtin_popcount(S);
	rep(i,1,n){
		std::sort(a[i],a[i]+m);
		got(i), insert(); // insert all
	}
	rep(i,1,n) id[i] = i;
	std::sort(id+1,id+n+1,cmp);
	int ans = 0x7fffffff;
	for(int *i=id+1,j=n; i!=id+n+1; ++i)
		for(got(*i); query()!=j; --j,got(*i)){
			ans = std::min(ans,a[id[j]][m]+a[*i][m]);
			got(id[j]), erase(); // remove
		}
	if(ans == 0x7fffffff) ans = -1;
	printf("%d\n",ans);
	return 0;
}