24.2.25周报(匈牙利 , 线段树)

星期一:

话接上回cf round927 div3, 很好的开局,但D题小模拟写瓦了一个很蠢的点,de了半小时没de出来,只能说还是粗糙了

cf round928 div4,打的一坨,止步于C

题意:自行理解

思路:一开始没注意时限只有半秒,写了个暴力,tle,简单优化下,又tle,开始疯狂数学推导,但推了一个多小时到最后发现好像推不出来, 真没想到是前缀和, 我的

因为答案有重复性(!!!!),加上1e4个询问很多,所以可以边跑边把前面的答案用前缀和存下来 ,算是被这题狠狠拿下了,好狠毒啊

星期三:

打了场牛客,学了个匈牙利算法,放个板子:

const int N=2e6+10;
const int mod=998244353;
ll n;
int m;
struct nod{
	int to,nex;
}e[N];
int h[N],cnt,ma[N];
bool vi[N];
bool dfs(int u){
	for(int i=h[u];i;i=e[i].nex){
		int v=e[i].to;
		if(vi[v]) continue;
		vi[v]=1;
		if(!ma[v] || dfs(ma[v])){
			ma[v]=u;
			return 1;
		}
	}
	return 0;
}
void solve(){
	cin >> n >> m;
	for(int i=1;i<=m;i++){
		int u,v;
		cin >> u >> v;
		add(u,v);              // 单向边
	}
	ll ans=0;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++) vi[j]=0;
		if(dfs(i)) ans++;
	}
	cout << ans;
}

挂道题:牛客寒假竞赛5 M

思路:想了下,觉得可能是dp,后面又感觉是并查集,再后面实在想不到,看题解居然是分类讨论

星期四:

学了下线段树,操作有构造线段树 O(n),单点修改,区间修改,区间查询均为O(logn)

板子如下:

#define lc p<<1
#define rc p<<1|1
const int N=2e6+10;
const int mod=998244353;
ll n;
ll w[N];
struct nod{
	ll l,r,sum,add;
}t[N*4];
void pushup(int p){
	t[p].sum=t[lc].sum+t[rc].sum;
}
void pushdn(int p){
	if(t[p].add){
		t[lc].sum+=(t[lc].r-t[lc].l+1)*t[p].add;
		t[rc].sum+=(t[rc].r-t[rc].l+1)*t[p].add;
		t[lc].add+=t[p].add;
		t[rc].add+=t[p].add;
		t[p].add=0;
	}
}
void bd(int p,int l,int r){
	t[p]={l,r,w[l],0};
	if(l==r) return ;
	int m=l+r>>1;
	bd(lc,l,m);
	bd(rc,m+1,r);
	pushup(p);
}
void updt(int p,int x,int y,int k){
	if(t[p].l>=x && t[p].r<=y){
		t[p].sum+=(t[p].r-t[p].l+1)*k;
		t[p].add+=k;
		return ;
	}
	int m=t[p].l+t[p].r>>1;
	pushdn(p);
	if(x<=m) updt(lc,x,y,k);
	if(y>m) updt(rc,x,y,k);
	pushup(p);
}
ll que(int p,int x,int y){
	if(t[p].l>=x && t[p].r<=y) return t[p].sum;
	int m=t[p].l+t[p].r>>1;
	pushdn(p);
	ll sum=0;
	if(x<=m) sum+=que(lc,x,y);
	if(y>m) sum+=que(rc,x,y);
	return sum;
}

星期五:

贴道断断续续做了一上午的模拟:挪象棋

  写了个模拟一直wa,找不到原因,写个暴力dij又t了, 然后嗯造了10个样例, 让模拟和dij各跑一遍,找出了错误样例,算是学到了个给模拟debug的法子吧

下午:牛客寒假训练营6

晚上cf, C题没敲出来

周日:

补道题

思路: 先打表 (main函数里打), 然后有同一素数因子的数放到一个并查集,最后若并查集只有一个即无解

不过有一点让我很困惑, iff数组的清零如果放solve最后就会wa

代码如下:

const int N=1e6+10;
const int mod=998244353;
ll n;
int a[N];
vector<int>ve,p[N];
int iff[N],fa[N];
int fnd(int x){
	return fa[x]==x?x:fa[x]=fnd(fa[x]);
}
void solve(){
	cin >> n;
	ve.clear();
	for(int i=1;i<=n;i++) fa[i]=i;
	for(int i=1;i<=n;i++){
		cin >> a[i];
		for(auto j:p[a[i]]){
			if(!iff[j]){
				iff[j]=i;
				ve.push_back(j);
			}else{
				int x=fnd(i),y=fnd(iff[j]);
				fa[x]=y;
			}
		}
	}
	for(auto i:ve) iff[i]=0;
	int num=0;
	for(int i=1;i<=n;i++)
		if(fnd(i)!=fnd(1)) num++;
	if(!num){cout << "-1 -1\n"; return ;}
	cout << num << " " << n-num << "\n";
	for(int i=1;i<=n;i++)
		if(fnd(i)!=fnd(1)) cout << a[i] << " ";
	cout << "\n";
	for(int i=1;i<=n;i++)
		if(fnd(i)==fnd(1)) cout << a[i] << " ";
	cout << "\n";
}
int main(){
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	p[1].push_back(1);
	for(int i=2;i<N;i++){
		if(p[i].empty()){
			for(int j=i;j<N;j+=i) p[j].push_back(i);
		}
	}
	int _=1;
	cin >> _;
	while(_--) solve();
	return 0;
}