hdu 5116 dp

/*题意： 给定一个点集合 问从中选出子集 如果存在(a,b) 对于x,y集合Q{(x,y),(x+1,y)...(x+a,y),(x,y+1)...(x,y+b)
并且gcd(a,b)==1 那么集合Q是good 问从给定集合S中找出集合对(A,B) A,B both good A,B都是S的子集  且A交B为空 */

/*
思路：
求出每个点向右，向上延伸的距离。
dp[d][n]: [1, n]中与 d 互质的数有多少个。
T[i][j] : x = i的经过(i, j)的L形个数

还需要再读深入
*/

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define prt(k) cout<<#k" = "<<k<<endl;
typedef long long ll;
const int N = 40020;
const int M = 202;
#include <vector>
bool mp[M+5][M+5];
int dp[M+5][M+5];
int cnt[M+5][M+5];

ll S; // all good sets
//P = {(x, y), (x + 1, y), ... , (x + a, y), (x, y + 1), ... , (x, y + b)}(a, b ≥ 1)
int gcd(int a, int b) { return b==0?a:gcd(b, a%b); }
int U[M][M], R[M][M]; //坐标轴：X轴向上，Y轴向右
ll T[M][M]; // T[i][j] : x = i的经过(i, j)的L形个数

int f(int i, int j)
{
    int ret = 0;
    for (int k=1;k<=j;k++) if (gcd(i,k)==1) ret++;
    return ret;
}

void init()
{
    memset(cnt, 0, sizeof cnt);
    for (int j=1;j<=200;j++)
        for (int i=1;i<=200;i++)
    {
        dp[i][j] = f(i, j);
        cnt[i][j] = cnt[i-1][j] + dp[i][j];//记录1~i中的数和1~j的数互质的数目之和
    }
}

int n;

int main()
{
    freopen("in.txt", "r", stdin);
    init();//预处理
    int re; scanf("%d", &re); int ca = 1;
    while (re--)
    {
        scanf("%d", &n);
        memset(mp, false, sizeof mp);
        memset(U, 0, sizeof U);
        memset(R, 0, sizeof R);
        memset(T, 0, sizeof T);
        S = 0;
        for (int i=0;i<n;i++)
        {
            int x, y; scanf("%d%d", &x, &y);
            mp[x][y] = true;
        }
        n = 200;
        for (int i=200;i>=1;i--)
            for (int j=200;j>=1;j--)
            if (mp[i][j]) {
                if (mp[i+1][j]) U[i][j]=U[i+1][j]+1;//向上延伸的距离
                if (mp[i][j+1]) R[i][j]=R[i][j+1]+1;//向右延伸的距离
            }
        for (int i=1;i<=200;i++)
            for (int j=1;j<=200;j++)
        if (mp[i][j]) {
            ll tmp[M];
            memset(tmp, 0, sizeof tmp);
            for (int k=1; k<=U[i][j]; k++)//a b
                tmp[k] = dp[k][R[i][j]];// [1, R[i][j]]中与 k 互质的数有多少个。
            for (int k = U[i][j]; k; k--)
                tmp[k-1] += tmp[k];
            S += tmp[0];//
            for (int k=0; k<=U[i][j]; k++)
                T[i+k][j] += tmp[k];//x = i的经过(i, j)的L形个数
        }
        ll ans = 0;
        for (int i=1;i<200;i++)
        {
            for (int j=1;j<200;j++)
            if (U[i][j]&&R[i][j]) {
                ll f = T[i][j];
                ll pp, p;
                pp = cnt[U[i][j]][R[i][j]];
                f -= pp;
                for (int k=1;k<=R[i][j];k++)
                {
                    f += T[i][j+k];
                    p = dp[k][U[i][j]];
                    ans += f * p * 2;
                }
                ans += pp * pp;
            }
        }
        ll ret = S * S - ans;
        printf("Case #%d: %I64d\n", ca++ , ret);
    }
    return 0;
}