csp-s模拟测试54x,y,z题解

题面:https://www.cnblogs.com/Juve/articles/11606834.html

x:

并差集,把不能分到两个集合里的元素和并到一起,设连通块个数为cnt,则答案为:$2^cnt-2$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#include<vector>
#define re register
#define int long long
using namespace std;
const int MAXN=1e5+5;
const int mod=1e9+7;
int n,a[MAXN],T,ans,num;
inline int gcd(re int a,re int b){
	return b==0?a:gcd(b,a%b);
}
inline int q_pow(re int a,re int b,re int p){
	re int res=1;
	while(b){
		if(b&1) res=res*a%p;
		a=a*a%p;
		b>>=1;
	}
	return res;
}
int fa[MAXN];
inline int find(re int x){
	return fa[x]=(fa[x]==x?x:find(fa[x]));
}
vector<int>v[MAXN*10];
int prime[MAXN*10],vis[MAXN*10],tot=0,pri[MAXN*10];
void get_prime(int N){
    vis[1]=1;
    for(int i=2;i<=N;i++){
        if(!vis[i]) prime[++tot]=i,pri[i]=tot;
        for(int j=1;j<=tot&&i*prime[j]<=N;j++){
            vis[i*prime[j]]=1;
            if(!(i%prime[j])) break;
        }
    }
}
inline void divi(re int x,re int pos){
	for(re int i=1;i<=tot&&prime[i]*prime[i]<=x;++i){
		if(x%prime[i]==0){
			v[i].push_back(pos);
			while(x%prime[i]==0) x/=prime[i];
		}
	}
	if(x>1) v[pri[x]].push_back(pos);
}
int maxx=0;
signed main(){
	get_prime(1e6);
	scanf("%lld",&T);
	while(T--){
		scanf("%lld",&n);
		ans=num=maxx=0;
		for(re int i=1;i<=n;++i){
			fa[i]=i;
			scanf("%lld",&a[i]);
			divi(a[i],i);
			maxx=max(maxx,a[i]);
		}
		for(int i=1;i<=tot&&prime[i]<=maxx;++i){
			int N=v[i].size();
			if(N==0) continue;
			int p=find(v[i][0]);
			for(int j=1;j<N;++j){
				int q=find(v[i][j]);
				if(p!=q) fa[q]=p;
			}
			v[i].clear();
		}
		for(int i=1;i<=n;++i){
			if(fa[i]==i) ++num;
		}
		printf("%lld\n",((q_pow(2,num,mod)-2)%mod+mod)%mod);
	}
	return 0;
}

y:

定义dp:f[i][j][k]表示走了i步,当前点是j,状态为k的一个bool数组,表示是否有该状态存在

为节省空间我们分两部分转移,这样第一维只有$\frac{d}{2}$,第三维只有$2^{\frac{d}{2}}$,

初始状态:$f1[0][1][0]=1,f2[0][i][0]=1(i \in n)$,

我们用连通性转移,最后统计答案是把两个dp数组的状态和在一起

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int MAXN=95;
int n,m,d,len1,len2,ans=0;
bool cal1[13][MAXN][(1<<13)+5],cal2[13][MAXN][(1<<13)+5];
vector< pair<int,int> >mp[MAXN];
signed main(){
	scanf("%d%d%d",&n,&m,&d);
	len1=d/2,len2=d-len1;
	cal1[0][1][0]=1;
	for(int i=1;i<=n;++i) cal2[0][i][0]=1;
	for(int i=1,u,v,c;i<=m;++i){
		scanf("%d%d%d",&u,&v,&c);
		mp[u].push_back(make_pair(v,c));
		mp[v].push_back(make_pair(u,c));
	}
	for(int i=0;i<len1;++i){
		int s=(1<<i);
		for(int j=0;j<s;++j)
			for(int k=1;k<=n;++k){
				int N=mp[k].size();
				for(int p=0;p<N;++p)
					cal1[i+1][mp[k][p].first][(j<<1)+mp[k][p].second]|=cal1[i][k][j];
			}
	}
	for(int i=0;i<len2;++i){
		int s=(1<<i);
		for(int j=0;j<s;++j)
			for(int k=1;k<=n;++k){
				int N=mp[k].size();
				for(int p=0;p<N;++p)
					cal2[i+1][mp[k][p].first][(j<<1)+mp[k][p].second]|=cal2[i][k][j];
			}
	}
	for(int i=0;i<(1<<d);++i)
		for(int j=1;j<=n;++j)
			if(cal1[len1][j][i>>len2]&&cal2[len2][j][i&((1<<len2)-1)]){
				++ans;
				break;
			}
	printf("%d\n",ans);
	return 0;
}

z:

全机房没几个作出来的,只能放标算了,去题面里找

原文地址:https://www.cnblogs.com/Juve/p/11606844.html

时间: 2024-07-31 19:11:02

csp-s模拟测试54x,y,z题解的相关文章

[CSP-S模拟测试]:小Y的图(最小生成树+LCA)

题目传送门(内部题131) 输入格式 第一行三个整数$n$.$m$和$Q$. 接下来$m$行每行三个整数$x$.$y$.$z$($1\leqslant x,y\leqslant n,1\leqslant z\leqslant 1,000,000$),表示有一条连接$x$和$y$长度为$z$的边. 接下来$Q$行每行两个整数$x$.$y$($x\neq y$),表示一组询问. 输出格式 $Q$行每行一个整数,表示一组询问的答案. 样例 样例输入: 5 5 41 2 31 3 23 2 11 4 5

2018-10-25 模拟测试题解

目录 问题 A: 魏传之长坂逆袭 题目描述 输入 输出 样例输入 样例输出 题解 问题 B: 蜀传之单刀赴会 题目描述 [问题描述] 输入 输出 样例输入 样例输出 题解 问题 C: 吴传之火烧连营 [题目背景] [问题描述] 输入 输出 样例输入 样例输出 [样例解释] [数据规模和约定] 题解 本篇题解也发表于zwcblog作者是同一个人 问题 A: 魏传之长坂逆袭 题目描述 众所周知,刘备在长坂坡上与他的一众将领各种开挂,硬生生从曹操手中逃了出去,随后与孙权一起火烧赤壁.占有荆益.成就霸业

[CSP-S模拟测试53]题解

A.u 只涉及到区间修改可以考虑差分,然而如果每一行都差分复杂度还是过高.我们发现差分标记也是连续的(一行横着的一行斜着的),所以可以维护两个 差分的差分,扫两遍统计即可. #include<cstdio> #include<iostream> #include<cstring> using namespace std; typedef long long ll; const int N=2005; int read() { int x=0,f=1;char ch=ge

2018冬令营模拟测试赛(三)

2018冬令营模拟测试赛(三) [Problem A]摧毁图状树 试题描述 输入 见"试题描述" 输出 见"试题描述" 输入示例 见"试题描述" 输出示例 见"试题描述" 数据规模及约定 见"试题描述" 题解 这题没想到贪心 QwQ,那就没戏了-- 贪心就是每次选择一个最深的且没有被覆盖的点向上覆盖 \(k\) 层,因为这个"最深的没有被覆盖的点"不可能再有其它点引出的链覆盖它了,而它又

2018冬令营模拟测试赛(五)

2018冬令营模拟测试赛(五) [Problem A][UOJ#154]列队 试题描述 picks 博士通过实验成功地得到了排列 \(A\),并根据这个回到了正确的过去.他在金星凌日之前顺利地与丘比签订了契约,成为了一名马猴烧酒. picks 博士可以使用魔法召唤很多很多的猴子与他一起战斗,但是当猴子的数目 \(n\) 太大的时候,训练猴子就变成了一个繁重的任务. 历经千辛万苦,猴子们终于学会了按照顺序排成一排.为了进一步训练,picks 博士打算设定一系列的指令,每一条指令 \(i\) 的效果

2018冬令营模拟测试赛(十七)

2018冬令营模拟测试赛(十七) [Problem A]Tree 试题描述 输入 见"试题描述" 输出 见"试题描述" 输入示例 见"试题描述" 输出示例 见"试题描述" 数据规模及约定 见"试题描述" 题解 这个数据范围肯定是树上背包了. 令 \(f(i, j, k)\) 表示子树 \(i\) 中选择了 \(j\) 个节点,路径与根的连接情况为 \(k\),具体地: \(k = 0\) 时,路径的两个端点

2018冬令营模拟测试赛(十九)

2018冬令营模拟测试赛(十九) [Problem A]小Y 试题描述 输入 见"试题描述" 输出 见"试题描述" 输入示例 见"试题描述" 输出示例 见"试题描述" 数据规模及约定 见"试题描述" 题解 目前未知. 这题目前就能做到 \(O(n \sqrt{M} \log n)\),其中 \(M\) 是逆序对数,然而会被卡 \(T\):当然这题暴力可以拿到和左边那个算法一样的分数,只要暴力加一个剪枝:当左

noip模拟测试11

T1:string 第一眼秒出思路,这不就是排序那道题的加强版吗? 然而歪?解复杂度虽然是对的,但常数过大,竟被卡到70 歪?解:(实际上std写的就是这个,但据说std被卡掉了 OAO) 因为字符集很小,所以我们可以把区间排序改为区间查询和覆盖 即:先查询区间内所有字符的个数,再从左端点开始按照大小关系依次将长度为字符个数的区间修改为该字符. 期望复杂度O ( 26*mlogn ),实际复杂度O ( 26*mlogn*(巨大的常数) ) 所以需要一(feng)定(kuang)的卡常 正?解:

2019.9.28 csp-s模拟测试54 反思总结

咕咕咕的冲动如此强烈x T1x: 看完题目想了想,感觉把gcd不为1的强行放在一组,看作一个连通块,最后考虑连通块之间的组合方式就可以了. 然后维护这个连通块可以写并查集可以连边跑dfs怎么着都行… 然而我在处理数字分解质因数这里T掉了,原因是一个很显然的优化写法我基本没怎么写过.线性筛的时候记录每个数是被哪个质数标记过的,分解一个数的时候直接处理记录下来的质数就可以. #include<iostream> #include<cstdio> #include<cmath>