New Year and the Factorisation Collaboration

题意翻译

给定 $n$（$21\leq n \leq 2^{1024}$），求 $n$ 的质因数。保证 $n$ 的所有质因数均为 $4x+3$ （$x$ 为自然数）的形式且至多由 $10$ 个质因数相乘得到。 IO 交互提供了高精度计算。 `+ x y` 给定 $x, y$（$0\leq x < n$），返回 $(x+y)\bmod n$。 `- x y` 给定 $x, y$（$0\leq x < n$），返回 $(x-y)\bmod n$。 `* x y` 给定 $x, y$（$0\leq x < n$），返回 $(x\cdot y)\bmod n$。 `/ x y` 给定 $x, y$（$0\leq x < n$），返回 $(x\cdot y^{-1})\bmod n$。假如 $y, n$ 不互质，返回 `-1`。 `sqrt x` 给定 $x$（$0\leq x < n$），返回 $y$（$0\leq y < n$）满足 $y^2\bmod n=x$，若不存在返回 `-1`。 `^ x y` 给定 $x, y$（$0\leq x < n$），返回 $(x^y)\bmod n$。 你至多可以使用 100 次交互库。运算花费时间表如下： `+ x y` 至多 1 ms。 `- x y` — 至多 1 ms。 `* x y` — 至多 1 ms。 `/ x y` — 至多 350 ms。 `sqrt x` — 至多 80 ms。 `^ x y` — 至多 350 ms。

题目描述

Integer factorisation is hard. The RSA Factoring Challenge offered $ $100\,000 for factoring RSA- $ 1024 $ , a $ 1024 $ -bit long product of two prime numbers. To this date, nobody was able to claim the prize. We want you to factorise a $ 1024 $ -bit number. Since your programming language of choice might not offer facilities for handling large integers, we will provide you with a very simple calculator. To use this calculator, you can print queries on the standard output and retrieve the results from the standard input. The operations are as follows: - `+ x y` where $ x $ and $ y $ are integers between $ 0 $ and $ n-1 $ . Returns $ (x+y) \bmod n $ . - `- x y` where $ x $ and $ y $ are integers between $ 0 $ and $ n-1 $ . Returns $ (x-y) \bmod n $ . - `* x y` where $ x $ and $ y $ are integers between $ 0 $ and $ n-1 $ . Returns $ (x \cdot y) \bmod n $ . - `/ x y` where $ x $ and $ y $ are integers between $ 0 $ and $ n-1 $ and $ y $ is coprime with $ n $ . Returns $ (x \cdot y^{-1}) \bmod n $ where $ y^{-1} $ is multiplicative inverse of $ y $ modulo $ n $ . If $ y $ is not coprime with $ n $ , then $ -1 $ is returned instead. - `sqrt x` where $ x $ is integer between $ 0 $ and $ n-1 $ coprime with $ n $ . Returns $ y $ such that $ y^2 \bmod n = x $ . If there are multiple such integers, only one of them is returned. If there are none, $ -1 $ is returned instead. - `^ x y` where $ x $ and $ y $ are integers between $ 0 $ and $ n-1 $ . Returns $ {x^y \bmod n} $ . Find the factorisation of $ n $ that is a product of between $ 2 $ and $ 10 $ distinct prime numbers, all of form $ 4x + 3 $ for some integer $ x $ . Because of technical issues, we restrict number of requests to $ 100 $ .

输入输出格式

输入格式

The only line contains a single integer $ n $ ( $ 21 \leq n \leq 2^{1024} $ ). It is guaranteed that $ n $ is a product of between $ 2 $ and $ 10 $ distinct prime numbers, all of form $ 4x + 3 $ for some integer $ x $ .

输出格式

You can print as many queries as you wish, adhering to the time limit (see the Interaction section for more details). When you think you know the answer, output a single line of form ! k p\_1 p\_2 ... p\_k, where $ k $ is the number of prime factors of $ n $ , and $ p_i $ are the distinct prime factors. You may print the factors in any order. Hacks input For hacks, use the following format:. The first should contain $ k $ ( $ 2 \leq k \leq 10 $ ) — the number of prime factors of $ n $ . The second should contain $ k $ space separated integers $ p_1, p_2, \dots, p_k $ ( $ 21 \leq n \leq 2^{1024} $ ) — the prime factors of $ n $ . All prime factors have to be of form $ 4x + 3 $ for some integer $ x $ . They all have to be distinct. Interaction After printing a query do not forget to output end of line and flush the output. Otherwise you will get Idleness limit exceeded. To do this, use: - fflush(stdout) or cout.flush() in C++; - System.out.flush() in Java; - flush(output) in Pascal; - stdout.flush() in Python; - see documentation for other languages. The number of queries is not limited. However, your program must (as always) fit in the time limit. The run time of the interactor is also counted towards the time limit. The maximum runtime of each query is given below. - `+ x y` — up to $ 1 $ ms. - `- x y` — up to $ 1 $ ms. - `* x y` — up to $ 1 $ ms. - `/ x y` — up to $ 350 $ ms. - `sqrt x` — up to $ 80 $ ms. - `^ x y` — up to $ 350 $ ms. Note that the sample input contains extra empty lines so that it easier to read. The real input will not contain any empty lines and you do not need to output extra empty lines.

输入输出样例

输入样例 #1

21

7

17

15

17

11

-1

15

输出样例 #1

+ 12 16

- 6 10

* 8 15

/ 5 4

sqrt 16

sqrt 5

^ 6 12

! 2 3 7

说明

We start by reading the first line containing the integer $ n = 21 $ . Then, we ask for: 1. $ (12 + 16) \bmod 21 = 28 \bmod 21 = 7 $ . 2. $ (6 - 10) \bmod 21 = -4 \bmod 21 = 17 $ . 3. $ (8 \cdot 15) \bmod 21 = 120 \bmod 21 = 15 $ . 4. $ (5 \cdot 4^{-1}) \bmod 21 = (5 \cdot 16) \bmod 21 = 80 \bmod 21 = 17 $ . 5. Square root of $ 16 $ . The answer is $ 11 $ , as $ (11 \cdot 11) \bmod 21 = 121 \bmod 21 = 16 $ . Note that the answer may as well be $ 10 $ . 6. Square root of $ 5 $ . There is no $ x $ such that $ x^2 \bmod 21 = 5 $ , so the output is $ -1 $ . 7. $ (6^{12}) \bmod 21 = 2176782336 \bmod 21 = 15 $ . We conclude that our calculator is working, stop fooling around and realise that $ 21 = 3 \cdot 7 $ .