【比赛报告】2018.8.10集训 NOIP练习赛卷三十五

A.细胞分裂 费马小定理降幂+快速幂取模

题目链接
问题描述

小 A 养了一大坨细胞。
最初小 A 只有 1 个细胞。 每秒， 小 A 的每个细胞都会分裂成 2 个细胞。
已知： 现在离“最初”已经过去了 x 秒， 那么现在的细胞数当然是可以计算的。
小 A 想知道的当然不是当前的细胞数。 小 A 知道他养的细胞的习性： 每 y 个细胞
会聚成一团。 经常会有剩下的细胞， 那么我们称这些细胞是孤独的。
小 A 想知道的就是孤独的细胞个数

输入

输入文件共一行， 为两个整数 x y， 以空格隔开

输出

输出文件共一行， 为一个整数， 即孤独的细胞个数。

输入样例

3 3

输出样例

2

数据规模和约定

对于 10%的数据， x<2^6。
对于 20%的数据， x<2^17。
对于 40%的数据， x<2^64。
对于 70%的数据， x<2^2333。
对于 100%的数据， 0≤x<2^233333， y 是 3 到 1000 之间（含两端） 的质数

---

题解

比赛时看到这题毫无头绪（我太弱了），首先打了个线性筛把素数表找到，然后开始打表找规律，发现对于每个素数，答案是循环的。于是把每个素数的循环长度和每一位的值存下来，再对x模上每一个素数的循环长度（我不知道大概有多少位所以都只读取每一位的数），本来想着能拿点分就好，结果AC了。正解是费马小定理降幂，对于指数x模上(p-1)，可以把x降得较小，然后快速幂取模。
正解

#include<cstdio>
#include<cstring>
const int N=1e5+10;
char num[N];
int qpow(int a,int x,int p)
{
	int ret=1;
	while(x)
	{
		if(x&1)ret=(ret*a)%p;
		a=(a*a)%p;
		x>>=1;
	}
	return ret%p;
}
int main()
{
	int p;
	scanf("%s%d",num,&p);
	int nowx=0;
	for(int i=0;i<strlen(num);i++)
	nowx=(nowx*10+num[i]-'0')%(p-1);
	printf("%d\n",qpow(2,nowx,p));
	return 0;
}

比赛AC代码

#include<cstdio>
#include<vector>
#include<cstring>
#define PB(v) push_back(v)
using namespace std;
typedef long long ll;
int prime[1010],p;
bool iscomp[1010];
vector<int>v[1010];
char x;
int xh[1010];
void primetable()
{
	for(int i=2;i<=1000;i++)
	{
		if(!iscomp[i])prime[p++]=i;
		for(int j=0;j<p&&i*prime[j]<=1000;j++)
		{
			iscomp[i*prime[j]]=1;
			if(i%prime[j]==0)break;			
		}
	}
}
int main()
{
	//freopen("cell.in","r",stdin);
	//freopen("cell.out","w",stdout);
	memset(xh,0,sizeof(xh));
	memset(iscomp,0,sizeof(iscomp));
	primetable();
	int y;
	for(int i=1;i<p;i++)
	{
		int base=1;
		for(int j=1;j<=1000;j++)
		{
			base*=2;
			base%=prime[i];
			if(v[prime[i]].size()>=1&&base==v[prime[i]][0])break;
			v[prime[i]].PB(base);
		}
	}
	while(scanf("%c",&x)&&x!=' ')
	{
		for(int i=1;i<p;i++)
		{
			xh[prime[i]]=xh[prime[i]]*10+x-'0';
			xh[prime[i]]%=v[prime[i]].size();
		}
	}
	scanf("%d",&y);
	if(xh[y]==0)printf("%d\n",v[y][v[y].size()-1]);
	else printf("%d\n",v[y][xh[y]-1]);
	return 0;
}

---

B.黑白图 spfa+kruskal

题目链接

---

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e5+10,M=1e6+10;
template<typename tp>inline void read(tp &x)
{
	int f=0;x=0;char ch=getchar();
	while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
	while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	if(f)x=-x;
}
int n,m,hd[N],tot,fa[N],vis[N],cnt,num;
ll dis[N],ans;
struct Edge{
    int v,w,nx;
}e[M];
struct EDGE{
	int u,v,w;
	bool operator <(const EDGE&rhs)const{return w<rhs.w;}
}edge[M];
inline int Find(int x){return fa[x]==x?x:fa[x]=Find(fa[x]);}
inline void add(int u,int v,int w)
{
	e[++tot].v=v;e[tot].w=w;e[tot].nx=hd[u];hd[u]=tot;
}
inline void spfa()
{
	memset(dis,127,sizeof(dis));
	memset(vis,0,sizeof(vis));
	dis[0]=0;queue<int>q;while(q.size())q.pop();q.push(0);
	while(q.size())
	{
		int u=q.front();q.pop();vis[u]=0;
		for(int i=hd[u];i;i=e[i].nx)
		{
			int v=e[i].v;
			if(dis[v]>dis[u]+e[i].w)
			{
				dis[v]=dis[u]+e[i].w;
				if(!vis[v])vis[v]=1,q.push(v);
			}
		}
	}
}
int main()
{
	//freopen("in.txt","r",stdin);
	read(n);read(m);
	for(int i=1;i<=n;++i)
	{
		int col;read(col);if(col==1)add(0,i,0);
	}
	for(int i=1;i<=m;++i)
	{
		int u,v,w;read(u);read(v);read(w);add(u,v,w);add(v,u,w);
	}
	spfa();bool flag=1;
	for(int i=1;i<=n&&flag;++i)
	    if(dis[i]>1ll<<60ll)flag=0;
	if(!flag)
	{
		puts("impossible");return 0;
	}
	for(int i=1;i<=n;++i)
	    for(int j=hd[i];j;j=e[j].nx)
	        if(dis[e[j].v]+e[j].w==dis[i])
	            edge[++cnt].u=i,edge[cnt].v=e[j].v,edge[cnt].w=e[j].w;
	sort(edge+1,edge+cnt+1);num=n;
	for(int i=0;i<=n;++i)
	    fa[i]=i;
	for(int i=1;i<=cnt&&num;++i)
	    if(Find(edge[i].u)!=Find(edge[i].v))
	        ans+=edge[i].w,fa[Find(edge[i].u)]=Find(edge[i].v),--num;
	printf("%lld\n",ans);
	return 0;
}

总结

建立超级源连接所有黑点很巧妙，这样保证白点和黑点间的连接是最短路。最后在最小连接间的新图求一个最小生成树。

---

C.连通能力 树的直径+树形DP

题目链接

---

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+100;
inline int Read()
{
	int x=0,f=0;char ch=getchar();
	while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
	while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	if(f)x=-x;return x;
}
inline void write(int x)
{
	int buf[40],p=0;
	if(x<0)putchar('-'),x=-x;
	do{
		buf[p++]=x%10;x/=10;
	}while(x);
	for(int i=p-1;i+1;i--)putchar(buf[i]+48);
	putchar('\n');
}
int n,tot=1,hd[N],fa[N],dis[N];
int d[N]/*直径*/,mx[N]/*最长链*/,mi[N]/*次长*/,ed1[N]/*最长链端点*/,ed2[N]/*次长端点*/,r[N]/*半径*/;
struct Edge{
	int v,w,nx;
}e[N<<1];
inline void add(int u,int v,int w)
{
	e[++tot].v=v;e[tot].w=w;e[tot].nx=hd[u];hd[u]=tot;
}
void dfs(const int x)
{
	mx[x]=mi[x]=d[x]=r[x]=0;ed1[x]=ed2[x]=x;
	for(int i=hd[x];i;i=e[i].nx)
	{
		const int &v=e[i].v,w=e[i].w;
		if(v==fa[x])continue;
		fa[v]=x;dis[v]=dis[x]+e[i].w;dfs(v);
		//更新直径
		if(d[v]>d[x])d[x]=d[v],r[x]=r[v];
		if(mx[v]+w>mx[x])
		{
			mi[x]=mx[x];ed2[x]=ed1[x];
			mx[x]=mx[v]+w;ed1[x]=ed1[v];
		} 
		else if(mx[v]+w>mi[x])mi[x]=mx[v]+w,ed2[x]=ed1[v];
	}
	if(mx[x]+mi[x]>d[x])//如果直径经过根 
	{
		d[x]=mx[x]+mi[x];
		while(ed1[x]!=x&&mi[x]+dis[fa[ed1[x]]]-dis[x]>=(d[x]+1)/2)
		//在最长路上找一个点使得该点到次长路端点距离<d/2
		ed1[x]=fa[ed1[x]];
		r[x]=mi[x]+dis[ed1[x]]-dis[x];//更新半径
		if(ed1[x]!=x&&d[x]-mi[x]-(dis[fa[ed1[x]]]-dis[x])<r[x])
		    r[x]=d[x]-mi[x]-(dis[fa[ed1[x]]]-dis[x]); 
	}
}
int main()
{
	//freopen("in.txt","r",stdin);
	int size=40<<20;//40M
	//__asm("movl %0,%%esp\n"::"r"((char*)malloc(size)+size));//调试用这个 
	__asm__("movq %0,%%rsp\n"::"r"((char*)malloc(size)+size));//提交用这个 
	n=Read();
	for(int i=1;i<n;++i)
	{
		int u,v,w;u=Read();v=Read();w=Read();
		add(u,v,w);add(v,u,w);
	}
	dis[1]=fa[1]=0;
	dfs(1);
	for(int i=1;i<=n;i++)
	    write(r[i]);
	exit(0);return 0;
} 

总结

又是贼复杂的图论题。