2021牛客暑假多校训练#3 J.Counting Triangles

//编译器级优化
// #pragma GCC optimize("O2")
// #include <bits/stdc++.h>
#include <string>
#include <iostream>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <time.h>
#include <cstring>
#include <cmath>
#include <set>
#include <unordered_map>
#include <algorithm>
#include <cstdio>
#include <cassert>
#include <bitset>
// #include <bits/extc++.h>
//外部包含文件
// #include "D:/program_learning/vs_code/workplace/ACMtemplate/data_structure/SparseTable.hpp"
// #include "D:/program_learning/vs_code/workplace/ACMtemplate/computational_geometry/Convex.hpp"
//类型别名
using namespace std;
// using B_ll = __int128_t;
// using B_ull = __uint128_t;
using f32 = float;
using f64 = double;
using ll = long long;
//排序用
#define htl(v, n) (v),(v)+(n), greater<decltype((v)[1])>()
#define lth(v, n) (v),(v)+(n), less<decltype((v)[1])>()
//宏函数
#define Mt(args...) make_tuple(args)
#define Mp(a, b) make_pair((a), (b))
#define Gt(a, b) std::get<(b)>((a))
#define Add(a, b) accumulate((a), (b), 0)
#define Abs(a)  (((a) > 0) ? (a) : -(a))
#define Max(a, b) (((a) > (b)) ? (a) : (b))
#define Min(a, b) (((a) < (b)) ? (a) : (b))
#define Ceil(a, b) (((a) + (b) - 1) / (b))
#define Floor(a, b) ((a) / (b))
#define Pow2(x) ((x) * (x))
#define Pow3(x) ((x) * (x) * (x))
//浮点数关系运算
#define isZero(x) ((x)>0?(x):-(x)<eps)
#define equal(x, y) (abs((x)-(y))<eps)
#define lthan(x, y) ((y)-(x)>eps)
#define bthan(x, y) ((x)-(y)>eps)
//简化循环
#define FOR(i, a, b) for(int i = (ll)(a); i <= (ll)(b); i++)
#define rFOR(i, a, b) for(int i = (ll)(a); i >= (ll)(b); i--)
#define gFOR(i, a) for(int i = h[(a)] ; ~i ; i = ne[i])
#define bFOR(i, a, b) for(int i = (a); i <= (b); i <<= 1)
#define jFOR(i, a, b, step) for(int i = (a); i <= (b); i += (step))
#define iFOR(it, begin, end) for(auto it = (begin); it!=(end);it++)
#define aFOR(a, s) for(auto (a) : (s))
//debug用
#define _var(x) cout << #x << " = " << (x) << endl
#define _arr1(a) for(ll i = 0 ; i < (a).size() ; i ++)  cout << (a)[i] << ((i == (a).size() - 1) ? endl : space)
#define _arr2(a) for(auto s : (a) )for(ll j = 0 ; j < s.size() ; j ++)   cout << s[j] << ((j == s.size() - 1) ? endl : space)
#define checktime() cerr << "Time : " << 1000 * ((double)clock()) / (double)CLOCKS_PER_SEC << "ms\n" ; system("pause")
//加快io
#define endl "\n"
#define __faster_io__  std::ios::sync_with_stdio(false),std::cin.tie(0),std::cout.tie(0);
//操作输出
#define space " "
#define setd(n) fixed << setprecision(n)
//初始化设置
#define case(t,i) ll t;cin >> t;for(ll i=1;i<=t;i++)
#define _set(b, c) memset((b), (c), sizeof(b))
//模板函数
template<typename T>
int read(T& result) { // fast read
    char ch;int w = 1;
    result = 0;
    
    ch = getchar();
    while(ch != '-' && !(ch >= '0' && ch <= '9'))
        ch = getchar();

    if(ch == '-')
        w = -1;

    while(isdigit(ch)){
        result = result * 10 + (ch - '0');
        ch = getchar();
    }

    result *= w;
    return 1;
}
template<typename T>
void write (T x) { // fast output
    string tmp;string tail = (x < 0) ? "-" : "";
    while(x){
        tmp += x % 10 + '0';
        x /= 10;
    }
    tmp += tail;
    reverse(tmp.begin(), tmp.end());
    cout << tmp;
}
template<typename T>
T qpow(T x, T n, const T& mod){
    ll ans = 1;
    while(n){
        if(n & 1)
            ans = ans % mod * x % mod;
        x = x % mod * x % mod;
        n >>= 1;
    }
    return ans;
}
template<typename T>
inline T lowbit(T& x){
    return x&-x;
}
template<typename T>
T m_plus(T a, T b, const T& mod){
    return (a % mod + b % mod) % mod;
}
template<typename T>
T m_minus(T a, T b, const T& mod){
    return ((a % mod - b % mod) % mod + mod) % mod;
}
template<typename T>
T m_mult(T a, T b, const T& mod){
    return a % mod * b % mod;
}
template<typename T>
T Mod(T x, T y)
{
    return (x % y + y) % y;
}
//常数
constexpr f32 fINF32 = 1e30;
constexpr f64 fINF64 = 1e300;
constexpr int INF32 = 0x3f3f3f3f;
constexpr ll INF64 = 0x3f3f3f3f3f3f3f3f;
constexpr f32 _fINF32 = -1e30;
constexpr f64 _fINF64 = -1e300;
constexpr int _INF32 = 0xcfcfcfcf;
constexpr ll _INF64 = 0xcfcfcfcfcfcfcfcf;
constexpr f64 eps = 1e-6;
const f64 pi = acos(-1);
constexpr ll mod = 9 + 7;
constexpr ll Bit = 1 << 20;
constexpr ll N = 1e1 + 5;
/*定义全局变量*/
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;
    }
}
ll n, seed; // 因为seed范围最大是1e9，如果用int就会WA，只能用long long。
bool edge1[8005][8005];
/*辅助函数*/
void solve(){
    cin >> n >> seed;
    GenHelper::srand(seed);
    FOR(i, 0, n - 1)FOR(j, i + 1, n - 1){ // 根据给定的随机数种子建图
        bool tmp = GenHelper::read();
        edge1[i][j] = edge[j][i] = tmp; // true表示black，false表示white
    }
    ll all = n * (n - 1) * (n - 2) / 6; // 一个完全图中所有能构成的三角形数目
    ll sum = 0; // 不满足要求的三角形的数目
    FOR(i, 0, n - 1){
        ll w = 0, b = 0; // w表示从某个点连出来的边是白边的数目，b表示从某个点连出来的边是黑边的数目。
        FOR(j, 0, n - 1){ // 考虑任何一个不满足要求的三角形的边的构成，要么是黑黑白形，要么是白白黑形。
                          // 我们对于每一个点，枚举一下所有连出去的黑边和白边的数目，其中任意一条黑边搭配上任意一条白边都一定是不满足要求的三角形上的一对边。
                          // 又因为我们是枚举的异色边，而一个不满足要求的三角形上有两对异色边和一对同色边，即，对于每个不满足要求的三角形它对我们计数结果的贡献是2。
                          // 所以最后将我们计数出来的异色边对的数目除以2就是所有不满足要求的三角形的个数，用总三角形数减去即可。
            if(i == j)  continue; // 图中不存在自环，所以枚举到相同的点要跳过。
            if(edge1[i][j] == true) b ++;
            else w ++;
        }
        sum += w * b;
    }
    cout << res - (sum / 2) << endl;
}
/*程序入口*/
int main()
{
#ifdef  LOCAL
    freopen("D:/program_learning/vs_code/workplace/ACMtemplate/DataIn.txt", "r", stdin);
    freopen("D:/program_learning/vs_code/workplace/ACMtemplate/DataOut.txt", "w", stdout);
#endif
	
    __faster_io__

    case(t, i){
        solve();
    }

#ifdef  LOCAL
    cout.flush();
    freopen("CON", "r", stdin);
    freopen("CON", "w", stdout);
#endif

    checktime();
    return 0;
}