这个题目的难点在于(自己设计中间变量对结果进行统计)对结果的统计:
思路大概如此:
首先统计出所有的符合条件的L形有多少个记为 x,然后统计出集合相交的A,B有多少个记为y,那么结果就是x*x - y;
对于所有的L形的统计比较简单,递推计算出每个点可以往上往下延伸多长,然后计数统计。
定义cross(i,j) 为记过点i,j的所有L形,cnt[ i ][ j ] 为x 从1 -> i 并且 y 从 1 - >j 中有多少满足条件的L形, c[ i ][ j ]为x = i并且 y 从 1 - >j 中有多少满足条件的L形, u[ i ][ j ] 为向上延伸的
长度,r[i][j] 为向右延伸的长度。
对于重复的集合,要这样计算,分别统计经过一个特定转角i,j的所有重复集合对,
1 , 那么转角同为 点 i,j , 那么 总数为cnt[ i ][ j] * cnt[ i ][j];
2 , 经过点i,j 但转角不是点i,j与转交是i,j的之间的重复问题, 总数为cnt[ i ][ j ] * ( cross[i][j] - cnt[ i ][ j ]) *2
3 , 不经过点i,j 与转角是i,j点所有重复情况,那么就要沿着y轴,从i,j出发想右走,一一统计,具体见代码。
//#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <map>
#include <set>
#include <vector>
#include <cctype>
#include <cmath>
#include <queue>
#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mem(a,n) memset(a,n,sizeof(a))
#define rep(i,n) for(int i=0;i<(int)n;i++)
#define rep1(i,x,y) for(int i=x;i<=(int)y;i++)
using namespace std;
typedef long long ll;
const ll oo = 1e12;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
const int N = 205;
int gcd(int a,int b){return b==0 ? a:gcd(b,a%b);}
int cnt[N][N],c[N][N];
void init(){
mem(cnt,0);
rep1(i,1,N-1) rep1(j,1,N-1){
int te = 0;
rep1(k,1,j) te+=(gcd(i,k)==1);
cnt[i][j] = cnt[i-1][j]+(c[i][j] = te);
}
}
bool vis[N][N];
int u[N][N],r[N][N],n,cr[N][N];
int main()
{
init();
int T,kase=1;
scanf("%d",&T);
while(T--){
scanf("%d",&n);
mem(vis,false);
rep1(i,1,n) {
int x,y; scanf("%d %d",&x,&y);
vis[x][y] = true;
}
ll ans1 = 0;
mem(u,0); mem(r,0); mem(cr,0);
int tem[N];
for(int i=200;i>=1;i--)
for(int j=200;j>=1;j--)if(vis[i][j]){
u[i][j]=u[i+1][j]+vis[i+1][j];
r[i][j]=r[i][j+1]+vis[i][j+1];
ans1+=cnt[u[i][j]][r[i][j]];
mem(tem, 0);
for(int k=u[i][j];k>=1;k--)
tem[k] = c[k][r[i][j]];
for(int k=u[i][j];k>=1;k--)
tem[k-1]+=tem[k];
for(int k=0;k<=u[i][j];k++)
cr[i+k][j] += tem[k];
}
ll ans2 = 0;
for(int i=200;i>=1;i--)
for(int j=200;j>=1;j--)if(vis[i][j]){
ll te = cnt[u[i][j]][r[i][j]];
ll yu = 0;
for(int k = 1; k<=r[i][j];k++){
yu += cr[i][j+k];
ans2 += yu*c[k][u[i][j]]*2;
}
ans2+=te*te + (ll)(cr[i][j]- te)*te*2;
}
printf("Case #%d: %I64d\n",kase++,ans1*ans1-ans2);
}
return 0;
}