题目#

[NOIP2015 普及组] 求和#

题目背景#

NOIP2015 普及组 T3

题目描述#

一条狭长的纸带被均匀划分出了 $n$ 个格子，格子编号从 $1$ 到 $n$ 。每个格子上都染了一种颜色 $color_i$ 用 $[1,m]$ 当中的一个整数表示），并且写了一个数字 $number_i$ 。

定义一种特殊的三元组： $(x,y,z)$ ，其中 $x,y,z$ 都代表纸带上格子的编号，这里的三元组要求满足以下两个条件：

$xyz$ 是整数, $x<y<z,y-x=z-y$
$colorx=colorz$

满足上述条件的三元组的分数规定为 $(x+z) \times (number_x+number_z)$ 。整个纸带的分数规定为所有满足条件的三元组的分数的和。这个分数可能会很大，你只要输出整个纸带的分数除以 $10,007$ 所得的余数即可。

（ $n, m <= 100000$ ）

题解#

首先一眼看出结果和y没有半毛钱关系，y唯一限制的地方是x，z的奇偶性必须相同

所以我们可以得出要想统计进入答案的格子之间必须同奇偶且同颜色，我们把有相同奇偶性和颜色的格子分成一组。

我们考虑其中的一组，考虑差值法，看一下每加入一个新的数会对答案产生多少的贡献

画出图表容易知道

number\编号	1	2	3	4
n_1	1			1
n_2		1		1
n_3			1	1
n_4	1	1	1	1+1+1

例如加入了该组的第四个，产生的贡献如下然后前缀和统计，计算即可

代码#

1
#include<bits/stdc++.h>
2
#define int long long
3
using namespace std;
4

5
int n, m;
6

7
const int N = 100100;
8
const int mod = 10007;
9

10
int counta[N][2], sum_n[N][2], sum_f[N][2], sum_fn[N][2];//count i,j 表示颜色为i，编号为j的数量，sum i,j 表示颜色为i，编号为j的
11
int a[N], b[N];
12

13
int read()
14
{
15
  int f = 1, x = 0;
16
  char ch = getchar();
17
  while(ch < '0' || ch > '9')
18
  {
19
    if(ch == '-') f = -1;
20
    ch = getchar();
21
  }
22

23
  while(ch >= '0' && ch <= '9')
24
  {
25
    x = x * 10 + ch - '0';
26
    ch = getchar();
27
  }
28
  return f * x;
29
}
30

31
signed main()
32
{
33
  n = read(), m = read();
34
  long long ans = 0;
35
  for(int i = 1;i <= n;i++)
36
  {
37
    a[i] = read();//数字——n
38
  }
39

40
  for(int i = 1;i <= n;i++)
41
  {
42
    b[i] = read();//颜色——b
43
  }
44

45
  for(int i = 1;i <= n;i++)
46
  {
47
    ans = (ans + sum_fn[b[i]][i % 2]) % mod;
48
    counta[b[i]][i % 2] ++;
49
    sum_f[b[i]][i % 2] = (sum_f[b[i]][i % 2] + i) % mod;
50
    sum_n[b[i]][i % 2] = (sum_n[b[i]][i % 2] + a[i]) % mod;
51
    sum_fn[b[i]][i % 2] = (sum_fn[b[i]][i % 2] + (i * a[i]) % mod) % mod;
52
    ans = (ans + (i * sum_n[b[i]][i % 2]) % mod) % mod;
53
    ans = (ans + (a[i] * sum_f[b[i]][i % 2]) % mod) % mod;
54
    ans = (ans + ((  (counta[b[i]][i % 2] - 3) * i * a[i]) % mod + mod ) % mod) % mod;
55
  }
56
  printf("%lld", ((ans % mod) + mod) % mod);
57
  return 0;
58
}