11.6 正睿停课训练 Day17
2018.11.6 正睿停课训练 Day17
虽然题都改了但还是咕了无数天的博客...
期望得分:忘了
实际得分:55+(80->60)+20
A chinese(思路 计数)
考虑\(\sum i*f[i]\)的意义,即所有方案中炼字的个数和。
考虑枚举每个数\(i∈[1,k]\)作为炼字时出现在哪些方案中(不是算一种方案有哪些炼字,不会算重啊)。
值\(i\)在某一个位置作为炼字时,它的方案数为\((i-1)^{n-1+m-1}*k^{(n-1)*(m-1)}\)。同时每个位置是等价的,最后乘\(n*m\)就好了。
枚举\(i\)就行了。复杂度\(O(k\log v)\)。
//587ms	512kb
#include <cstdio>
#define mod 1000000007
typedef long long LL;
inline int FP(int x,int k)
{
	int t=1;
	for(; k; k>>=1,x=1ll*x*x%mod)
		if(k&1) t=1ll*t*x%mod;
	return t;
}
int main()
{
	int n,m,K; scanf("%d%d%d",&n,&m,&K);
	if(n==1&&m==1) return printf("%d\n",K),0;
	LL ans=0;
	for(int i=2; i<=K; ++i) ans+=FP(i-1,n+m-2);
	printf("%lld\n",ans%mod*FP(K,1ll*(n-1)*(m-1)%(mod-1))%mod*n%mod*m%mod);
	return 0;
}
B physics(单调队列/剪枝 DP)
带修改(只将1修改为0)最大全1正方形。
由于数据比较水,每次修改时判一下修改点是否在当前最大矩形内,如果不在就不用DP了。然后你就A了。。
正解:
对于一次询问的情况,我们可以预处理up[i][j]表示从(i,j)往上连续的1有多少。令L[i]表示当前行在第i列左边第一个up比它小的位置,R[i]同理,可以用单调栈预处理。那么答案就是\(max{min{up[row][i],R[i]-L[i]-1}}\)。
如果有多次修改,首先可以倒序做,这样答案就是单调不降的了(不需要每次在整个矩形里求max)。
修改一个点\((x,y)\)只会影响这一列的up,down(down定义类似up),可以暴力修改。
现在只需要对\(x\)这一行算一遍就行了。但是有up和down两个值,直接求最大正方形不好算,但是可以求是否有大小为\(k\)的最大全1正方形。
存在大小为\(k\)的正方形,也就是存在一个\(i\)满足\(min{i,min{up[i-k+1,i]}+min{down[i-k+1,k]}-1}>=k\),这就是滑动窗口了。
而且答案最多只会变\(n\)次,所以对每次修改不断尝试更新答案就可以了。
复杂度\(O(n^2+nq)\)。因为数据弱跑的比暴力+剪枝慢。
确实有分治做法,但是感觉常数更大的一匹啊?(别问我怎么知道的)
//6558ms	47884kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;
int A[N][N],up[N][N],down[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
int preSolve(int n,int m)
{
	static int L[N],R[N],sk[N];
	for(int i=1; i<=n; ++i)
		for(int j=1; j<=m; ++j) up[i][j]=A[i][j]?up[i-1][j]+1:0;
	for(int i=n; i; --i)
		for(int j=1; j<=m; ++j) down[i][j]=A[i][j]?down[i+1][j]+1:0;
	int ans=0;
	for(int i=1; i<=n; ++i)
	{
		up[i][0]=up[i][m+1]=-1;
		for(int j=1,top=0; j<=m+1; ++j)
		{
			while(up[i][sk[top]]>up[i][j]) R[sk[top--]]=j;
			sk[++top]=j;
		}
		for(int j=m,top=0; ~j; --j)
		{
			while(up[i][sk[top]]>up[i][j]) L[sk[top--]]=j;
			sk[++top]=j;
		}
		for(int j=1; j<=m; ++j) ans=std::max(ans,std::min(up[i][j],R[j]-L[j]-1));
	}
	return ans;
}
void Modify(int n,int c)
{
	for(int i=1; i<=n; ++i) up[i][c]=A[i][c]?up[i-1][c]+1:0;
	for(int i=n; i; --i) down[i][c]=A[i][c]?down[i+1][c]+1:0;
}
bool Update(int m,int r,int K)
{
	static int q1[N],v1[N],q2[N],v2[N];
	int *up=::up[r],*down=::down[r];
	for(int i=1,h1=1,t1=0,h2=1,t2=0; i<=m; ++i)
	{
		while(h1<=t1 && i-q1[h1]>=K) ++h1;
		while(h1<=t1 && up[i]<=v1[t1]) --t1;
		q1[++t1]=i, v1[t1]=up[i];
		while(h2<=t2 && i-q2[h2]>=K) ++h2;
		while(h2<=t2 && down[i]<=v2[t2]) --t2;
		q2[++t2]=i, v2[t2]=down[i];
		if(std::min(i,v1[h1]+v2[h2]-1)>=K) return 1;
	}
	return 0;
}
int main()
{
	static int Ans[10005],qx[10005],qy[10005];//not N...
	int n=read(),m=read(),Q=read();
	for(int i=1; i<=n; ++i)
	{
		register char c=gc(); while(c!='+'&&c!='-') c=gc();
		A[i][1]=c=='+';
		for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
	}
	for(int i=1; i<=Q; ++i) A[qx[i]=read()][qy[i]=read()]=0;
	int ans=preSolve(n,m);
	for(int i=Q; i; --i)
	{
		Ans[i]=ans, A[qx[i]][qy[i]]=1, Modify(n,qy[i]);
		while(Update(m,qx[i],ans+1)) ++ans;
	}
	for(int i=1; i<=Q; ++i) printf("%d\n",Ans[i]);
	return 0;
}
暴力+剪枝:
//1750ms	32068kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;
int A[N][N],f[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
int main()
{
	int n=read(),m=read(),Q=read();
	for(int i=1; i<=n; ++i)
	{
		register char c=gc(); while(c!='+'&&c!='-') c=gc();
		A[i][1]=c=='+';
		for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
	}
	int ansx,ansy,ans=0;
	for(int i=1; i<=n; ++i)
		for(int j=1; j<=m; ++j)
			f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0,
			f[i][j]>ans&&(ans=f[i][j],ansx=i,ansy=j);
	for(int x,y; Q--; )
	{
		A[x=read()][y=read()]=0;
		if(x>=ansx-ans&&x<=ansx&&y>=ansy-ans&&y<=ansy)
		{
			ans=0;
			for(int i=1; i<=n; ++i)
				for(int j=1; j<=m; ++j)
					f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0,
					f[i][j]>ans&&(ans=f[i][j],ansx=i,ansy=j);
		}
		printf("%d\n",ans);
	}
	return 0;
}
C chemistry(期望 DP)
之前怎么没写思路啊。。
过了这么久不想写了,粘solution了。
感觉挺nb的。。








//364ms	17200kb
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mod 1000000007
#define Mod(x) x>=mod&&(x-=mod)
typedef long long LL;
const int N=2e5+5;
int n,K,P,A[N],C[13][13],Enum,H[N],nxt[N<<1],to[N<<1],f[N],g[N][11];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
inline int FP(int x,int k)
{
	int t=1;
	for(; k; k>>=1,x=1ll*x*x%mod)
		if(k&1) t=1ll*t*x%mod;
	return t;
}
inline void AE(int u,int v)
{
	to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
	to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
void DFS(int x,int fa)
{
	int pw=P;
	for(int i=1; i<=K; ++i) g[x][i]=pw=1ll*pw*A[x]%mod;//f[x][i]=g[x][i]=P*A[x]^i
	f[x]=pw;
	for(int i=H[x],v; i; i=nxt[i])
		if((v=to[i])!=fa)
		{
			DFS(v,x);
			LL tmp=f[x]+f[v];
			for(int j=1; j<K; ++j) tmp+=1ll*C[K][j]*g[x][j]%mod*g[v][K-j]%mod;
			f[x]=tmp%mod;
			for(int j=K; j; --j)
			{
				tmp=1ll*P*g[v][j]+g[x][j];
				for(int k=1; k<j; ++k)
					tmp+=1ll*C[j][k]*g[x][k]%mod*g[v][j-k]%mod;
				g[x][j]=tmp%mod;
			}
		}
}
int main()
{
	n=read(),K=read(),P=1ll*read()*FP(read(),mod-2)%mod;
	for(int i=1; i<=n; ++i) A[i]=read();
	for(int i=1; i<n; ++i) AE(read(),read());
	C[0][0]=1;
	for(int i=1; i<=K; ++i)
	{
		C[i][0]=C[i][i]=1;
		for(int j=1; j<i; ++j) C[i][j]=C[i-1][j]+C[i-1][j-1], Mod(C[i][j]);
	}
	DFS(1,1);
	printf("%d\n",f[1]);
	return 0;
}
考试代码
A
#include <cstdio>
#include <cstring>
#include <algorithm>
#define mod 1000000007
typedef long long LL;
const int N=5;
int n,m,K,f[N*N],A[N][N];
bool OK(int x,int y)
{
	int a=A[x][y];
	for(int i=1; i<=n; ++i) if(A[i][y]>=a&&i!=x) return 0;
	for(int i=1; i<=m; ++i) if(A[x][i]>=a&&i!=y) return 0;
	return 1;
}
void Calc()
{
	static bool vis[N];
	memset(vis,0,sizeof vis);
	int ans=0;
	for(int i=1; i<=n; ++i)
		for(int j=1; j<=m; ++j)
			if(!vis[j]&&OK(i,j))
				{vis[j]=1, ++ans; break;}
	++f[ans];
}
void DFS(int x,int y)
{
	if(y>m) ++x,y=1;
	if(x>n)
	{
		Calc();
		return;
	}
	for(int i=1; i<=K; ++i) A[x][y]=i, DFS(x,y+1);
}
void Solve()
{
	DFS(1,1);
	LL ans=0;
	for(int i=0; i<=n*m; ++i) ans+=1ll*i*f[i]%mod;
	for(int i=0; i<=n*m; ++i) printf("%d ",f[i]); puts("");
	printf("n:%d m:%d K:%d ans:",n,m,K);
	printf("%d\n",(int)(ans%mod));
//	printf("%d,",(int)(ans%mod));
}
int main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
//	for(n=2; n<=2; ++n)
//		for(m=2; m<=2; ++m)
//			for(K=1; K<=15; ++K) Solve();
//	return 0;
	int n,m,K; scanf("%d%d%d",&n,&m,&K); ::n=n, ::m=m, ::K=K;
	if(n==4&&m==4&&K>=3)
	{
		if(K==3) printf("%d\n",20470320);
		else if(K==4) printf("%d\n",330277355);
//		else if(K==5) printf("%d\n",);
		return 0;
	}
	DFS(1,1);
	LL ans=0;
	for(int i=0; i<=n*m; ++i) ans+=1ll*i*f[i]%mod;
	printf("%d\n",(int)(ans%mod));
	return 0;
}
B
写分治T掉暴力能得的20分
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2003;
int A[N][N],f[N][N],mxu[N][N],mxd[N][N],mxl[N][N],mxr[N][N],ans[10000005][2],Ans;
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
int Calc(int n,int m)
{
	int ans=0;
	for(int i=1; i<=n; ++i)
		for(int j=1; j<=m; ++j)
			ans=std::max(ans,f[i][j]=A[i][j]?std::min(f[i-1][j-1],std::min(f[i-1][j],f[i][j-1]))+1:0);
	return ans;
}
void Violence(int n,int m,int Q)
{
	while(Q--) A[read()][read()]=0, printf("%d\n",Calc(n,m));
}
inline int Go_left(int x,int y)
{
	int ans=0;
	while(A[x][y-ans]) ++ans;
	return ans;
}
inline int Go_right(int x,int y)
{
	int ans=0;
	while(A[x][y+ans]) ++ans;
	return ans;
}
inline int Go_up(int x,int y)
{
	int ans=0;
	while(A[x-ans][y]) ++ans;
	return ans;
}
inline int Go_down(int x,int y)
{
	int ans=0;
	while(A[x+ans][y]) ++ans;
	return ans;
}
inline bool OK(int *A,int *A2,int l,int r,int x)
{
	static int q[N],p[N],q2[N],p2[N];
	int h=1,t=0,h2=1,t2=0;
	for(int i=l; i<=r; ++i)
	{
		while(h<=t&&q[t]>=A[i]) --t;
		q[++t]=A[i], p[t]=i;
		while(h2<=t2&&q2[t2]>=A2[i]) --t2;
		q2[++t2]=A2[i], p2[t2]=i;
		while(p[h]<=i-x) ++h;
		while(p2[h2]<=i-x) ++h2;
		if(i>=x&&q[h]+q2[h2]-1>=x) return 1;
	}
	return 0;
}
int Row(int row,int l,int r)
{
	int L=2,R=r-l+1,mid,ans=1;
	while(L<=R)
		if(OK(mxu[row],mxd[row],l,r,mid=L+R>>1)) ans=mid,L=mid+1;
		else R=mid-1;
	return ans;
}
inline bool OK2(int col,int l,int r,int x)
{
	static int q[N],p[N],q2[N],p2[N];
	int h=1,t=0,h2=1,t2=0;
	for(int i=l; i<=r; ++i)
	{
		while(h<=t&&q[t]>=mxl[i][col]) --t;
		q[++t]=mxl[i][col], p[t]=i;
		while(h2<=t2&&q2[t2]>=mxr[i][col]) --t2;
		q2[++t2]=mxr[i][col], p2[t2]=i;
		while(p[h]<=i-x) ++h;
		while(p2[h2]<=i-x) ++h2;
		if(i>=x&&q[h]+q2[h2]-1>=x) return 1;
	}
	return 0;
}
int Colu(int col,int u,int d)
{
	int L=2,R=d-u+1,mid,ans=1;
	while(L<=R)
		if(OK2(col,u,d,mid=L+R>>1)) ans=mid,L=mid+1;
		else R=mid-1;
	return ans;
}
int Pre(int u,int d,int l,int r,int rt)
{
	if(u>d||l>r) return 0;
//	printf("Pre(%d~%d %d~%d %d)\n",u,d,l,r,rt);
	if(d-u>=r-l)
	{
		int mid=u+d>>1;
//		printf("mid:%d\n",mid);
		for(int i=l; i<=r; ++i) mxu[mid][i]=Go_up(mid,i), mxd[mid][i]=Go_down(mid,i);
//		for(int i=l; i<=r; ++i) printf("mxu[%d][%d]=%d mxd[%d][%d]=%d\n",mid,i,mxu[mid][i],mid,i,mxd[mid][i]);
		int res=Row(mid,l,r);
		ans[rt][0]=Pre(u,mid-1,l,r,rt<<1), ans[rt][1]=Pre(mid+1,d,l,r,rt<<1|1);
//		printf("1 mid:%d Pre(%d~%d %d~%d %d) res:%d %d %d\n",mid,u,d,l,r,rt,res,ans[rt][0],ans[rt][1]);
		return std::max(res,std::max(ans[rt][0],ans[rt][1]));
	}
	else
	{
		int mid=l+r>>1;
		for(int i=u; i<=d; ++i) mxl[i][mid]=Go_left(i,mid), mxr[i][mid]=Go_right(i,mid);
		int res=Colu(mid,u,d);
		ans[rt][0]=Pre(u,d,l,mid-1,rt<<1), ans[rt][1]=Pre(u,d,mid+1,r,rt<<1|1);
//		printf("2 mid:%d Pre(%d~%d %d~%d %d) res:%d %d %d\n",mid,u,d,l,r,rt,res,ans[rt][0],ans[rt][1]);
		return std::max(res,std::max(ans[rt][0],ans[rt][1]));
	}
}
void Solve(int u,int d,int l,int r,int x,int y,int rt)
{
	if(u>d||l>r) return;
	if(d-u>=r-l)
	{
		int mid=u+d>>1;
		if(x>mid)
		{
			for(int i=l; i<=r; ++i) mxd[mid][i]=Go_down(mid,i);
			Ans=std::max(Ans,std::max(Row(mid,l,r),ans[rt][0])), Solve(mid+1,d,l,r,x,y,rt<<1|1);
		}
		else if(x<mid)
		{
			for(int i=l; i<=r; ++i) mxu[mid][i]=Go_up(mid,i);
			Ans=std::max(Ans,std::max(Row(mid,l,r),ans[rt][1])), Solve(u,mid-1,l,r,x,y,rt<<1);
		}
		else mxu[x][y]=mxd[x][y]=0, Ans=std::max(Ans,std::max(ans[rt][0],ans[rt][1]));
	}
	else
	{
		int mid=l+r>>1;
		if(y>mid)
		{
			for(int i=u; i<=d; ++i) mxl[i][mid]=Go_left(i,mid);
			Ans=std::max(Ans,std::max(Colu(mid,u,d),ans[rt][0])), Solve(u,d,mid+1,r,x,y,rt<<1|1);
		}
		else if(y<mid)
		{
			for(int i=u; i<=d; ++i) mxr[i][mid]=Go_right(i,mid);
			Ans=std::max(Ans,std::max(Colu(mid,u,d),ans[rt][1])), Solve(u,d,l,mid-1,x,y,rt<<1);
		}
		else mxl[x][y]=mxr[x][y]=0, Ans=std::max(Ans,std::max(ans[rt][0],ans[rt][1]));
	}
}
int main()
{
//	freopen("ex_physics1.in","r",stdin);
//	freopen(".out","w",stdout);
	int n=read(),m=read(),Q=read();
	for(int i=1; i<=n; ++i)
	{
		register char c=gc(); while(c!='+'&&c!='-') c=gc();
		A[i][1]=c=='+';
		for(int j=2; j<=m; ++j) A[i][j]=gc()=='+';
	}
	if(1ll*n*m*Q<=2e8) return Violence(n,m,Q),0;
	Pre(1,n,1,m,1);
	for(int x,y; Q--; )
		x=read(),y=read(),A[x][y]=0,Ans=1,Solve(1,n,1,m,x,y,1),printf("%d\n",Ans);
	return 0;
}
C
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mod 1000000007
#define Mod(x) x>=mod&&(x-=mod)
typedef long long LL;
const int N=2e5+5;
int n,K,P,A[N],pw[N],pw2[N],Enum,H[N],nxt[N<<1],to[N<<1],f[105][105],g[105][105],sz[N];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
inline int FP(int x,int k)
{
	int t=1;
	for(; k; k>>=1,x=1ll*x*x%mod)
		if(k&1) t=1ll*t*x%mod;
	return t;
}
inline void AE(int u,int v)
{
	to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
	to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
namespace Subtask1
{
	const int N=204;
	int n,K,f[N];
	LL Ans,Sum;
	bool chose[N];
	void DFS2(int x,int fa)
	{
		if(chose[x])
		{
			f[x]=0;
			for(int i=H[x],v; i; i=nxt[i])
				if((v=to[i])!=fa) DFS2(v,x);
		}
		else
		{
			f[x]=A[x];
			for(int i=H[x],v; i; i=nxt[i])
				if((v=to[i])!=fa) DFS2(v,x), f[x]+=f[v], Mod(f[x]);
			if(chose[fa]) Sum+=FP(f[x],K);
		}
	}
	void Calc(int t)
	{
		chose[0]=1, Sum=0, DFS2(1,0), Ans+=Sum%mod*pw[t]%mod*pw2[n-t]%mod;
	}
	void DFS(int x,int t)
	{
		if(x>n) {Calc(t); return;}
		DFS(x+1,t), chose[x]=1, DFS(x+1,t+1), chose[x]=0;
	}
	void Main()
	{
		n=::n, K=::K, DFS(1,0), printf("%d",(int)(Ans%mod));
	}
}
void DFS(int x,int fa)
{
	static int tmp[105];
	f[x][1]=A[x], sz[x]=1;
	for(int i=H[x],v; i; i=nxt[i])
		if((v=to[i])!=fa)
		{
			DFS(v,x);
			memset(tmp,0,sizeof tmp);
			for(int j=0; j<=sz[x]; ++j)
				for(int k=1; k<=sz[v]; ++k)
					tmp[j+k]+=((g[x][j]+g[v][k])%mod+FP(f[v][k],K))%mod, Mod(tmp[j+k]);
			memcpy(g[x],tmp,sizeof tmp);
			memset(tmp,0,sizeof tmp);
			for(int j=1; j<=sz[x]; ++j)
				for(int k=0; k<=sz[v]; ++k)
					printf("tmp[%d]=f[x][%d]+f[v][%d]=%d\n",j+k,j,k,f[x][j]+f[v][k]),
					tmp[j+k]+=f[x][j]+f[v][k], Mod(tmp[j+k]);
			memcpy(f[x],tmp,sizeof tmp);
			sz[x]+=sz[v];
			printf("%d->%d\n",x,v);
			for(int i=1; i<=n; ++i) printf("f[%d][%d]=%d g[%d][%d]=%d\n",x,i,f[x][i],x,i,g[x][i]);
		}
}
int main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	n=read(),K=read(),P=1ll*read()*FP(read(),mod-2)%mod;
	pw[0]=pw2[0]=1, pw[1]=mod+1-P,pw2[1]=P;//pw:被活化的概率
	for(int i=1; i<=n; ++i) pw[i]=1ll*pw[i-1]*(mod+1-P)%mod, pw2[i]=1ll*pw2[i-1]*P%mod;
	for(int i=1; i<=n; ++i) A[i]=read();
	for(int i=1; i<n; ++i) AE(read(),read());
	if(n<=200) return Subtask1::Main(),0;
//	DFS(1,1);
//	LL ans=0;
//	for(int i=1; i<=n; ++i) printf("%d:%d %d\n",i,f[1][i],g[1][i]);
//	for(int i=1; i<=n; ++i) ans+=1ll*pw2[i]*pw[n-i]%mod*(FP(f[1][i],K)+g[1][i])%mod;
//	printf("%d\n",(int)(ans%mod));
	return 0;
}
11.6 正睿停课训练 Day17的更多相关文章
- 11.5 正睿停课训练 Day16
		目录 2018.11.5 正睿停课训练 Day16 A 道路规划(思路) B 逻辑判断(枚举 位运算/DP 高维前缀和) C 区间(贪心/树状数组) 考试代码 A B C 2018.11.5 正睿停课 ... 
- 11.2 正睿停课训练 Day15
		目录 2018.11.2 正睿停课训练 Day15 A 郁闷的小G(二分) B 小G的树(树形DP) C 数的距离(思路) 考试代码 B C 2018.11.2 正睿停课训练 Day15 时间:3.5 ... 
- 11.1 正睿停课训练 Day14
		目录 2018.11.1 正睿停课训练 Day14 A 字符串 B 取数游戏(贪心) C 魔方(模拟) 考试代码 B C 2018.11.1 正睿停课训练 Day14 时间:3.5h 期望得分:100 ... 
- 10.29 正睿停课训练 Day11
		目录 2018.10.29 正睿停课训练 Day11 A 线段树什么的最讨厌了(思路 DFS) B 已经没有什么好害怕的了(差分 前缀和) C 我才不是萝莉控呢(DP 贪心 哈夫曼树) 考试代码 A ... 
- 10.30 正睿停课训练 Day12
		目录 2018.10.30 正睿停课训练 Day12 A 强军战歌(DP 树状数组 容斥) B 当那一天来临(思路) C 假如战争今天爆发(贪心) 考试代码 B C 2018.10.30 正睿停课训练 ... 
- 10.31 正睿停课训练 Day13
		目录 2018.10.31 正睿停课训练 Day13 A Poker(期望) B Label(高斯消元) C Coin(二分图染色 博弈) 考试代码 A(打表) B 2018.10.31 正睿停课训练 ... 
- 10.25 正睿停课训练 Day9
		目录 2018.10.25 正睿停课训练 Day9 A 数独(思路 DP) B 红绿灯(最短路Dijkstra) C 轰炸(计算几何 圆并) 考试代码 B C 2018.10.25 正睿停课训练 Da ... 
- 10.24 正睿停课训练 Day8 AM
		目录 2018.10.24 正睿停课训练 Day8 AM A 棒棒糖(组合) B 彩虹糖(思路 博弈) C 泡泡糖(DP) 考试代码 A B C 2018.10.24 正睿停课训练 Day8 AM 期 ... 
- 10.23 正睿停课训练 Day7
		目录 2018.10.23 正睿停课训练 Day7 A 矩形(组合) B 翻转(思路) C 求和(思路 三元环计数) 考试代码 B1 B2 C 2018.10.23 正睿停课训练 Day7 期望得分: ... 
随机推荐
- Imperial roads 非严格次小生成树
			cf测评姬比uva快了五倍... /* 不管这条边是不是在mst上,直接跑lca求出路径上的最大边w即可 ans=mst-w+dist(u,v) */ #include<bits/stdc++. ... 
- Python变量的作用域
			局部变量 局部变量是指在函数内部定义并使用的变量,他只在函数内部有效.即函数内部的名字只在函数运行时才会创建,在函数运行之前或者运行完毕之后,所有的名字就都不存在了.所以,如果在函数外部使用函数内部定 ... 
- shell 批量替换内容
			一.查找/data/www/vhosts/WOStest2_YiiFramework/console/sh/下所有.sh文件 把/data/www/vhosts/WOS_YiiFramework 替换 ... 
- Python+selenium自动化测试中Windows窗口跳转方法
			Python+selenium自动化测试中Windows窗口跳转方法 #第一种方法 #获得当前窗口 nowhandle=driver.current_window_handle #打开弹窗 drive ... 
- 论文阅读笔记二十二:End-to-End Instance Segmentation with Recurrent Attention(CVPR2017)
			论文源址:https://arxiv.org/abs/1605.09410 tensorflow 代码:https://github.com/renmengye/rec-attend-public 摘 ... 
- 集腋成裘-02-css基础-01
			CSS 层叠样式表(Cascading Style Sheets)(级联样式表) 1 选择器 1.1 写法 选择器是一个选择标签的过程. 选择器{属性:值;...} 1.2 基础选择器 标签选择器 类 ... 
- Apache Kafka 概述
			kafka教程,完全参照w3school: https://www.w3cschool.cn/apache_kafka/apache_kafka-dac11yot.html 以下是入门学习过程中摘录的 ... 
- PAT Basic 1069. 微博转发抽奖(20)
			小明PAT考了满分,高兴之余决定发起微博转发抽奖活动,从转发的网友中按顺序每隔N个人就发出一个红包.请你编写程序帮助他确定中奖名单. 输入格式: 输入第一行给出三个正整数M(<= 1000).N ... 
- SQL Server表关联
			表关联:Hash.Nested Loops.Merge.这是实际算法,不是T-SQL中的inner/left/right/full/cross join.优化器会把这些T-SQL写法转换成上面的3种算 ... 
- Vue axios 返回数据绑定到vue对象问题
			在项目中需要用到后台的数据对前端渲染,使用到了vue整合的axios,使用vue中的钩子函数在页面组件挂载完成之后向后台发送一个get请求然后将返回后的数据赋值data()中定义的属性: 执行后前端报 ... 
