Highways UVA - 1393 容斥原理

问题

分析

参考：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/82942861
https://www.cnblogs.com/wafish/p/10465193.html
只考虑‘/’形状的线条，然后乘2
点[i,j]互素，那么和(0,0)之间的连线不经过任何点，所以用dp[i,j]表示在矩形[i,j]的互素点的个数，也就是(0,0)点连接的直线数量
由于对称性，点(0,0)和点(i,j)上的直线数目一样，所以可以用dp(i,j)来表示点(i,j)连接的不经过其他点的线段数量。如3*3的点内部，dp[1,1]=1,dp[1,2]=dp[2,1]=2,dp[3,3]=3
ans[i][j]是[i][j]矩形内的所有的‘/’方向的直线数量，ans[i][j]=ans[i-1][j]+ans[i][j-1]-ans[i-1][j-1]+dp[i][j]-dp[i/2][j/2]，前面都好理解，最后减去的dp[i/2][j/2]是减去由[i/2,j/2]延伸出来的的线段，这样就可以消除重复的部分

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <vector>
#include <utility>
using namespace std;
typedef long long ll;
const int maxn=300+1;
int dp[maxn][maxn],ans[maxn][maxn];

int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}

int main(void){
    for(int i=1;i<maxn;++i){
        for(int j=1;j<maxn;++j){
            dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+(gcd(i,j)==1?1:0);
        }
    }
    for(int i=1;i<maxn;++i){
        for(int j=1;j<maxn;++j){
            ans[i][j]=ans[i-1][j]+ans[i][j-1]-ans[i-1][j-1]+dp[i][j]-dp[i/2][j/2];
        }
    }
    int n,m;
    while(cin>>n>>m && n){
        printf("%d\n",ans[n-1][m-1]<<1);  //结果要乘以2
    }
}