HDU4828 Grids 百度之星初赛1002解题报告

搞ACM起步比较晚,但是对这方面兴趣很大,第一次参加百度之星,资格赛水过一道题进入初赛,结果初赛被一道题难倒了。

初赛分两场,每场筛选500人(赛程规定),而之后的结果上看,只有约250人A出题,而大多数人都是通过这道题拿到复赛入场券。这道题我一直觉得非常接近结果了,但是在最重要的环节,也就是取模上出了很大的问题,导致没有A出该题,下面是该题的分析。

题干上,大多数人可以看出看出规律,即求卡特兰数,HDOJ上有相类似的题目,但是这道题的数据量非常大,有100万之多,因而题目上要求对1000000007取模。

卡特兰数有很多递推公式,我选择的是h(n)=h(n-1)*(4*n-2)/(n+1);这条。

由于数据量比较大,因而我选择用打表的方式,但是这里出现了一个问题,卡特兰数的数据在不断的增大,因而一定要在取模的基础上进行递推,然而简单地进行取模会造成递推公式出错,即不能在原结果上直接取模。

比赛时就卡在这点上了,怎么想都想不出改进的方法。

直到后来做到HDOJ的拓展欧几里德的题目(HDOJ 1576),才直到这道题的解法。

除法的取模需要用extgcd的算法进行取逆,这是该题的关键。

该题代码如下:


#include
using namespace std;
__int64 h[1000001];
const int  MOD=1000000007;
int extgcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int r=extgcd(b,a%b,x,y);
    int t=x;
    x=y;
    y=t-a/b*y;
    return r;
}

int main()
{
	int T;
	int n,i,j,x,y;
	h[0]=1;
	for(i=1;i<=1000000;i++)
	{
		h[i]=(4*i-2)%MOD;
		h[i]=(h[i]*h[i-1])%MOD;
		extgcd(i+1,MOD,x,y);
		h[i]=(h[i]*(x+MOD)%MOD)%MOD;
	}
	cin>>T;
	for(i=1;i<=T;i++)
	{
		cin>>n;
		cout<<"Case #"<


你可能感兴趣的:(hdu,ACM)