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

---

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

(File IO): input:loop.in output:loop.out

Time Limits: 1000 ms Memory Limits: 524288 KB

Description

Input

Output

Sample Input

167 198

Sample Output

906462341

Data Constraint

---

解题思路(我的)

首先我们可以将题目所求写成这样的形式：

$$\sum_{i=1}^n{\sum_{j=1}^m{\sum_{k=1}^n{\sum_{l=1}^m{\sum_{x=i}^k{\sum_{y=j}^l1}}}}}$$

20分:call题目上的标O(n^6)

40分:
首先我们可以知道：

$$\sum_{i=x1}^{y1}{\sum_{j=x2}^{y2}f(i)*g(j)}=(\sum_{i=x1}^{y1}f(i))*(\sum_{j=x2}^{y2}g(j))$$ （这很显然啊）
于是公式化为
$$\sum_{i=1}^n{\sum_{j=1}^m{\sum_{k=1}^{n-i+1}{\sum_{l=1}^{m-j+1}{\sum_{x=1}^k{\sum_{y=1}^l1}}}}}$$
—–>
$$\sum_{i=1}^n{\sum_{j=1}^m{\sum_{k=1}^{n-i+1}{\sum_{l=1}^{m-j+1}kl}}}$$

O(n^4)枚举

60分:
你要会逆元
很明显公式又可以化为

$$\sum_{i=1}^n{\sum_{j=1}^m{({\sum_{k=1}^{n-i+1}k})*({\sum_{l=1}^{m-j+1}l})}}$$
—–>
$$\sum_{i=1}^n{\sum_{j=1}^m{\frac{(n-i+2)*(n-i+1)}2*\frac{(m-j+2)*(m-j+1)}2}}$$
—–>
$$\sum_{i=1}^n{\sum_{j=1}^m{\frac{(n-i+2)*(n-i+1)*(m-j+2)*(m-j+1)}4}}$$
—–>
$$\frac{\sum_{i=1}^n{\sum_{j=1}^m{(n-i+2)*(n-i+1)*(m-j+2)*(m-j+1)}}}4$$
—–>
$$\frac{\sum_{i=1}^n{\sum_{j=1}^m{(n^2-2ni+3n+i^2-3i+2)*(m^2-2mj+3m+j^2-3j+2)}}}4$$

O(n^2)解决

80分:
再将公示转化得到

$$\frac{(\sum_{i=1}^n{(n^2-2ni+3n+i^2-3i+2)})*(\sum_{j=1}^m{(m^2-2mj+3m+j^2-3j+2)})}4$$

O(n)解决

100分:
有八十分的思路还想不到100的？除非你不会用公式……
有

$$\sum_{i=1}^ni=(1+i)*i/2$$
$$\sum_{i=1}^ni^2=n*(n+1)*(2n+1)/6$$

原式化为

$$\frac{(n^3-n^2(n+1)+3n^2+\frac{(2n+1)(n+1)n}{6}-\frac32(n+1)n+2n)*(m^3-m^2(m+1)+3m^2+\frac{(2m+1)(m+1)m}{6}- \frac32(m+1)m+2m)}4$$

—–>

$$\frac{(n^3-n^3-n^2+3n^2+\frac13n^3+\frac12n^2+\frac16n-frac32n^2-frac32n+2n)*(m^3-m^3-m^2+3m^2+\frac13m^3+\frac12m^2+\frac16m-frac32m^2-frac32m+2m)}4$$

—–>

$$\frac{(\frac13n^3+n^2+\frac23n)*(\frac13m^3+m^2+\frac23m)}4$$

—–>

$$\frac{(n^3+3n^2+2n)*(m^3+3m^2+2m)}{36}$$

O(1)解决

codes:

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const long long maxn=1e9+7;
long long n,m;

long long get(long long x,long long y)
{
    long long ans=1;
    while(y)
    {
        if(y&1)
        {
            ans=(ans*x)%maxn;
        }
        x=(x*x)%maxn;
        y>>=1;
    }
    return ans;
}

long long Get(long long x)
{
    return((x*x)%maxn*x%maxn+3*x*x%maxn+2*x%maxn)%maxn;
}

int main()
{
    freopen("loop.in","r",stdin);
    freopen("loop.out","w",stdout);
    scanf("%lld%lld",&n,&m);
    long long x=Get(n%maxn),y=Get(m%maxn);
    printf("%lld",(x*y%maxn)*get(36,maxn-2)%maxn);
}