Gym - 101981G The 2018 ICPC Asia Nanjing Regional Contest—— G.Pyramid（找规律。。）

问题虫洞： G - Pyramid

 

黑洞内窥：

输入一个数n，求n阶等边三角形中有多少个（可以是任意规格的，只要是等边）等边三角形

例如三阶等边三角形中的其中一个等边三角形的亚子：

 

思维光年：

数呗，总的三角形数目sum = 规则的三角形数目a + 不规则的三角形数目b （上图a为规则的，b为不规则的）

规则的三角形数目a = 正三角形的数目a1 + 倒三角形的数目a2

递推可得：a1[1] = 1, a2[1] = 0

                 a1[n] = a1[n-1] + C(n+1, 2)

                 a2[n] = a2[n-1] + C((n+1)/2, 2) + C((n+2)/2, 2)

                 ( C(n, 2) = n * (n-1) / 2 )

所以：    a[n] = a1[n] + a2[n];

不规则三角形的数目b = 正六边形的数目*2 + 对称的三边拉长后的六边形的数目*2 + 倒对称的三边拉长后的六边形的数目*2

由于BZ水平有限，未能在比赛中推出此b的递推公式，，，导致我队该题凉凉。。。。。

（后来想了一想，数据都开到1e9了，递推你妹啊  ， 这个方法确实不太靠谱，，，，）

 

理性的求解：

打表找规律。。。。。

（鬼知道答案就是C((n+3), 4)%1e9+7 哦~~~~~）

这里推荐一篇博客吧：Gym - 101981G The 2018 ICPC Asia Nanjing Regional Contest G.Pyramid 找规律

sum = (n)*(n+1)*(n+2)*(n+3)/24 = C(n+3,4)

 

ACcode：

#include<stdio.h>
#include<iostream>
#include<map>
#include<algorithm>
#include<cstring>
#include<string.h>
#include<math.h>
#include<vector>
#include<map>
using namespace std;
typedef long long ll;
#define MAXN 1000005
#define INF 0x3f3f3f3f//将近int类型最大数的一半，而且乘2不会爆int
#define MOD 1000000007 // MOD%4 = 3
const double pi = acos(-1.0);
const double eps = 1e-6;

ll qpow(ll a, ll b, ll c)
{
    ll ans = 1;
    while(b)
    {
        if(b&1)
            ans = ans*a%c;
        a = a*a%c;
        b >>= 1;
    }
    return ans;
}

int main()
{
    ll t, n;
    scanf("%lld", &t);
    ll inv24 = qpow(24, MOD-2, MOD);
    while(t--)
    {
        scanf("%lld", &n);
        ll ans = n*(n+1)%MOD*(n+2)%MOD*(n+3)%MOD*inv24%MOD;
        cout << ans << '\n';
    }
   return 0;
}