排列组合

先补充一个其它的东西—–二项式定理：

排列：

定义：从n个不同元素中，任取m(m≤n,m与n均为自然数,下同）个元素按照一定的顺序排成一列，叫做从n个不同元素中取出m个元素的一个排列；从n个不同元素中取出m(m≤n）个元素的所有排列的个数，叫做从n个不同元素中取出m个元素的排列数，用符号 A(n,m）表示。
计算公式：
此外规定0！==1

其满足：A(n,m)=n×A(n-1,m-1);(编程时可用此递推)

组合：

定义：从n个不同元素中，任取m(m≤n）个元素并成一组，叫做从n个不同元素中取出m个元素的一个组合；从n个不同元素中取出m(m≤n）个元素的所有组合的个数，叫做从n个不同元素中取出m个元素的组合数。用符号 C(n,m) 表示。
计算公式：
两个性质：
1.C(n,m)=C(n,n-m)
2.C(n,m)=C(n-1,m)+C(n-1,m-1)；(编程时可用此递推)

对于性质1，
当m>n/2时，就可转换成后一个式子，使运算简化。[当C(n,x)==C(n,y)时，要么x==y要么x+y==n]

对于性质2，
可这样理解:c(n,m)即为从n件物品中选m件的方案数。如果第m件物品不选，方案数就变为c(n-1,m)，如果选第m件物品，方案数就变为c(n-1,m-1)，总方案数就为两种情况的方案数之和。

最后对于排列组合，补充一个几何问题：
某区有7条南北向街道，5条东西向街道（如图）

⑴图中共有多少个矩形？
⑵从A点到B点最近的走法有多少种？

分析：⑴在7条竖线中任选2条，5条横线中任选2条，这样4条线
可组成1个矩形，故可组成矩形C（7,2）·C（5,2）=210个
⑵每条东西向的街道被分成6段，每条南北向的街道被分成4段，从A到B最短的走法，无论怎样走，一定包括10步，其中6步方向相同，另外4步方向相同，每种走法，即是从10步中选出6步，这6步是走东西方向的，共有C（10,6）=210种走法（同样可以从10段中选出4段走南北方向，每一种选法即是1种走法）。所以共有210种走法。（鉴于某个xbq(dalao),我不得不将“段”字改为“步”字，它说这样生动…）

来到题：
bzoj3505
题意：
给定一个nxm的网格，请计算三点都在格点上的三角形共有多少个。注意三角形的三点不能共线。
输入一行，包含两个空格分隔的正整数m和n。
输出一个正整数，为所求三角形数量。

题解：
首先答案就是所有取出三个点的方案数减去会三点共线的方案数。
显然n*m的网格上有(n+1)(m+1)个整点，然后令t=(n+1)(m+1)，那么取三个点的方案数就是C(t,3)。
接下来考虑怎么算三点共线的方案数：

有一个结论是在(a,b) (x,y)两点构成的线段上有gcd(a-x,b-y)-1个整点(a>x,b>y)

我们固定(a,b) (x,y)为共线的三点中左边两个，那么第三个点的方案数就是gcd(a-x,b-y)-1
但是这样枚举abxy的复杂度是O(n^2*m^2)
优化是把这线段平移到原点处，平移，平移，平移啊，那么会发现其实只要枚举(0,0) (x-a+1,y-b+1)，其他的线段平移就可以了。
画完很容易发现这样(0,0) (x-a+1,y-b+1)的线段可以平移出(n-i+1)*(m-j+1)种不同方案。

so 代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9'){ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
inline int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
int n,m;
ll c[1005005][4];
ll ans,tmp;
void getc()
{
     c[0][0]=1;
     for(int i=1;i<=n*m;i++)
     {
         c[i][0]=1;
         for(int j=1;j<=3;j++)
             c[i][j]=c[i-1][j-1]+c[i-1][j];
         }
}
void solve()
{
     ans=c[n*m][3];
     ans-=n*c[m][3];
     ans-=m*c[n][3];
     for(int i=1;i<n;i++)
         for(int j=1;j<m;j++)
         {
             tmp=gcd(i,j)+1;
             if(tmp>2)
                 ans-=(tmp-2)*2*(n-i)*(m-j);
         }
}
int main()
{
    n=read();m=read();
    n++;m++;
    getc();
    solve();
    printf("%lld",ans);
    return 0;
}