牛客多校3：Counting Triangles（找规律）

传送门
题目：

namespace GenHelper
{
    unsigned z1,z2,z3,z4,b,u;
    unsigned get()
    {
        b=((z1<<6)^z1)>>13;
        z1=((z1&4294967294U)<<18)^b;
        b=((z2<<2)^z2)>>27;
        z2=((z2&4294967288U)<<2)^b;
        b=((z3<<13)^z3)>>21;
        z3=((z3&4294967280U)<<7)^b;
        b=((z4<<3)^z4)>>12;
        z4=((z4&4294967168U)<<13)^b;
        return (z1^z2^z3^z4);
    }
    bool read() {
      while (!u) u = get();
      bool res = u & 1;
      u >>= 1; return res;
    }
    void srand(int x)
    {
        z1=x;
        z2=(~x)^0x233333333U;
        z3=x^0x1234598766U;
        z4=(~x)+51;
      	u = 0;
    }
}
using namespace GenHelper;
bool edge[8005][8005];
int main() {
  int n, seed;
  cin >> n >> seed;
  srand(seed);
  for (int i = 0; i < n; i++)
    	for (int j = i + 1; j < n; j++)
        	edge[j][i] = edge[i][j] = read();
 	return 0;
}

示例1
输入
10 114514
输出
35
说明
There’re 35 triangles that all three edges have the same color.
大致翻译：

分析：
由已知得到一个无向图，每两点之间必有连线，连线的值为0或1，任意三点必组成一个三角形；
三角形的边值可分为两种情况：
1.三边都相同（为1 1 1 或0 0 0）;
2.两边相同，一边不同（为1 1 0 或 0 0 1）；
由于三重循环必爆，所以我们可以通过计算总三角形的个数 – 不符合三角形的个数即答案；
三角形的总个数计算公式：
sum=n*(n-1)*(n-2)/6;(n>=2时有意义)
下面主要是分析如何求不符合三角形的总个数（即情况2的三角形）：

由题目所给的例子做一个简要的分析：
题目所给例子的邻接矩阵：

举0所连接的三角形的一个例子：
这里我们是通过求异边对数来求解情况2的三角形的个数

由于情况2的三角形为两边相同一边不同，如：

我们在计算时 a 和 c 统计了一次， b 和 c 也统计了一次，所以情况二的三角形的个数应为（1的边数）*（0的边数）/ 2才为情况二所统计三角形的总个数；
代码：

#include<bits/stdc++.h>
using namespace std;
namespace GenHelper
{
unsigned z1,z2,z3,z4,b,u;
unsigned get()
{
    b=((z1<<6)^z1)>>13;
    z1=((z1&4294967294U)<<18)^b;
    b=((z2<<2)^z2)>>27;
    z2=((z2&4294967288U)<<2)^b;
    b=((z3<<13)^z3)>>21;
    z3=((z3&4294967280U)<<7)^b;
    b=((z4<<3)^z4)>>12;
    z4=((z4&4294967168U)<<13)^b;
    return (z1^z2^z3^z4);
}
bool read()
{
    while (!u) u = get();
    bool res = u & 1;
    u >>= 1;
    return res;
}
void srand(int x)
{
    z1=x;
    z2=(~x)^0x233333333U;
    z3=x^0x1234598766U;
    z4=(~x)+51;
    u = 0;
}
}
using namespace GenHelper;
bool edge[8005][8005];
int main() {
  int n, seed;
  cin >> n >> seed;
  srand(seed);
  for (int i = 0; i < n; i++)
    	for (int j = i + 1; j < n; j++)
        	edge[j][i] = edge[i][j] = read();
    if(n<3){
        cout<<0;
        return 0;
    }
    long long sum=0,s1,s0;
    for(int i=0;i<n;i++){
        s1=0,s0=0;
        for(int j=0;j<n;j++){
            if(i==j)continue;
            if(edge[i][j])s1++;//统计该点连接的黑色边数
            else s0++;//统计该点连接的白色边数
        }
        sum+=s1*s0;
    }
    cout<<1ll*n*(n-1)*(n-2)/6-sum/2;
 	return 0;
}