动态规划（DP计数）：HDU 5116 Everlasting L

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N=210;
typedef long long LL;
LL map[N][N],st[N],S,M;
LL f[N][N],dp[N][N],sum[N][N];
/*f[i][j]:[1,j]中与i互质的数的个数*/
LL dwn[N][N],rht[N][N];int T,cas,k,x,y;
LL Gcd(LL a,LL b){return b?Gcd(b,a%b):a;}
void Prepare(){
    for(int i=1;i<=200;i++)
        for(int j=1;j<=200;j++){
            f[i][j]=f[i][j-1]+(Gcd(i,j)==1);
            dp[i][j]=dp[i-1][j]+f[i][j];
        }
}
void Init(){S=M=0;
    memset(map,0,sizeof(map));
    memset(sum,0,sizeof(sum));
    memset(dwn,0,sizeof(dwn));
    memset(rht,0,sizeof(rht));
}
int main(){
    Prepare();
    scanf("%d",&T);
    while(T--){
        Init();scanf("%d",&k);
        while(k--){scanf("%d%d",&x,&y);map[x][y]=1;}
        for(int i=200;i>=1;i--)
            for(int j=200;j>=1;j--){
                if(!map[i][j])continue;
                if(map[i+1][j])dwn[i][j]=dwn[i+1][j]+1;
                if(map[i][j+1])rht[i][j]=rht[i][j+1]+1;
            }
        
        for(int i=1;i<=200;i++)
            for(int j=1;j<=200;j++){
                if(!map[i][j])continue;
                memset(st,0,sizeof(st));
                for(int k=1;k<=dwn[i][j];k++)
                    st[k]=f[k][rht[i][j]];
                for(int k=dwn[i][j];k>=1;k--)
                    st[k-1]+=st[k];
                for(int k=0;k<=dwn[i][j];k++)
                    sum[i+k][j]+=st[k];
                S+=st[0];            
            }
        
        for(int i=1;i<=200;i++)
            for(int j=1;j<=200;j++){
                if(!dwn[i][j])continue;
                if(!rht[i][j])continue;
                LL tot=dp[dwn[i][j]][rht[i][j]];
                LL calc=sum[i][j]-tot;
                for(int k=1;k<=rht[i][j];k++){
                    calc+=sum[i][j+k];
                    M+=2*calc*f[k][dwn[i][j]];
                }
                M+=tot*tot;
            }    
        printf("Case #%d: %lld\n",++cas,S*S-M);
    }
    return 0;
}