CF-round#599-1242

最新推荐文章于 2022-07-13 02:01:45 发布

蠕虫君

最新推荐文章于 2022-07-13 02:01:45 发布

阅读量179

点赞数

分类专栏：学习笔记 ACM 文章标签：整场笔记 CF

本文链接：https://blog.csdn.net/qq_41046771/article/details/102959626

版权

ACM 同时被 2 个专栏收录

68 篇文章 0 订阅

订阅专栏

学习笔记

64 篇文章 0 订阅

订阅专栏

这篇博客主要分享了在Codeforces Round #599中的解题心得，涉及图论问题的BFS求解0-1 MST的特殊情况，以及高级DFS结合状压DP解决Sum Balance问题。在BFS部分，讨论了如何处理非完全图的连通分量；在DFS部分，探讨了寻找数字平衡方案的思路，包括特判、环的寻找和状压DP的应用。

摘要由CSDN通过智能技术生成

图论-BFS

题目:0-1 MST
链接:https://codeforces.com/contest/1242/problem/b
大意:给你一张完全图,n<=1e5个结点,其中最多有m<=1e5条边权为1的边，剩余的边边权全为0
分析:求无向补图最少构成几个连通分量:如果不是补图，那我们就直接BFS了，。。。。然而不是补图还是只能BFS,但BFS的方法要发生变化，原本用u拓展v的充要条件是:存在边u->v,现在的条件是:v还未被拓展&&u和v之间不存在权为1的边

#include <bits/stdc++.h>
using namespace std;
#define F(i, a, b) for (int i = int(a); i <= int(b); ++i)
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int maxn = 1e5+9;
set<int>bad[maxn],lef;
int n,m;
int comp;
int q[maxn],h;
int main()
{
	ios::sync_with_stdio(0);
	int T(1);
	while (T--)
	{
		cin>>n>>m;
		F(i,1,m)
		{
			int f,t;cin>>f>>t;
			bad[f].insert(t);bad[t].insert(f);
		}
		F(i,1,n)lef.insert(i);
		F(i,1,n)if(lef.count(i))//如果lef集合中有i这个元素
		{
			++comp;
			h=0;q[++h]=i;lef.erase(i);//lef中擦除i,i入队
			while(h)
			{
				int a=q[h--];
				vector<int>cool;
				for(auto&t:lef)if(!bad[a].count(t))cool.push_back(t);//新加入的点: lef& !bad[a]
				for(auto&t:cool){lef.erase(t);q[++h]=t;}
			}
		}
		cout<<comp-1<<'\n';
	}
}

//What to Debug
/*
-1.最好把全部warning都X掉,否则:https://vjudge.net/solution/19887176
0.看看自己是否有可能需要快读，禁endl
1.数组越界，爆int,浮点精度(查看精度是否达到题目要求,看有没有浮点数比较:eps)，取模操作，初始化数组，边缘数据,输出格式(cas),强制在线是否更新了las
2.通读代码，代码无逻辑错误
3.读题，找到题意理解失误或算法错误
4.放弃
*/

高级DFS-状压DP

题目:Sum Balance
链接:https://codeforces.com/contest/1242/problem/C
大意:告诉你有最多15个口袋，每个口袋里面最多有5000个数字，每个数字最大为1e9,且所有数字各不相同,你要从每个口袋里面拿出一个数字，然后把这些数字分别放回这些口袋(数字可以先随意打乱顺序再放回口袋),问你要一个方案，使得放回去之后每个口袋里的数字之和相等，或者说没有方案
分析:我太菜了,DIV1的C题做不出来就别想黄了。
- 首先肯定是要特判掉，如果sum%k!=0，直接返回没有方案
- 每一袋数字要达到最终状态，就肯定要使得袋中的数字+=(ave-sum[i])。为了使得袋中数字变为ave，就肯定要失去一个数字，再拿回来一个数字，由于所有的数字各不相同，所以每个数字都相当于有一个[换回指针],用它去换那个数字就能达到ave。
- 那最终方案是怎样的?最终方案肯定是若干个环，这个环的特点是每个口袋都用到了一个元素，环通过[换回指针]首尾相连，可能不止一个环，也可能存在单元素的环
- 第四步:使用map<ll,pair<int,int>>计数我们找出所有的[换回指针]，但是要注意有几个坑:
  - 不允许[换回指针]指向同行（口袋）的其他元素
  - 允许[换回指针]指向元素自己
- 第五:要找出所有的环，我们可以通过DFS(或者tarjan-用tarjan其实没那么好，因为我们DFS的过程中还有各种自定义的操作)。找环的时候又有几个坑:
  - 我自己不是很擅长找环----好吧菜了没得洗
  - 每个点有三种情况:没有出边，出边指向自己，出边指向别人(如果前面没有处理好出边还有可能指向同行元素)
- 好的，我们终于找出了所有的环，接下来想到与我们有n(15)个元素，有一些子集(环),我们要用环来不重复地cover掉所有元素,状压dp就可以了，这里我又犯了一个错，现在默写一遍:for(auto sub=tot:sub;sub=(sub-1)&tot)

#include <bits/stdc++.h>
using namespace std;
#include<unordered_map>
#define PII pair<int,int>
#define F(i, a, b) for (int i = int(a); i <= int(b); ++i)
typedef long long ll;
typedef unsigned long long ull;
const int inf = 0x3f3f3f3f;
const int maxn = 20;
const int maxm = 5e3 + 9;
int num[maxn], n;
ll sum[maxn], a[maxn][maxm];
int myans[maxn][2];
bool deng(PII A, PII B) { return A.first == B.first&&A.second == B.second; }

unordered_map<ll, PII>ma;
PII go[maxn][maxm];

vector<PII>ve;
int vis[maxn][maxm], vid(0);

int f[2 << maxn];PII ans[2 << maxn];
PII pre[2 << maxn];
void dfs3(int x, int y)
{
	ve.push_back({ x,y });
	vis[x][y] = vid;
	if (!go[x][y].first) return;
	int p = go[x][y].first, q = go[x][y].second;
	if (vis[p][q] == vid)
	{
		int sta = 1<<(p-1);
		int use[maxn] = { 0 };
		use[p] = 1;
		while (!deng(ve.back() , {p, q}))
		{
			int x = ve.back().first;
			if (use[x]) return;
			sta += 1<<(x-1); use[x] = 1; ve.pop_back();
		}
		pre[sta] = { x,y };
		f[sta] = 1; ans[sta] = { p,q };
	}
	else if (!vis[p][q]) dfs3(p, q);
}
int dfs(int x, int y)
{
	if (!vis[x][y])
	{
		ve.push_back({ x,y });
		vis[x][y] = vid;
		if (deng(go[x][y], { 0,0 }))
		{
			return 0;
		}
		else 
			return dfs(go[x][y].first, go[x][y].second);
	}
	else
	{
		if (vis[x][y] != vid)return 0;
		else
		{
			int con(0);
			ve.push_back({ x,y });
			vector<int>use(maxn, 0);
			do
			{
				auto v = ve.back();
				if (use[v.first])return 0;
				use[v.first] = 1;
				con |= 1 << (v.first - 1);
				ve.pop_back();
			} while (ve.size()&&!deng(ve.back(), { x, y }));
			pre[con] = { x,y };
			return 1;
		}
	}

}

void dfs2(int x)
{
	if (pre[x].second)
	{
		auto beg = pre[x];
		auto to = beg;
		do
		{
			auto from= go[to.first][to.second];
			myans[from.first][0] = a[from.first][from.second];
			myans[from.first][1] = to.first;
			to = go[to.first][to.second];
		} while (!deng(to, beg));
	}
	else dfs2(pre[x].first), dfs2(x-pre[x].first);
}
bool solve()
{
	if (sum[0] % n)return 0;
	ll ave = sum[0] / n;
	//建图
	F(i, 1, n) F(j, 1, num[i])
	{
		ll wan = ave - sum[i] + a[i][j];
		if (ma.find(wan) == ma.end())continue;
		auto ret = ma[wan];
		if (ret.first == i && ret.second != j)continue;
		go[i][j] = ma[wan];
	}
	F(i, 1, n)F(j, 1, num[i])if (!vis[i][j])
	{
		++vid;
		ve.clear();
		dfs(i, j);
	}
	int lim = (1 << n) - 1;
	F(i, 1, lim)if (deng(pre[i], { 0,0, }))
	{
		for (int j = i; j; j = (j - 1)&i)if (!deng({ 0,0 }, pre[j]) && !deng({ 0,0 },pre[i-j]))
		{
			pre[i] = { j,0 };
			break;
		}
	}

	if (deng(pre[lim], { 0,0 }))return 0;
	dfs2(lim); return 1;
}
int main()
{
	ios::sync_with_stdio(0);
	int T(1);
	while (T--)
	{
		cin >> n;
		F(i, 1, n)
		{
			cin >> num[i];
			F(j, 1, num[i])
			{
				cin >> a[i][j];
				sum[i] += a[i][j];
				sum[0] += a[i][j];
				ma[a[i][j]] = { i,j };
			}
		}
		if (solve())
		{
			cout << "Yes" << '\n';
			F(i, 1, n)cout << myans[i][0] << ' ' << myans[i][1] << '\n';
		}
		else cout << "No" << '\n';
	}
}

//What to Debug
/*
-1.最好把全部warning都X掉,否则:https://vjudge.net/solution/19887176
0.看看自己是否有可能需要快读，禁endl
1.数组越界，爆int,浮点精度(查看精度是否达到题目要求,看有没有浮点数比较:eps)，取模操作，初始化数组，边缘数据,输出格式(cas),强制在线是否更新了las
2.通读代码，代码无逻辑错误
3.读题，找到题意理解失误或算法错误
4.放弃
*/