『算法-ACM竞赛-数学-数论』POJ1365——Prime Lan

『算法-ACM 竞赛-数学-数论』POJ1365——Prime Lan

Description

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
Everybody in the Prime Land is using a prime base number system. In this system, each positive integer x is represented as follows: Let {pi}i=0,1,2,... denote the increasing sequence of all prime numbers. We know that x > 1 can be represented in only one way in the form of product of powers of prime factors. This implies that there is an integer kx and uniquely determined integers ekx, ekx-1, ..., e1, e0, (ekx > 0), that  The sequence

(ekx, ekx-1, ... ,e1, e0)



is considered to be the representation of x in prime base number system.

It is really true that all numerical calculations in prime base number system can seem to us a little bit unusual, or even hard. In fact, the children in Prime Land learn to add to subtract numbers several years. On the other hand, multiplication and division is very simple.

Recently, somebody has returned from a holiday in the Computer Land where small smart things called computers have been used. It has turned out that they could be used to make addition and subtraction in prime base number system much easier. It has been decided to make an experiment and let a computer to do the operation ``minus one''.

Help people in the Prime Land and write a corresponding program.

For practical reasons we will write here the prime base representation as a sequence of such pi and ei from the prime base representation above for which ei > 0. We will keep decreasing order with regard to pi.

Input

1
The input consists of lines (at least one) each of which except the last contains prime base representation of just one positive integer greater than 2 and less or equal 32767. All numbers in the line are separated by one space. The last line contains number 0.

Output

1
The output contains one line for each but the last line of the input. If x is a positive integer contained in a line of the input, the line in the output will contain x - 1 in prime base representation. All numbers in the line are separated by one space. There is no line in the output corresponding to the last ``null'' line of the input.

Sample Input

1
2
3
4
17 1
5 1 2 1
509 1 59 1
0

Sample Output

1
2
3
2 4
3 2
13 1 11 1 7 1 5 1 3 1 2 1

给你 n 的因数分解，乘起来然后，分解 n-1

直接上重型武器

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
#include <stack>
#include <map>
#include <vector>
#include <string>
using namespace std;
typedef long long ll;
ll pr;
ll pmod(ll a, ll b, ll p) { return (a * b - (ll)((long double)a / p * b) * p + p) % p; } //普通的快速乘会T
ll gmod(ll a, ll b, ll p)
{
	ll res = 1;
	while (b)
	{
		if (b & 1)
			res = pmod(res, a, p);
		a = pmod(a, a, p);
		b >>= 1;
	}
	return res;
}
inline ll gcd(ll a, ll b)
{ //听说二进制算法特快
	if (!a)
		return b;
	if (!b)
		return a;
	int t = __builtin_ctzll(a | b);
	a >>= __builtin_ctzll(a);
	do
	{
		b >>= __builtin_ctzll(b);
		if (a > b)
		{
			ll t = b;
			b = a, a = t;
		}
		b -= a;
	} while (b);
	return a << t;
}
bool Miller_Rabin(ll n)
{
	if (n == 46856248255981ll || n < 2)
		return false; //强伪素数
	if (n == 2 || n == 3 || n == 7 || n == 61 || n == 24251)
		return true;
	if (!(n & 1) || !(n % 3) || !(n % 61) || !(n % 24251))
		return false;
	ll m = n - 1, k = 0;
	while (!(m & 1))
		k++, m >>= 1;
	for (int i = 1; i <= 20; ++i) // 20为Miller-Rabin测试的迭代次数
	{
		ll a = rand() % (n - 1) + 1, x = gmod(a, m, n), y;
		for (int j = 1; j <= k; ++j)
		{
			y = pmod(x, x, n);
			if (y == 1 && x != 1 && x != n - 1)
				return 0;
			x = y;
		}
		if (y != 1)
			return 0;
	}
	return 1;
}
ll _abs(ll a){
	if(a>=0) return a;
	else return -a;
}
ll Pollard_Rho(ll x)
{
	ll n = 0, m = 0, t = 1, q = 1, c = rand() % (x - 1) + 1;
	for (ll k = 2;; k <<= 1, m = n, q = 1)
	{
		for (ll i = 1; i <= k; ++i)
		{
			n = (pmod(n, n, x) + c) % x;
			q = pmod(q, _abs(m - n), x);
		}
		t = gcd(x, q);
		if (t > 1)
			return t;
	}
}
map<long long, int> m;
void fid(ll n)
{
	if (n == 1)
		return;
	if (Miller_Rabin(n))
	{
		pr = max(pr, n);
		m[n]++;
		return;
	}
	ll p = n;
	while (p >= n)
		p = Pollard_Rho(p);
	fid(p);
	fid(n / p);
}
int main()
{
	ll n, a;
    int b;
    while (1)
    {
        n = 1;
        while (true)
        {
            scanf("%lld", &a);
            if (a == 0)
                return 0;
            char c = getchar();
            if (c == '\n')
                break;
            scanf("%d", &b);
            c = getchar();
            for (int i = 1; i <= b; i++)
                n *= a;
            if (c == '\n')
                break;
        }
        n--;
		m.clear();
		pr = 0;
		fid(n);
		for (map<long long, int>::iterator c = m.end(); c != m.begin();)
		{
			--c;
			printf("%lld %d ", c->first, c->second);
		}
		puts(" ");
	}
	return 0;
}
算法 > ACM竞赛 > 数学 > 数论
『算法-ACM竞赛-数学-数论』POJ1365——Prime Lan
https://chiamzhang.github.io/2024/06/29/『算法-ACM竞赛-数学-数论』POJ1365——Prime Land/
Author
Chiam
Posted on
June 29, 2024
Licensed under