Odd Country

最新推荐文章于 2023-12-28 19:53:46 发布

转载最新推荐文章于 2023-12-28 19:53:46 发布 · 111 阅读

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/lightning34/p/4591761.html

Description

在一片美丽的大陆上有 $100000$ 个国家，记为 $1$ 到 $100000$。这里经济发达，有数不尽的账房，并且每个国家有一个银行。某大公司的领袖在这 $100000$ 个银行开户时都存了 $3$ 大洋，他惜财如命，因此会不时地派小弟 GFS 清点一些银行的存款或者让 GFS 改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算，也就是说领袖开户时的存款总和为 $3^{100000}$。这里发行的软妹面额是最小的 $60$ 个素数 $（p_1 = 2, p_2 = 3, \dots, p_{60} = 281）$，任何人的财产都只能由这 $60$ 个基本面额表示，即设某个人的财产为 $\text{fortune}$（正整数），则 $\text{fortune} = p_1^{k_1} \cdot p_2^{k_2} \cdot \dots \cdot p_{60}^{k_{60}}$。

领袖习惯将一段编号连续的银行里的存款拿到一个账房去清点，为了避免 GFS 串通账房叛变，所以他不会每次都选择同一个账房。GFS 跟随领袖多年已经摸清了门路,知道领袖选择账房的方式。如果领袖选择清点编号在 $[a, b]$ 内的银行财产，他会先对 $[a, b]$ 的财产求和（计为 $\text{product}$），然后在编号属于 $[1, \text{product}]$ 的账房中选择一个去清点存款，检验自己计算是否正确同时也检验账房与 GFS 是否有勾结。GFS 发现如果某个账房的编号 $\text{number}$ 与 $\text{product}$ 相冲，领袖绝对不会选择这个账房。怎样才算与 $\text{product}$ 不相冲呢？若存在整数 $x,y$ 使得 $\text{number} \cdot x + \text{product} \cdot y = 1$，那么我们称 $\text{number}$ 与 $\text{product}$ 不相冲，即该账房有可能被领袖相中。当领袖又赚大钱了的时候，他会在某个银行改变存款，这样一来相同区间的银行在不同的时候算出来的 $\text{product}$ 可能是不一样的，而且领袖不会在某个银行的存款总数超过 $1000000$。

现在 GFS 预先知道了领袖的清点存款与变动存款的计划，想请你告诉他，每次清点存款时领袖有多少个账房可以供他选择，当然这个值可能非常大，GFS 只想知道对 $19961993$ 取模后的答案。

Input Format

第一行一个整数 $x$ 表示领袖清点和变动存款的总次数。

接下来 $x$ 行，每行 $3$ 个整数 $a_i, b_i, c_i$。$a_i$ 为 $0$ 时表示该条记录是清点计划，领袖会清点 $b_i$ 到 $c_i$ 的银行存款，你需要对该条记录计算出 GFS 想要的答案。$a_i$ 为 $1$ 时表示该条记录是存款变动，你要把银行 $b_i$ 的存款改为 $c_i$，不需要对该记录进行计算。

Output Format

输出若干行，每行一个数，表示那些年的答案。

Sample Input

Sample Output

Hint

$20\%$ 的数据满足: $x \leq 10000$，当 $a_i = 0$ 时 $0 \leq c_i − b_i \leq 100$; 每个 $\text{product} \leq 10^{18}$;

$30\%$ 的数据满足：$x \leq 50000$，当 $a_i = 0$ 时 $0 \leq c_i − b_i \leq 10000$;

$50\%$ 的数据满足：$x \leq 100000$，当 $a_i = 0$ 时 $0 \leq c_i − b_i \leq 100000$;

以上数据不重合，$20\% +30\% + 50\% = 100\%$，且保证满足题干要求。

  1 #include <cstdio>
  2 #include <cstring>
  3 
  4 #define MOD 19961993ll
  5 #define MID (left + right >> 1)
  6 
  7 int f[300], p[100], pn;
  8 int n, q1, q2, q3;
  9 
 10 long long prod[400000], has[400000], rev[100], erci[100];
 11 
 12 long long tmp, tmphas;
 13 
 14 long long pow(long long num, long long mi)
 15 {
 16     if (mi == 0) return 1;
 17     if (mi == 1) return num;
 18     if (mi == 2) return num * num % MOD;
 19     return pow(pow(num, mi / 2), 2) * pow(num, mi % 2) % MOD;
 20 }
 21 
 22 void build(int i, int left, int right)
 23 {
 24     has[i] = 2;
 25     if (left != right)
 26     {
 27         build(i << 1, left, MID);
 28         build(i << 1 | 1, MID + 1, right);
 29         prod[i] = prod[i << 1] * prod[i << 1 | 1] % MOD;
 30         return;
 31     }
 32     prod[i] = 3;
 33 }
 34 
 35 void modify(int i, int left, int right, int pos, int data)
 36 {
 37     if (left == right)
 38     {
 39         prod[i] = data % MOD;
 40         has[i] = 0;
 41         for (int j = 0; j < 60; j++)
 42         {
 43             if (data % p[j] == 0)
 44             {
 45                 has[i] |= erci[j];
 46             }
 47         }
 48         return;
 49     }
 50     if (pos <= MID)
 51     {
 52         modify(i << 1, left, MID, pos, data);
 53     }
 54     else
 55     {
 56         modify(i << 1 | 1, MID + 1, right, pos, data);
 57     }
 58     prod[i] = prod[i << 1] * prod[i << 1 | 1] % MOD;
 59     has[i] = has[i << 1] | has[i << 1 | 1];
 60 }
 61 
 62 void query(int i, int left, int right, int ql, int qr)
 63 {
 64     if (ql <= left && qr >= right)
 65     {
 66         tmp = tmp * prod[i] % MOD;
 67         tmphas |= has[i];
 68         return;
 69     }
 70     if (qr < left || ql > right)
 71     {
 72         return;
 73     }
 74     if (ql <= MID)
 75     {
 76         query(i << 1, left, MID, ql, qr);
 77     }
 78     if (qr > MID)
 79     {
 80         query(i << 1 | 1, MID + 1, right, ql, qr);
 81     }
 82 }
 83 
 84 int main()
 85 {
 86     pn = 0;
 87     for (int i = 2; i < 300; i++)
 88     {
 89         if (!f[i])
 90         {
 91             rev[pn] = (i - 1) * pow(i, MOD - 2) % MOD;
 92             erci[pn] = (pn == 0 ? 1 : erci[pn - 1] << 1);
 93             p[pn] = i;
 94             pn++;
 95             for (int j = i + i; j < 300; j += i)
 96             {
 97                 f[j] = i;
 98             }
 99         }
100     }
101     build(1, 1, 100000);
102     scanf("%d", &n);
103     for (int i = 0; i < n; i++)
104     {
105         scanf("%d%d%d", &q1, &q2, &q3);
106         if (q1 == 0)
107         {
108             tmp = 1;
109             tmphas = 0;
110             query(1, 1, 100000, q2, q3);
111             for (int j = 0; j < 60; j++)
112             {
113                 if (tmphas & erci[j])
114                 {
115                     tmp = tmp * rev[j] % MOD;
116                 }
117             }
118             printf("%lld\n", tmp);
119         }
120         else
121         {
122             modify(1, 1, 100000, q2, q3);
123         }
124     }
125 }

转载于:https://www.cnblogs.com/lightning34/p/4591761.html