20171024测试 爆炸 最低成绩


得分:40 + 0 + 0

第一题 tarjan求双连通分量 + 错误结论 哈 哈 哈
二三题真是暴力都舍不得打啊

T1: Graph
题意:给你一个图 求最少添多少条边 能使任意两点间有两条通路
我自己想的是 先在非双连通分量中 统计叶子节点 然后统计有多少颗子树和入度为1的双连通分量
其实真的想多了 
标算:tarjan缩点后 把双连通分量看作一个点 再整个图连边 统计叶子节点的数量
其实标算也麻烦了 我发现tarjan后 直接统计入度为1的分量 整个双连通分量或一个点 都看做一个分量
因为连边后的节点就是一个分量 连边后叶子节点就是入度为1的分量
直接统计有多少个分量 ans=(int)(num+1)/2;
#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,m,top,tot,tim,anum,s;
int low[100005],dfn[100005],first[100005],q[100005],num[100005],belong[100005];
int x[200005],y[200005],out[100005];
bool ex[100005];
struct node{
	int next,to;
}e[400005];
int read(){
	int i=0,f=1;char ch=getchar();
	for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') f=-1;
	for(;ch<='9'&&ch>='0';ch=getchar()) i=(i<<3)+(i<<1)+ch-'0';
	return i*f;
}
void w(int x){
	if(x>=10) w(x/10);putchar(x%10+'0');
}
void inser(int x,int y){
	tot++;
	e[tot].next=first[x];
	first[x]=tot;
	e[tot].to=y;
	tot++;
	e[tot].next=first[y];
	first[y]=tot;
	e[tot].to=x;
}
void dfs(int x,int fa){
	dfn[x]=low[x]=++tim;
	q[++top]=x;ex[x]=true;
	for(int u=first[x];u;u=e[u].next){
		int to=e[u].to;
		if(to!=fa){
			if(!dfn[to]) dfs(to,x),low[x]=min(low[x],low[to]);
			else if(ex[to]) low[x]=min(low[x],dfn[to]);
		}
	}
	if(dfn[x]==low[x]){
		int v=-1;tot++;
		while(x!=v){
			v=q[top--];
			num[tot]++;
			belong[v]=tot;
			ex[v]=false;
		}
	}
}
int main(){
	n=read();m=read();tot=1;
	for(int i=1;i<=m;++i){
		x[i]=read();y[i]=read();
		inser(x[i],y[i]);
	}	
	tot=0;
	dfs(1,0);
	for(int i=1;i<=m;++i){
		if(belong[x[i]]!=belong[y[i]]) out[belong[x[i]]]++,out[belong[y[i]]]++;
	}
	for(int i=1;i<=tot;++i)
		if(out[i]==1) anum++;
	w((anum+1)/2);
}

T2: Mul

题意:定义无平方因子数:不能被任意一个质数的平方整除
T组数据 每组数据给出一个n和k
求出 选择不超过k个n以内的正整数相乘 使得乘积为一个无平方因子数
求有多少种取法 答案对1e9+7取模
这道题暴力可以打70分
来说下标算 
状压DP
每个数只会含有一个比 sqrt(N) 大的质因数,我们把所有数字按含有比sqrt(N)大的分组
sart(N)内的质数只有8个
dp[i][j] 表示选了 i 个数 状态为 j 的方案数
#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int p[8]={2,3,5,7,11,13,17,19};
int T,n,k,a[505],sta[505],dp[505][260],temp[505];
vector<int>v[505];
inline void inc(register int &a,register int &b){
	if(a+b>=mod) a=a+b-mod;
	else a=a+b;
}
int main(){
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&n,&k);
		for(int i=1;i<=n;++i) v[i].clear(),sta[i]=0;
		for(int i=1;i<=n;++i){	//先处理1
			a[i]=i;
			for(int j=0;j<8;++j)
				if(a[i]%(p[j]*p[j])==0) sta[i]=-1;	// -1的补码为1111 1111
				/*在计算机系统中,数值一律用补码来表示和存储。
				原因在于,使用补码,可以将符号位和数值域统一处理;
				同时,加法和减法也可以统一处理。
				补码与原码相互转换,其运算过程是相同的,不需要额外的硬件电路。*/ 
				else if(a[i]%p[j]==0) a[i]/=p[j],sta[i]|=(1<<j);	//状态:是否含有 p[j]因子
		}
		for(int i=1;i<=n;++i){
			if(~sta[i]){		// ~表示按位取反 -1取反返回0
				if(a[i]==1) v[i].push_back(i);
				else v[a[i]].push_back(i);
			}
		}
		//每个数只会含有一个比sqrt(N)大的质因数,我们把所有数字按含有比sqrt(N)大的数分组
		memset(dp,0,sizeof(dp));
		dp[0][0]=1;
		for(int i=1;i<=n;++i)
			if(~sta[i]&&v[i].size()){
				int sz=v[i].size();
				for(int j=0;j<sz;++j) temp[j]=sta[v[i][j]];
				for(int j=k-1;j>=0;--j)			//转移 
					for(int l=0;l<sz;++l)
						for(int v=temp[l],s=255^v,x=255^v;;s=(s-1)&x){
						//255^v和~v都可以 
						/*一共8个质因数 1111 1111为255
						s为所有不含有 v状态的情况
						证:255^v 1111 1111^v=x v为0的位上x为1 v为1的位上x为0
						所以说255^v和~v都可以 因为 ~ 本来就是对八位取反
						那么 s&x  v为1的位上 x为0 s为0 所以s一定不含状态v*/ 
						if((s&temp[l])==0) inc(dp[j+1][s|v],dp[j][s]);	//加
						if(!s) break;
					}
			}
		long long ans=0;
		for(int i=1;i<=k;++i)
			for(int j=0;j<256;++j)
				ans+=dp[i][j];
		cout<<ans%mod<<endl;
	}
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值