【刷题】洛谷 P2675 《瞿葩的数字游戏》T3-三角圣地

题目背景

国王1带大家到了数字王国的中心：三角圣地。

题目描述

不是说三角形是最稳定的图形嘛，数字王国的中心便是由一个倒三角构成。这个倒三角的顶端有一排数字，分别是1 ~ N。1 ~ N可以交换位置。之后的每一行的数字都是上一行相邻两个数字相加得到的。这样下来，最底端就是一个比较大的数字啦！数字王国称这个数字为“基”。国王1希望“基”越大越好，可是每次都自己去做加法太繁琐了，他希望你能帮他通过编程计算出这个数的最大值。但是这个值可能很大，所以请你输出它mod 10007 的结果。

任务：给定N，求三角形1~N的基的最大值再去 mod 10007。

输入输出格式

输入格式：

一个整数N

输出格式：

一个整数，表示1~N构成的三角形的最大的“基”

输入输出样例

输入样例#1：

输出样例#1：

输入样例#2：

1125

输出样例#2：

700

说明

数据：

20% 0<=N<=100

50% 0<=N<=3000

100% 0<=N<=1000000

样例解释：

1 3 4 2

4 7 6

11 13

24 是N=4的时候的最大值，当然还有别的构成形式。

PS：它叫做三角圣地，其实它就是个三角形~

本题数据已经更新，目前全部正确无误！

不要面向数据编程！

题解

把三角形画出来之后，发现，数列里一个数在“基”被加的次数就是从它的位置走向底端的方案数

然后，\(ans=C_{n-1}^0a_0+C_{n-1}^1a_1+C_{n-1}^2a_2+...+C_{n-1}^{n-1}a_{n-1}=\sum_{i=0}^{n-1}C_{n-1}^ia_i\)

所以越往中间走，被加的次数就越多，那么贪心地把数列拍好，然后 \(O(n)\) 用Lucas算就好了

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXN=1000000+10,Mod=1e4+7;
int n,A[MAXN];
ll res,fac[MAXN],inv[MAXN];
template<typename T> inline void read(T &x)
{
	T data=0,w=1;
	char ch=0;
	while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
	if(ch=='-')w=-1,ch=getchar();
	while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
	x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
	if(x<0)putchar('-'),x=-x;
	if(x>9)write(x/10);
	putchar(x%10+'0');
	if(ch!='\0')putchar(ch);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline ll qexp(ll a,ll b)
{
	ll res=1;
	while(b)
	{
		if(b&1)res=res*a%Mod;
		a=a*a%Mod;
		b>>=1;
	}
	return res;
}
inline void init()
{
	fac[0]=1;
	for(register int i=1;i<Mod;++i)fac[i]=fac[i-1]*i%Mod;
	inv[Mod-1]=qexp(fac[Mod-1],Mod-2);
	for(register int i=Mod-2;i>=0;--i)inv[i]=inv[i+1]*(i+1)%Mod;
}
inline ll C(ll n,ll m)
{
	if(n<m)return 0;
	if(n<Mod&&m<Mod)return fac[n]*inv[m]%Mod*inv[n-m]%Mod;
	else return C(n/Mod,m/Mod)*C(n%Mod,m%Mod)%Mod;
}
int main()
{
	read(n);
	init();
	if(n&1)
	{
		A[(n>>1)+1]=n;A[(n>>1)+2]=n-1;
		for(register int i=(n>>1);i>=1;--i)A[i]=A[i+1]-2;
		for(register int i=(n>>1)+3;i<=n;++i)A[i]=A[i-1]-2;
	}
	else
	{
		A[1]=1;
		for(register int i=2;i<=(n>>1);++i)A[i]=A[i-1]+2;
		A[(n>>1)+1]=n;
		for(register int i=(n>>1)+2;i<=n;++i)A[i]=A[i-1]-2;
	}
	for(register int i=1;i<=n;++i)(res+=C(n-1,i-1)*A[i]%Mod)%=Mod;
	write(res,'\n');
	return 0;
}