ural 2062 Ambitious Experiment

2062. Ambitious Experiment

Time limit: 3.0 second
Memory limit: 128 MB

During several decades, scientists from planet Nibiru are working to create an engine that would allow spacecrafts to fall into hyperspace and move there with superluminal velocity. To check whether their understanding of properties of hyperspace is right, scientists have developed the following experiment.

A chain of n particles is placed in hyperspace. Positions of particles in the chain are numbered from 1 to n. Initially, ith particle has charge ai.

According to the current theory, if particle number i got special radiation with power d, oscillations would spread by hyperspace and increase by d charge of particles with numbers i, 2i, 3iand so on (i.e. with numbers divisible by i).

Using a special device, scientists can direct the radiation of the same power at a segment of adjacent particles. For example, suppose that initially there were 6 particles with zero charges, and scientists have sent radiation with power five to particles with numbers 2 and 3. Then charge of 2nd, 3rd, and 4th particles will increase to five, and charge of 6th particle will increase to ten (the oscillations will reach it twice). Charge of other particles won’t change.

Charge of particles can’t change without impact of the device.

During the experiment, the scientists plan to perform actions of the following types:

1. Measure current charge of the particle number i.
2. Direct radiation with power d at particles with numbers from l to r inclusive.

Your program will be given a list of performed actions. For every action of the first type the program should output value of expected charge of the particle calculated in accordance with the current theory described above.

If the expected charges of the particles coincide with charges measured during the experiment, it will turn out that scientists’ understanding of hyperspace is right, and they will be able to start building of the hyperdrives. Then inhabitants of Nibiru will finally meet their brothers from Earth in just a few years!

Input

The first line contains a single integer n — number of particles (1 ≤ n ≤ 3 · 105).

The second line contains n integers ai separated by spaces — initial charges of the particles (0 ≤ ai ≤ 106).

The third line contains a single integer q — number of actions in the experiment (1 ≤ q ≤ 3 · 105).

Each of the following q lines contain two or four integers — a description of the next action in one of the following formats:

• 1 i — measure current charge of the particle number i (1 ≤ i ≤ n).
• 2 l r d — direct radiation with power d at particles with numbers from l to r inclusive (1 ≤l ≤ r ≤ n, 0 ≤ d ≤ 106).

Output

For each query output the expected charge of the ith particle.

Samples

input	output
3 1 2 3 2 2 1 3 5 1 2	12
6 1 2 1 4 5 6 5 2 2 4 2 1 3 1 4 2 3 5 1 1 5	3 8 6

Problem Author: Alexey Danilyuk (prepared by Nikita Sivukhin)
Problem Source: Ural Regional School Programming Contest 2015

Tags: data structures  (hide tags for unsolved problems)

Difficulty: 478

 

 

题意：给n个数，两种操作：1、询问x的值，2、修改l~r的值，对于每个l<=i<=r的i，都对 i,2i,3i....这些点加d

分析：首先很容易想到分块做法

对于每个点x，它对于其他点（2x,3x,....）的贡献都是一样的

那么对于每个查询的点x，对它有贡献的点就是是他的因数的那些点

那么修改就分段sqrt(n)暴力修改

查询就sqrt(x)枚举约数，查询到每一个贡献点的值加起来

 

 /**
 Create By yzx - stupidboy
 */
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <cmath>
 #include <deque>
 #include <vector>
 #include <queue>
 #include <iostream>
 #include <algorithm>
 #include <map>
 #include <set>
 #include <ctime>
 #include <iomanip>
 using namespace std;
 typedef long long LL;
 typedef double DB;
 #define MIT (2147483647)
 #define INF (1000000001)
 #define MLL (1000000000000000001LL)
 #define sz(x) ((int) (x).size())
 #define clr(x, y) memset(x, y, sizeof(x))
 #define puf push_front
 #define pub push_back
 #define pof pop_front
 #define pob pop_back
 #define mk make_pair

 inline int getInt()
 {
     int ret = ;
     char ch = ' ';
     bool flag = ;
     while(!(ch >= '' && ch <= ''))
     {
         if(ch == '-') flag ^= ;
         ch = getchar();
     }
     while(ch >= '' && ch <= '')
     {
         ret = ret *  + ch - '';
         ch = getchar();
     }
     return flag ? -ret : ret;
 }

 const int N = , M = ;
 int n, m;
 LL block[N / M + ][M], tag[N / M + ], arr[N];

 inline int getBlockIndex(int x)
 {
     return x / M;
 }

 inline int getIndex(int x)
 {
     return x % M;
 }

 inline void input()
 {
     cin >> n;
     for(int i = ; i < n; i++) cin >> arr[i];
 }

 inline LL work(int x)
 {
     int b = getBlockIndex(x), idx = getIndex(x);
     return block[b][idx] + tag[b];
 }

 inline LL query(int x)
 {
     x++;
     LL ret = ;
     for(int i = ; i * i <= x; i++)
         if(x % i == )
         {
             ret += work(i - );
             if(x / i != i) ret += work(x / i - );
         }
     return ret;
 }

 inline void change(int l, int r, int d)
 {
     int left = getBlockIndex(l), right = getBlockIndex(r);
     for(int i = left + ; i <= right - ; i++) tag[i] += d;
     if(left < right)
     {
         for(int i = getIndex(l); i < M; i++) block[left][i] += d;
         for(int i = ; i <= getIndex(r); i++) block[right][i] += d;
     }
     else
     {
         for(int i = getIndex(l); i <= getIndex(r); i++)
             block[left][i] += d;
     }
 }

 inline void solve()
 {
     for(cin >> m; m--; )
     {
         int opt, l, r, x;
         cin >> opt;
         if(opt == )
         {
             cin >> x;
             x--;
             cout << query(x) + arr[x] << "\n";
         }
         else
         {
             cin >> l >> r >> x;
             l--, r--;
             change(l, r, x);
         }
     }
 }

 int main()
 {
     ios::sync_with_stdio();
     input();
     solve();
     return ;
 }