【Codeforces】1119 Global Round 2 C-H简要题解

最新推荐文章于 2020-01-28 00:32:45 发布

ccosi

最新推荐文章于 2020-01-28 00:32:45 发布

阅读量1.8k

点赞数 5

分类专栏： -------DP------- 树形DP -------数据结构------- 堆 -------算法------- 构造 -------MATH------- ---多项式---

本文链接：https://blog.csdn.net/corsica6/article/details/89073379

版权

-------数据结构------- 同时被 3 个专栏收录

127 篇文章 0 订阅

订阅专栏

-------DP-------

114 篇文章 1 订阅

订阅专栏

-------算法-------

95 篇文章 0 订阅

订阅专栏

C.Ramesses and Corner Inversion

将 $A[i][j]\ xor \ B[i][j]$ ，转成了判断 $A$ 是否能构成0矩阵。

发现无论如何操作，每行每列1的个数的奇偶性不变。所以若每行每列均有偶数个1，则必然有解(感性证明)。

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first 
#define sc second
using namespace std;
typedef double db;
typedef long long ll;
const int N=505;

int n,m;
int a[N][N],b[N][N];


inline void fl()
{
	puts("No");exit(0);
}

int main(){
	int i,j,x,y;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;++i)
	 for(j=1;j<=m;++j)
	  scanf("%d",&a[i][j]);
	for(i=1;i<=n;++i)
	 for(j=1;j<=m;++j)
	  scanf("%d",&b[i][j]);
	for(i=1;i<=n;++i){
		for(x=0,j=1;j<=m;++j)
		 if(a[i][j]^b[i][j]) x^=1;
		if(x) fl();
	} 
	for(i=1;i<=m;++i){
		for(x=0,j=1;j<=n;++j)
		 if(a[j][i]^b[j][i]) x^=1;
		if(x) fl();
	}
	puts("Yes");
	return 0;
}

D.Frets On Fire

先把 $s_i$ 去重，设去重后共有 $t$ 个数。

设 $l e n = r - l + 1$ ， $l e n = 1$ 时，不同的 $s_i$ 所代表区间一定不交，所以贡献为 $t$ 。
随时间增加，不同段所代表区间会相交，如下图， $s_1+len>=2$ ，此时 $1$ 的贡献为常数 $s_2-s_1$ ，总贡献为 $t-1)(len)+s_2-s_1$ 。

依次类推，贡献可以表达成 $f (x) = k x + b$ 的一次函数，且关于 $x$ 分成了 $t$ 段。

将 $s_i$ 升序排序，处理出 $d_i=s_{i+1}-s_i$ ，再将 $d_i$ 升序排序，得到两两合并的时间，可以预处理出这个分段函数。询问是二分值域即可。

在这里插入图片描述

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first 
#define sc second
using namespace std;
typedef double db;
typedef long long ll;
const int N=2e5+10;

int n,m,cnt,res;
ll a[N],cs;

struct P{
   ll k,b,t;
}t[N];

struct pr{
	ll d;int id;
	bool operator<(const pr&ky)const{
	   return d<ky.d;
	}
}q[N];

inline int fd(ll x)
{
	int re,l=1,r=cnt,mid;
	for(;l<=r;){
		mid=(l+r)>>1;
		if(t[mid].t<=x) l=(re=mid)+1;
		else r=mid-1;
	}
	return re;
}

int main(){
	int i,j;ll l,r,len;
    scanf("%d",&n);
    for(i=1;i<=n;++i) scanf("%I64d",&a[i]);
    sort(a+1,a+n+1);n=unique(a+1,a+n+1)-a-1; 
	if(n==1){
    
    for(scanf("%d",&m);m;--m){
    	scanf("%I64d%I64d",&l,&r);
    	printf("%I64d ",r-l+1);
	}
	return 0;	
	}
   
    for(i=1;i<n;++i) q[i].d=a[i+1]-a[i],q[i].id=i;
    sort(q+1,q+n);
    t[(cnt=1)].k=n;t[1].b=0;t[1].t=0;
    res=n;
    for(i=j=1;i<n;i=j){
    	for(;j<n && q[j].d==q[i].d;++j){
    	    res--;cs+=q[i].d;
		}
		t[++cnt].k=res;t[cnt].b=cs;t[cnt].t=q[i].d;
	}
    for(scanf("%d",&m);m;--m){
    	scanf("%I64d%I64d",&l,&r);
    	len=r-l+1;
        j=fd(len);
        printf("%I64d\n",(t[j].k*len+t[j].b));
	}
	return 0;
}

*E.Pavel and Triangles

每个三元组形如： $(2^j,2^i,2^i)(0\leq j\leq i)$

考虑 $a_j$ 中的 $2$ 个数能匹配上的 $1$ 个数的选择范围是 $a_i(i>j)$ 的子集 $[a_1,a_j]\subset[a_1,a_i]$ 。所以贪心让更前面的数匹配更多。

逐个考虑 $a_1-a_n$ ，每次用尽量多的2取匹配剩下的1，假设还剩下 $b_i$ 个数，那么再匹配 $\lfloor\frac{b_i}{3}\rfloor$ ，最后把 $b_i\%3$ 加入剩下的1中即可。

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first 
#define sc second
using namespace std;
typedef double db;
typedef long long ll;
const int N=3e5+10;

int n,a[N];
ll ans;int nw,res;

int main(){
	int i,j;ll l,r,len;
    scanf("%d",&n);
    for(i=1;i<=n;++i) {
    	scanf("%d",&a[i]);
    	j=min(a[i]>>1,nw);
    	a[i]-=(j<<1);nw-=j;ans+=j;
    	ans+=(a[i]/3);a[i]%=3;
    	nw+=a[i]; 
	}
	printf("%I64d",ans);
	return 0;
}

*F.Niyaz and Small Degrees

当 $x$ 是定值时，可以 $O(n\log n)$ DP求出答案：
设 $f_{i,j=(1/0)}$ 表示 $i$ 子树中所有点度数 $\leq x$ ，且不选/选 $i$ 的父边(选则 $i$ 的度数 $< x$ )时的最优答案。
转移时设 $tot=\sum \limits_{k\in son_i}f_{k,0}$ ，将点按 $f_{k,1}-f_{k,0}$ 排序， $f [i] [j]$ =前 $degree_i-x-j$ 小的 $f_{k,1}-f_{k,0}$ + $t o t$

发现此时 $degree_i\leq x$ 的点 $i$ 都是无用的， $(i,j)(degree_j\leq x)$ 的边一定会选，而 $i,j)(degree_i>x)$ 的边的选择仅仅由 $j$ 决定，所以只需要对所有 $d e g r e e > x$ 的点连通块进行DP。

而 $\sum\limits_{i=0}^{n-1}\sum\limits_{j=1}^n[degree_j>i]=\sum\limits_{i=1}^ndegree_i=2(n-1)$
总复杂度依旧是 $O(n\log n)$ 的

对于每个点 $i(degree_i>x)$ ，用堆维护所有它指向 $degree\leq x$ 的点的边的边权 $w$ ，和所有它指向的 $d e g r e e > x$ 的点 $j$ 的 $f_{j,1}-f_{j,0}$ ， $tot=\sum \limits_{j\in son_i,degree_j>x}f_{j,0}$ ，取堆中前 $degree_i-x$ 小的加入答案即可。

p.s

注意每次要动态插入 $f_{j,1}-f_{j,0}$ ，所以堆不够用了，可以在treap上二分。

若堆中第 $k>degree_i-i-j$ 小的值 $< 0$ ，则贪心把堆中值往 $t o t$ 加，直到堆空或者值 $\geq 0$ 。

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first
#define sc second
using namespace std;
typedef double db;
typedef long long ll;
const int N=3e5+10;
const ll inf=1e18;

int n,d[N],id[N],lim;
ll ans[N],f[N][2],nw;
set<pii>e[N];
vector<int>hv[N];
int vs[N],tim;

namespace treap{
#define lc(x) t[x].ch[0]
#define rc(x) t[x].ch[1]
#define F(x) t[x].fa
struct node{
	int sz,ch[2],fa,rnd;ll ss,v;
	inline void ini(ll x){v=ss=x;rnd=rand();sz=1;}
}t[N<<2];
int cnt,rt[N];

inline void pu(int x)
{
	if(!x) return;
	t[x].ss=t[lc(x)].ss+t[x].v+t[rc(x)].ss;
	t[x].sz=t[lc(x)].sz+1+t[rc(x)].sz;
}

inline void rot(int x)
{
	int y=F(x),z=F(y),dr=(rc(y)==x);
	t[y].ch[dr]=t[x].ch[dr^1];
	if(t[y].ch[dr]) F(t[y].ch[dr])=y;
	F(x)=z;if(z) t[z].ch[(rc(z)==y)]=x;
	t[x].ch[dr^1]=y;F(y)=x;pu(y);pu(x);
}

int ins(int x,ll v)
{
	if(!x) {t[++cnt].ini(v);return cnt;}
	int dr=v>t[x].v,y;t[x].sz++;t[x].ss+=v;
	t[x].ch[dr]=y=ins(t[x].ch[dr],v);F(y)=x;
	if(t[y].rnd>t[x].rnd) rot(y),x=y;
	return x;
}

int merge(int x,int y)
{
	F(x)=F(y)=0;//!!!
	if((!x) || (!y)) return (x|y);
	if(t[x].rnd>t[y].rnd){
		rc(x)=merge(rc(x),y);F(rc(x))=x;
		pu(x);return x;
	}else{
		lc(y)=merge(x,lc(y));F(lc(y))=y;
		pu(y);return y;
	}
} 

int del(int x,ll v)
{
	if(t[x].v==v) return merge(lc(x),rc(x));
	int dr=v>t[x].v,y;t[x].sz--;t[x].ss-=v;
	t[x].ch[dr]=y=del(t[x].ch[dr],v);F(y)=x;
	return x;
} 

inline int al(int x)
{
	int re=0;
	for(;x;){
		if(t[x].v<0) {re+=t[lc(x)].sz+1;x=rc(x);}
		else x=lc(x);
	}
	return re;
}

inline ll ask(int x,int kth)
{
	ll re=0,vl;int res;
	for(;kth;){
		res=t[lc(x)].sz+1;vl=t[x].ss-t[rc(x)].ss;
		if(kth<res) x=lc(x);
		else kth-=res,x=rc(x),re+=vl;
	}
	return re;
}


}
using namespace treap;

//f[x][0]不一定 <= f[x][1] 
void dfs(int x,int fr)
{
	vs[x]=tim;ll sum=0;
	for(pii y:e[x]) if(y.fi!=fr){
	   dfs(y.fi,x);sum+=f[y.fi][0];
	   rt[x]=ins(rt[x],f[y.fi][1]-f[y.fi][0]+y.sc);
	}
	int pos=al(rt[x]);
	f[x][0]=ask(rt[x],max(pos,d[x]-lim))+sum;f[x][1]=ask(rt[x],max(pos,d[x]-lim-1))+sum;
	for(pii y:e[x]) if(y.fi!=fr)
	  rt[x]=del(rt[x],f[y.fi][1]-f[y.fi][0]+y.sc);
}

inline bool cmp(int a,int b){return d[a]>d[b];}

int main(){
	int i,j,x,y,z;ll sum=0;
    scanf("%d",&n);
    for(i=1;i<n;++i){
    	scanf("%d%d%d",&x,&y,&z);sum+=z;d[x]++;d[y]++;
    	e[x].insert(mkp(y,z));e[y].insert(mkp(x,z));
	}ans[0]=sum;
	for(i=1;i<=n;++i) id[i]=i,hv[d[i]].pb(i);
	sort(id+1,id+n+1,cmp);
	for(lim=1;lim<n;++lim){
		tim++;nw=0;
		for(int x:hv[lim])
		 for(pii y:e[x]) if(d[y.fi]>lim)
		  rt[y.fi]=ins(rt[y.fi],y.sc),e[y.fi].erase(e[y.fi].lower_bound(mkp(x,0)));
		for(i=1;i<=n && d[id[i]]>lim;++i) if(vs[id[i]]!=tim){
			dfs(id[i],0);nw+=f[id[i]][0]; 
		}
		ans[lim]=nw;
	}
	for(i=0;i<n;++i) printf("%I64d ",ans[i]);
	return 0;
}

*G.Get Ready for the Battle

巧妙构造

设 $k_i=(\sum \limits_{j=1}^ihp_j)\%n(1\leq i<m)$ ， $k_0=0,k_m=n$ ，并将 $k_i$ 升序排序。

构造 $m$ 个组人数分别为 $sz_i=k_{i+1}-k_i(0\leq i<m)$ ，暴力模拟即可。

这样一定可以达到下界 $\lceil\frac{\sum\limits_{i=1}^mhp_i}{n}\rceil$

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first
#define sc second
using namespace std;
typedef double db;
typedef long long ll;
const int N=1e6+10;

int n,m,t[N],ss,a[N],cnt;
int v[N],sz[N],ans[N<<2],tot;

int main(){
	int i,j,x,y,z;
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;++i) {
    	scanf("%d",&t[i]);
    	ss+=t[i];v[i]=ss%n;
	}v[m]=n;
	sort(v,v+m+1);
	for(i=1;i<=m;++i) sz[i]=v[i]-v[i-1];

	for(i=1,j=1;i<=m;++i){
		x=t[i];
		for(;x>0;){
			ans[++tot]=i;
			x-=sz[j];
			j=(j==m)?1:(j+1); 
		}
	}
	for(;tot%m;) ans[++tot]=1;
	printf("%d\n",tot/m);	
	for(i=1;i<=m;++i) printf("%d ",sz[i]);
	puts("");
	for(i=1;i<=tot;++i){
		printf("%d ",ans[i]);	
		if(!(i%m)) puts("");
	}
	return 0;
}

H.Triple

暴力FWT：

对于第 $i$ 个数列，设为多项式 $f_i$ ，其中 $f_{i,a_i}=x,f_{i,b_i}=y,f_{i,c_i}=z$ ，其余项为0， $f_i$ 分别FWT后将 $\prod fwt(f_i)$ 逆变换就得到了答案。复杂度 $O(n2^kk)$

优化暴力：

将序列整体异或上 $a_i$ ，即 $f_{i,0}=x,f_{i,a_i\ xor \ b_i}=y,f_{i,a_i\ xor \ c_i}=z$ 。

发现FWT后每一项的系数只可能是以下4种之一：
$x + y + z, x + y - z, x - y + z, x - y - z$

单独处理每一位 $i(0\leq i<2^k)$ ，假设当前位出现了 $a$ 个 $x + y + z$ , $b$ 个 $x + y - z$ , $c$ 个 $x - y + z$ , $d$ 个 $x - y - z$ ，则乘积为 $x+y+z)^a(x+y-z)^b(x-y+z)^c(x-y-z)^d$

考虑将 $f_{a_i\ xor\ b_i}+=1(1\leq i\leq n)$ 和 $f_{a_i\ xor \ c_i}+=1$ 分别做FWT， $cnt_b$ 表示前者当前位的系数， $cnt_c$ 表示后者当前位的系数。

则
$a + b + c + d = n$
$a+b-c-d=cnt_b$
$a-b+c-d=cnt_c$

再需要一个等式就可以求出 $x, y, z, w$ ，所以加一个 $f_{b_i\ xor \ c_i}+=1$ ，即 $a - b - c + d = c n t$ (考虑分别将 $f_{b_i},f_{c_i}$ FWT后合并， $1)^2=(-1)^2=1,1(-1)=-1$ )

复杂度 $O(2^kk)$

#include<bits/stdc++.h>
#define gc getchar
#define pb push_back
#define mkp make_pair
#define pii pair<int,int>
#define fi first
#define sc second
const int mod=998244353;
using namespace std;
typedef double db;
typedef long long ll;
const int N=(1<<17)+100;

int n,m,a,b,c,S,nv2,nv4;
int f[N],g[N],h[N];

inline void ad(int &x,int y){x+=y;if(x>=mod) x-=mod;}
inline void dc(int &x,int y){x-=y;if(x<0) x+=mod;}
inline int inc(int x,int y){x+=y;return x>=mod?x-mod:x;}
inline int dec(int x,int y){x-=y;return x<0?x+mod:x;}

inline int fp(int x,int y)
{
	int re=1;
	for(;y;y>>=1,x=(ll)x*x%mod)
	 if(y&1) re=(ll)re*x%mod;
	return re; 
}

inline void fwt(int *e)
{
	int i,j,k,x,y;
	for(i=1;i<S;i<<=1){
		for(j=0;j<S;j+=(i<<1))
		 for(k=0;k<i;++k){
		 	x=e[j+k];y=e[i+j+k];
		 	e[j+k]=inc(x,y);e[i+j+k]=dec(x,y);
		 }
	}
}

inline void ifwt(int *e)
{
	int i,j,k,x,y;
	for(i=1;i<S;i<<=1){
		for(j=0;j<S;j+=(i<<1))
		 for(k=0;k<i;++k){
		 	x=e[j+k];y=e[i+j+k];
		 	e[j+k]=(ll)inc(x,y)*nv2%mod;
		    e[i+j+k]=(ll)dec(x,y)*nv2%mod;
		 }
	}
}

int main(){
	int i,j,x,y,z,w,A,B,C,D,tot=0;
	nv2=fp(2,mod-2);nv4=(ll)nv2*nv2%mod;
    scanf("%d%d%d%d%d",&n,&m,&a,&b,&c);
    S=1<<m;
    for(i=1;i<=n;++i){
	    scanf("%d%d%d",&x,&y,&z);
	    tot^=x;f[y^x]++;g[z^x]++;h[y^z]++; 
	}
	A=inc(a,inc(b,c));B=dec(inc(a,b),c);
	C=inc(dec(a,b),c);D=dec(a,inc(b,c));
	fwt(f);fwt(g);fwt(h);
	for(i=0;i<S;++i){
		x=(ll)inc(n,inc(f[i],inc(g[i],h[i])))*nv4%mod;
		y=dec((ll)inc(n,f[i])*nv2%mod,x);
		z=dec((ll)inc(n,g[i])*nv2%mod,x);
		w=dec((ll)inc(n,h[i])*nv2%mod,x);
		f[i]=fp(A,x)*(ll)fp(B,y)%mod*(ll)fp(C,z)%mod*(ll)fp(D,w)%mod; 
	}
	ifwt(f);
	for(i=0;i<S;++i) printf("%d ",f[i^tot]);
	return 0;
}