Codeforces 599 D Spongebob and Squares

传送门：
http://codeforces.com/problemset/problem/599/D

题意：
求出所有的n和m，使得在n*m的矩形里恰有k个不同的正方形

和明显的数学公式题嘛！！
n*m的矩形中的不同正方形的数量必定可以用公式去表示出来，必定是枚举n，然后可以在O（1）的时间内把m解出来，那么n要枚举到多大呢，1e18，估计是n3次方，那么就枚举到1e6吧 ，且最后答案一定是对称的，故存完之后再排个序貌似也是可以的！！！

是这个样子的，我们可以假设n<=m，那么此时由下面那个公式可得 ，n不会超过2e6，与我之前的猜测吻合，那么此时我只需要设置一个上限，然后判断当n>m的时候终止， 最后再去判断一下枚举的一半的临界值是不是一个正方形，不用排序，直接顺着反着输出就ok了！！！！

f(n,m) = sigma(k=1,k=min(n,m))(n-k+1)*(m-k+1)=x的解有多少个
化简之后，我们可以得到f(n,m) = n^2m+n^2+n*m+n-(n+1)*n/2*(n+m+2)+n*(n+1)*(2n+1)/6
k=6*n*n*m - 3*(n+m)*n*(n-1) + n*(n-1)*(2*n-1)。

code:

#include<bits/stdc++.h>
#define LL long long
using namespace std;
struct Node{
    LL x, y;
};
vector<Node> ans;
LL limit = 1000000*6;
void getans(LL x)
{
    for(LL i = 1; i <= limit; i++)
    {
        LL a = 6*x + i*i*i - i;
        LL b = 3*i*i + 3*i;
        if(a % b)
            continue;
        LL m = a / b;
        if(i > m)
            break;
        Node now;
        now.x = i;
        now.y = m;
        ans.push_back(now);
    }
}
int main()
{
    LL x; cin>>x; getans(x);
    int top = ans.size()-1;
    if(ans[top].x == ans[top].y)
        top = ans.size() * 2 - 1;
    else
        top = ans.size() * 2;
    printf("%d\n", top);
    int i;
    top = ans.size()-1;
    for(i = 0; i <= top; i++)
        printf("%lld %lld\n", ans[i].x, ans[i].y);
    if(ans[top].x == ans[top].y)
        i = top-1;
    else
        i = top;
    for(; i >= 0; i--)
        printf("%lld %lld\n", ans[i].y, ans[i].x);
    return 0;
}