【NOIP2016提高A组模拟10.15】算循环

Description

意思就是让你运行这个程序输出结果

Input

n和m（$n,m<=10^{18}$）

Output

答案，tmp

Sample Input

167 198

Sample Output

906462341

Solution

首先对题目中的程序搞一搞
改变顺序：有关n的放一起，有关m的放一起
拆开：分别算出有关n的答案和有关m的答案，然后乘起来
继续拆开，套用公式搞定

$$\sum_{i=1}^n i^2={{n(n+1)(2n+1)} \over6}$$
为了应对Cold_Chair对我的步步紧逼，这里写一个公式的简单证明过程：
假设公式成立，那么对于第n项的公式就是 ${{n(n+1)(2n+1)} \over6}$
然后加 $(n+1)^2$
然后展开，加进去，再化简，就变成了 ${{(n+1)(n+2)(2(n+1)+1)} \over6}$
又变回了原公式

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
#define fo(i,a,b) for(ll i=a;i<=b;i++)
#define mo 1000000007ll
#define ny2 500000004
#define ny6 166666668
using namespace std;
ll hh(ll n)
{
    ll an=0;
    an=(n*n)%mo;an=(an*n)%mo;
    an=(an-((((((2*n)*(n-2))%mo)*ny2)%mo)*n)%mo+mo)%mo;
    an=(an+(((n-1)*(n-2))%mo)*ny2)%mo;
    n=(n-2);
    an=(an+(((((n*(n+1))%mo)*(2*n+1))%mo)*ny6))%mo;
    an=(an*ny2)%mo;
    return an%mo;
}
int main()
{
    freopen("loop.in","r",stdin);freopen("loop.out","w",stdout);
    ll n,m,an=0;
    scanf("%lld%lld",&n,&m);
    an=(hh(n%mo)*hh(m%mo))%mo;
    printf("%lld",an);
    fclose(stdin);fclose(stdout);
}