UVA 1393 Highways，UVA 12075 Counting Triangles —— （组合数，dp）

　　先看第一题，有n*m个点，求在这些点中，有多少条直线，经过了至少两点，且不是水平的也不是竖直的。

　　分析：由于对称性，我们只要求一个方向的线即可。该题分成两个过程，第一个过程是求出n*m的矩形中，dp[i][j]代表在这个矩形中终点是到(i,j)这个点的满足题意的直线条数，那么，用dp的话就可以得出递推关系：由长和宽分别小1的左右两个矩形中满足题意的线的条数减去他们共有的矩形中满足的线的条数（容斥减去重复部分），之后还要判断从最左上角的点(1,1)到(i,j)是否可以组成一条线，这个条件是gcd(i,j)是否等于1。

　　之后第二个过程就是递推答案了，设ans[i][j]表示在这个矩形中满足题意的条数，那么同样的，可以由上面的容斥来递推，同时，还要加上这个矩形内到(i,j)这个点满足的条数，另外还要减去一半规模大小的到这个点的线的条数，因为如果(i,j)为(6,8)，那么一半规模下，（1,1）到（3,4）这个点的线和到（6,8）这条线是重复的。

　　这样就做完了题目（最后不要忘了乘以2）。具体见代码：

 #include <stdio.h>
 #include <algorithm>
 #include <string.h>
 #include <iostream>
 #include <vector>
 #include <queue>
 using namespace std;
 typedef long long ll;
 const int N = +;

 int dp[N][N];
 int ans[N][N];

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

 void init()
 {
     for(int i=;i<N;i++)
     {
         for(int j=;j<N;j++)
         {
             dp[i][j] = dp[i-][j] + dp[i][j-] - dp[i-][j-] + (gcd(i,j)==);
         }
     }
     for(int i=;i<N;i++)
     {
         for(int j=;j<N;j++)
         {
             ans[i][j] = ans[i-][j] + ans[i][j-] - ans[i-][j-] + dp[i][j] - dp[i>>][j>>];
         }
     }
 }

 int main()
 {
     init();
     int n,m;
     while(scanf("%d%d",&n,&m)==)
     {
         if(n== && m==) break;
         printf("%d\n",*ans[n-][m-]);
     }
 }

　　做出了这题，第二题就是类似的了。先在所有的点中枚举出选3个点的可能性，然后，减去一条水平或者竖直线上重复的，再减去同在一条斜线上重复的即可。相当类似，具体见代码吧：

 #include <stdio.h>
 #include <algorithm>
 #include <string.h>
 #include <iostream>
 #include <vector>
 #include <queue>
 using namespace std;
 typedef long long ll;
 const int N = +;

 ll dp[N][N];
 ll ans[N][N];

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

 void init()
 {
     for(int i=;i<N;i++)
     {
         for(int j=;j<N;j++)
         {
             dp[i][j] = dp[i-][j] + dp[i][j-] - dp[i-][j-] + (ll)(gcd(i,j)-);
         }
     }
     for(int i=;i<N;i++)
     {
         for(int j=;j<N;j++)
         {
             ans[i][j] = ans[i-][j] + ans[i][j-] - ans[i-][j-] + dp[i][j];
         }
     }
 }

 ll C(ll x) {return x*(x-)*(x-)/;}

 int main()
 {
     init();
     int n,m;
     int cnt = ;
     while(scanf("%d%d",&n,&m)==)
     {
         if(n== && m==) break;
         ll Ans = C((n+)*(m+)) - (n+)*C(m+) - (m+)*C(n+);
         Ans -= *ans[n][m];
         cout<<"Case "<<cnt++<<": "<<Ans<<endl;
     }
 }

　　

　　但是，这两题都要注意的地方是，递推dp时三个矩形都是在右下角的（因为向下或者向右平移一个单位的话条数是不变的），这样递推起来的话只要再考虑从（1,1）这个点到（i，j）这个点的情况即可；而递推ans的时候，矩形是偏左上方的，那么，只要再加上整个大矩形内到（i，j）这个点的情况即可。当然，纯属个人理解。