『树状数组』绕圈跑

显然我们可以预处理出每一个人在事件结束后跑步的圈数s_i

那么答案可以表示为：$$\sum _{i=1}^{n-1}\sum _{j=i+1}^n\lfloor s_j-s_i\rfloor$$

我们来思考一下$\lfloor s_1-s_2\rfloor$和$\lfloor s_1\rfloor -\lfloor s_2\rfloor$的区别。

显然，如果$s_1$的小数部分大于$s_2$的小数部分，那么答案没有影响。反之，答案便会差一个$1$.

因此，对于$i<j$,且 { s i } &gt; { s j } \{s_i\}&gt;\{s_j\} {si​}>{sj​}来说，答案就会差了$1$.此时，这个减去的个数就是小数部分的逆序对个数。

这个怎么求呢？我们只需要对小数离散化就可以。

• 注意离散化需要使用排序的方法来解决，且判断相等的精度要设的更精确。

代码如下：

#include <bits/stdc++.h>

#define int long long

using namespace std;
const int N = 200000;

int n, L, c, ans = 0, cnt = 0;
int a[N];
double b[N];

struct node {
	int num, id;
	double v;
	friend bool operator < (node p1,node p2) {
		return p1.v < p2.v;
	}
} t[N];

map<double,int>mp;

inline int read(void)
{
	int s = 0, w = 1;char c = getchar();
	while (c<'0' || c>'9') {if (c == '-') w = -1; c = getchar();}
	while (c>='0' && c<='9') s = s*10+c-48,c = getchar();
	return s*w;
}

struct TREE {
	int S[10000000] = {};
	#define lowbit(i) (i & -i)
	void add(int x,int v)
	{
		for (int i=x;i<=n;i+=lowbit(i))
		    S[i] += v;
		return;
	}
	int ask(int x)
	{
		int sum = 0;
		for (int i=x;i>=1;i-=lowbit(i))
		    sum += S[i];
		return sum;
	}
} tree;

void Read_init(void)
{
	n = read(), L = read(), c = read();
	for (int i=1;i<=n;++i) 
		a[i] = read();
	sort(a+1,a+n+1);
	double Time = 1.0*c*L/a[n];
	for (int i=1;i<=n;++i)
	    b[i] = 1.0*a[i]*Time/c;
	return;
}

bool cmp(node p1,node p2) {
	return p1.id < p2.id;
}

void Work_Lisan(void)
{
	for (int i=1;i<=n;++i)
	    t[i].v = b[i]-trunc(b[i]), t[i].id = i;
	sort(t+1,t+n+1);
	for (int i=1;i<=n;++i)
	{
	    if (i == 1 || fabs(t[i].v-t[i-1].v) > 0.0000001) 
	    	t[i].num = ++cnt;
	    else t[i].num = t[i-1].num;
	}
	sort(t+1,t+n+1,cmp);
	return;
}

void Get_answer(void)
{
	for (int i=1;i<=n;++i) {
		ans += (i-1)*trunc(b[i]);
		ans -= (n-i)*trunc(b[i]);
	}
	for (int i=n;i;--i)
	    ans -= tree.ask(t[i].num-1), 
	    tree.add(t[i].num,1);
	printf("%lld\n", ans);
	return;
}

signed main(void)
{
	freopen("running.in","r",stdin);
	freopen("running.out","w",stdout);
	Read_init();
	Work_Lisan();
	Get_answer();
	return 0;
}