hdu5116(dp + 计数）

题目链接

解题报告：
这是一个计数的问题，分析一下后可以发现，答案就等于合法的集合数的平方减去有相交的组合的个数。
可以先处理出几个关键数组
$dp[i][j] = 一条边为i，另一条边的范围是1～j，有多少个合法的L形 N$
$cnt[i][j] = 以(i,j)为拐点的合法的L型的数量 N$
$T[i][j] = 经过(i,j)的合法的L型的数量 N$
而我们分析一下有相交的组合的个数，可以分为三类。
1、两个L型的拐点都是（i，j），则这一类的数量为cnt ＊ cnt
2、一个拐点为（i，j），则这一类的数量为cnt ＊ （T － cnt） ＊ 2
3、不交于拐点，则可以枚举拐点为（i，j）的L型右边的长度，这样有dp[l][U[i][j]]个，然后可以与经过（i ＋ 1，j）到（i ＋ l，j）的L型交。（具体见代码）

//
//  main.cpp
//  hdu5116
//
//  Created by 黄宇凡 on 9/1/16.
//  Copyright © 2016 黄宇凡. All rights reserved.
//

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;

const int maxn = 205;
typedef long long ll;
int cnt[maxn][maxn];//cnt[i][j] 以(i,j)为拐点的L形合法数量
ll T[maxn][maxn];//过(i,j)点的L形合法数量
int dp[maxn][maxn];//i固定，j合法的数量
int U[maxn][maxn],R[maxn][maxn];//每格可以向上向右衍生多少

int n;
const int N = 40005;
bool mp[maxn][maxn];
int x[N],y[N];

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 i = 1;i <= 200;i++){
        for(int j = 1;j <= 200;j++){
            dp[i][j] = f(i,j);
            cnt[i][j] = cnt[i - 1][j] + dp[i][j];
        }
    }
}

ll S;

int main(int argc, const char * argv[]) {
    int t;
    init();
    int cas = 0;
    cin >> t;
    while(t--){
        cin >> 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 = 1;i <= n;i++){
            scanf("%d%d",&x[i],&y[i]);
            mp[x[i]][y[i]] = true;
        }
        for(int i = 200;i >= 1;i--){
            for(int j = 200;j >= 1;j--){
                if(mp[i][j]){
                    if(mp[i][j + 1]) R[i][j] = R[i][j + 1] + 1;
                    if(mp[i + 1][j]) U[i][j] = U[i + 1][j] + 1;
                }
            }
        }
        for(int i = 1;i <= 200;i++){
            for(int j = 1;j <= 200;j++){
                if(mp[i][j]){
                    ll A[maxn];
                    memset(A,0,sizeof(A));
                    for(int k = 1;k <= U[i][j];k++){
                        A[k] = dp[k][R[i][j]];
                    }
                    for(int k = U[i][j];k;k--){
                        A[k - 1] += A[k];
                    }
                    S += A[0];
                    for(int k = 0;k <= U[i][j];k++){
                        T[i + k][j] += A[k];
                    }
                }
            }
        }
        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 += p * f * 2;
                    }
                    ans += pp * pp;
                }
            }
        }
        ll ret = S * S - ans;
        printf("Case #%d: ",++cas);
        cout << ret << endl;
    }
    return 0;
}