Language: C++

#pragma region include
#include <iostream>
#include <iomanip>
#include <stdio.h>

#include <sstream>
#include <algorithm>
#include <cmath>
#include <complex>

#include <string>
#include <cstring>
#include <vector>
#include <tuple>
#include <bitset>

#include <queue>
#include <complex>
#include <set>
#include <map>
#include <stack>
#include <list>

#include <fstream>
#include <random>
//#include <time.h>
#include <ctime>
#pragma endregion //#include
/////////
#define REP(i, x, n) for(int i = x; i < n; ++i)
#define rep(i,n) REP(i,0,n)
/////////
#pragma region typedef
typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
#pragma endregion //typedef
////定数
const int INF = (int)1e9;
const LL MOD = (LL)1e9+7;
const LL LINF = (LL)1e18;
const double PI = acos(-1.0);
const double EPS = 1e-9;
/////////
using namespace::std;

/*
thx
http://phyllo-algo.hatenablog.com/entry/2017/01/29/050622
*/
const int MAX = 100005;
map<int,int> query[MAX];
void solve(){
	int N,M;
	cin >> N >> M;
	vector< vector<int> > pts(N,vector<int>(3));
	for(int i=0;i<N;++i){//N種類の区間,集まる種類がある
		cin >> pts[i][0] >> pts[i][1];//共に[1,M]
		pts[i][2] = 0;//種類
	}

	//d=1ではすべての区間を通る制約になっている=N
	vector<int> temp(3);//追加するクエリ用
	for(int d=1;d<=M;d++){
		for(int p=1;(LL)d*(LL)p<=M;p++){
			//Xを含まず、Yを含む区間では、
			//d毎に進んだ場合にその区間を通過する一番左になる。
			int X = d*(p-1);
			int Y = d*p;
			
			if( query[X].count(Y-1) == 0 ){
				query[X][Y-1] = 1;
				temp[0] = X;
				temp[1] = Y-1;
				temp[2] = 1;//種類
				pts.push_back( temp );
			}
			if( query[X].count(MAX) == 0 ){
				query[X][MAX] = 1;
				temp[0] = X;
				temp[1] = MAX;
				temp[2] = 1;
				pts.push_back( temp );
			}
			if( query[Y].count(Y-1) == 0 ){
				temp[0] = Y;
				temp[1] = Y-1;
				temp[2] = 1;
				pts.push_back( temp );
			}
			if( query[Y].count(MAX) == 0 ){
				temp[0] = Y;
				temp[1] = MAX;
				temp[2] = 1;
				pts.push_back( temp );
			}
		}
	}
	sort( pts.begin(), pts.end() );

	BITree bit(M+1);
	int size = pts.size();
	for(int i=0;i<size;++i){
		if( pts[i][2] == 0 ){
			bit.add( pts[i][1]+1, 1 );
		}else{
			int res = bit.sum( pts[i][1]+1 );
			query[ pts[i][0] ][ pts[i][1] ] = res; 
		}
	}

	//
	for(int d=1;d<=M;d++){
		int ret = 0;
		for(int p=1;(LL)d*p<=M;p++){
			int X = d*(p-1);
			int Y = d*p;
			ret += query[Y][MAX] - query[Y][Y-1] - query[X][MAX] + query[X][Y-1];
			/*
			区間[A,B]をプロットしている。
			X < A
			A <= Y <= B
			から
			X < A <= Y
			Y <= B
			この領域に含まれる点(区間)の個数を数える。
			query[_1][_2]は(_1,_2)と(0,0)で作る長方形内の個数にしてるので。
			query[Y][MAX]	:=(Y,MAX)の長方形
			query[Y][Y-1]	:=(Y,Y-1)の長方形 ( y<Yの長方形
			query[X][MAX]	:=(X,MAX)の長方形 (x<=Xの部分
			query[X][Y-1]	:=(X,Y-1)の長方形 (x<=Xかつ y<Yの分、引きすぎたのでたす
			これらの値を点を順番にプロットする時に、計算している。
			*/
		}
		cout << ret << endl;
	}
}

#pragma region main
signed main(void){
	std::cin.tie(0);
	std::ios::sync_with_stdio(false);
	std::cout << std::fixed;//小数を10進数表示
	cout << setprecision(16);//小数点以下の桁数を指定//coutとcerrで別	

	solve();
}
#pragma endregion //main()


Compiler output:
a.cc: In function 'void solve()':
a.cc:101:9: error: 'BITree' was not declared in this scope
  101 |         BITree bit(M+1);
      |         ^~~~~~
a.cc:105:25: error: 'bit' was not declared in this scope
  105 |                         bit.add( pts[i][1]+1, 1 );
      |                         ^~~
a.cc:107:35: error: 'bit' was not declared in this scope
  107 |                         int res = bit.sum( pts[i][1]+1 );
      |                                   ^~~