NOIP2019复赛模拟11-day2

T1

Trape
出题人真的是太有毒了！！！
题意在搞什么！！！QAQ
（梯形只存在横向放着，没有乱七八糟旋转）

题解：

枚举下底边，使左右边界不断右移即可。（可以使用单调队列）

#include <bits/stdc++.h>
using namespace std;
int a[5200][5200];
int main()
{
	int i,j,ans=0,n,m,x,y;
	scanf("%d%d",&n,&m);
	memset(a,0,sizeof(a));
	while(m--)
	{
		scanf("%d%d",&x,&y);
		x+=2500,y+=2500;
		for(i=0; i<=50; i++)
			for(j=(i+1)/2; j<=75-(i+1)/2; j++)
				ans=max(ans,++a[x+j][y-i]);
	}
	cout<<ans;
	return 0;
}

---

T2

Rcomb
一开始掉坑里了，满脑子数学思想。
冷静下来后。。。发现太水了。<(￣ ﹌ ￣)>
直接DP从i位开始，计算每个点的贡献。最后再乘以权值即可。
（要注意头尾，和中间是不等价的）

题解：

• 然后可以发现二维可以转成一维，（每一位只与前面一位有关）
• dalao们也可以$o(n)$找规律

#include <bits/stdc++.h>
using namespace std;
double a[1000010],b[1000010];
int main()
{
	int n,i,j;
	double ans=0;
	scanf("%d",&n);
	for(i=1; i<=n; i++) scanf("%lf",&a[i]);
	b[1]=0;
	for(i=2; i<=n; i++)
	{
		b[i]=b[i-1]+1.0/(i-1);
		for(j=i-1; j>=2; j--) b[j]=(b[j-1]*(j-1)+b[j]*(i-j)+2)/(i-1);
		b[1]=b[1]+1.0/(i-1);
	}
	for(i=1; i<=n; i++) ans+=a[i]*b[i];
	printf("%.5lf\n",ans);
	return 0;
}

T3

Game2
出题人可能觉得之前的题目太水了，出了一道其实很水，但很毒瘤的题
此题用到 SG ，anti-SG，SAM，倍增
<(￣ ﹌ ￣)>

题解：

• 先对S建一个SAM，然后计算后缀自动机上每一个状态点的SG值，然后倍增跑一遍
• 对于每一个l，r，先在SAM上找到他对应的状态点（拿这个子串在上面跑，最后到达的点），然后异或一下即可。
• SG（谁先走完最后一步谁胜）：异或值不为0，则先手必胜；为0，则后手必胜
• anti-SG（谁先走完最后一步谁输）：
  • 若异或值不为0，且有SG值大于1
  • 若异或值为0 ，且每个SG值都等于1（若有为0的可以忽略）
  • 以上两种先手必胜

#include <bits/stdc++.h>
using namespace std;
char ch[1000010];
int cnt=0,S,T,n,h[1000010],id[4000010];
int w[100],pos[1000010];
struct node
{
	int s[30];
	int max_len,SG,link[30];
} a[4000010];
int add(int x)
{
	int now=(++cnt),pre=T;
	T=cnt;
	a[now].max_len=a[pre].max_len+1;
	for(; pre&&a[pre].s[x]==0; pre=a[pre].link[0]) a[pre].s[x]=now;
	if(!pre) return (a[now].link[0]=S);
	int y=a[pre].s[x];
	if(a[y].max_len==a[pre].max_len+1) return (a[now].link[0]=y);
	int next=(++cnt);
	a[next]=a[y];
	a[next].max_len=a[pre].max_len+1;
	a[y].link[0]=a[now].link[0]=next;
	for(; pre&&a[pre].s[x]==y; pre=a[pre].link[0]) a[pre].s[x]=next;
}
int main()
{
	int n,m,i,j,x,l,r;
	scanf("%s%d",ch+1,&m);
	n=strlen(ch+1);
	S=T=(++cnt);
	for(i=1; i<=n; i++) add(ch[i]-'a'),pos[i]=T;
	for(i=1; i<20; i++)
		for(j=1; j<=cnt; j++)
			a[j].link[i]=a[a[j].link[i-1]].link[i-1];
	memset(h,0,sizeof(h));
	for(i=1; i<=cnt; i++) h[a[i].max_len]++;
	for(i=1; i<=n; i++) h[i]+=h[i-1];
	for(i=cnt; i>=1; i--) id[h[a[i].max_len]--]=i;
	for(i=cnt; i>=1; i--)
	{
		x=id[i];
		memset(w,0,sizeof(w));
		for(j=0; j<26; j++)
			if(a[x].s[j]) w[a[a[x].s[j]].SG]++;
		while(w[a[x].SG]) a[x].SG++;
	}
	int s=0,ok=1;
	while(m--)
	{
		scanf("%d%d",&l,&r);
		x=pos[r];
		for(i=19; i>=0; i--)
			if(a[a[x].link[i]].max_len>=(r-l+1)) x=a[x].link[i];
		s=s^a[x].SG;
		if(a[x].SG>1) ok=0;
		puts(s?"Fir":"Sec");
		if((s>0&&ok==0)||(s==0&&ok==1)) puts("Fir");
		else puts("Sec");
	}
	return 0;
}
/*
1. 全为 1 ，SG=0
2. 不全为 1，SG>0
*/