#P7487. 「Stoi2031」兰亭序

    ID: 6542 Type: RemoteJudge 2000ms 128MiB Tried: 0 Accepted: 0 Difficulty: 6 Uploaded By: Tags>数学递推递归O2优化素数判断,质数,筛法

「Stoi2031」兰亭序

题目背景

无关风月 我题序等你回 悬笔一绝 那岸边浪千叠 情字何解 怎落笔都不对 而我独缺 你一生的了解 ——《兰亭序》

题目描述

月非常喜欢复数,尤其喜欢形如 e2πite^{2\pi it} 的复数。她选择了两个正整数 n,kn,k,并将 1+e2πix1xkn1+e^{\frac{2\pi i x_1 \dots x_k}{n}} 称为 (x1,,xk)(x_1,\dots,x_k)绝对度,所有满足 1xin1 \le x_i \le n (i{1,2,,k})(i \in \{1,2,\dots,k\})(x1,,xk)(x_1,\dots,x_k)绝对度 之积称为 (n,k)(n,k)无关度。现在她想请你帮她对 t{1,2,,k}t \in \{1,2,\dots,k\} 求出 (n,t)(n,t)无关度 ansmod335544323ans \bmod{335544323}。由于回答太多的数太麻烦,你只要回答她所有答案进行异或运算后的结果。

输入格式

一行两个正整数 n,kn,k

输出格式

一行一个自然数表示答案。

15 2

201012023

1 3

2

521 6

262795752

6546546546546543 22211

388124125

提示

简述版题意:

给定 n,kn,k,对 1tk1 \le t \le k

$$\prod_{x_1=1}^{n}\prod_{x_2=1}^{n}\dots\prod_{x_t=1}^{n}\left( 1+e^{\frac{2\pi ix_1x_2\dots x_t}{n}} \right) \bmod{335544323} $$

输出所有 kk 个答案的异或和。

其中 eit=cost+isinte^{it}=\cos{t}+i\sin{t} 对所有 tRt \in \mathbb{R} 成立,ii 为虚数单位,满足 i2=1i^2=-1

样例解释:

对于第一组样例,t=1,2t=1,2 时答案分别为 2,351843720888322,35184372088832,取模后为 2,2010120212,201012021,异或和为 201012023201012023

对于第二组样例,t=1,2,3t=1,2,3 时答案均为 22,异或和为 22

限于篇幅,剩下的样例不作解释说明。

数据范围:

本题采用捆绑测试,各个 Subtask 的限制与分值如下:

Subtask No. nn \le kk \le 特殊限制 分值
11 1010 11 77
22 11 1010
33 1010 22
44 101810^{18} 10510^5 nn 为偶数
55 1010 nk730n^k \le 730 1616
66 10910^9 10310^3 1919
77 101810^{18} 10510^5 3737

对于 100%100\% 的数据,1n1018,1k1051 \le n \le 10^{18},1 \le k \le 10^5