13af6ab5fSopenharmony_ci/*
23af6ab5fSopenharmony_ci * Copyright (c) 2024 Huawei Device Co., Ltd.
33af6ab5fSopenharmony_ci * Licensed under the Apache License, Version 2.0 (the "License");
43af6ab5fSopenharmony_ci * you may not use this file except in compliance with the License.
53af6ab5fSopenharmony_ci * You may obtain a copy of the License at
63af6ab5fSopenharmony_ci *
73af6ab5fSopenharmony_ci * http://www.apache.org/licenses/LICENSE-2.0
83af6ab5fSopenharmony_ci *
93af6ab5fSopenharmony_ci * Unless required by applicable law or agreed to in writing, software
103af6ab5fSopenharmony_ci * distributed under the License is distributed on an "AS IS" BASIS,
113af6ab5fSopenharmony_ci * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
123af6ab5fSopenharmony_ci * See the License for the specific language governing permissions and
133af6ab5fSopenharmony_ci * limitations under the License.
143af6ab5fSopenharmony_ci */
153af6ab5fSopenharmony_ci
163af6ab5fSopenharmony_ciexport class BitopsNSieveBits {
173af6ab5fSopenharmony_ci  static primes(isPrime: int[], n1: int, n2: int): void {
183af6ab5fSopenharmony_ci    let i: int;
193af6ab5fSopenharmony_ci    let m: int = n2 << n1;
203af6ab5fSopenharmony_ci    let size: int = m + 31 >> 5;
213af6ab5fSopenharmony_ci
223af6ab5fSopenharmony_ci    for (i = 0; i < size; i++) {
233af6ab5fSopenharmony_ci      isPrime[i] = -1;
243af6ab5fSopenharmony_ci    }
253af6ab5fSopenharmony_ci
263af6ab5fSopenharmony_ci    for (i = 2; i < m; i++) {
273af6ab5fSopenharmony_ci      if ((isPrime[i >> 5] & 1 << (i & 31)) != 0) {
283af6ab5fSopenharmony_ci        for (let j: int = i + i; j < m; j += i) {
293af6ab5fSopenharmony_ci          isPrime[j >> 5] &= ~(1 << (j & 31));
303af6ab5fSopenharmony_ci        }
313af6ab5fSopenharmony_ci      }
323af6ab5fSopenharmony_ci    }
333af6ab5fSopenharmony_ci  }
343af6ab5fSopenharmony_ci
353af6ab5fSopenharmony_ci  private static sieve(n1: int, n2: int): int[] {
363af6ab5fSopenharmony_ci    // Not parsed new int[...]
373af6ab5fSopenharmony_ci    let isPrime: int[] = new int[(n2 << n1) + 31 >> 5];
383af6ab5fSopenharmony_ci    BitopsNSieveBits.primes(isPrime, n1, n2);
393af6ab5fSopenharmony_ci    return isPrime;
403af6ab5fSopenharmony_ci  }
413af6ab5fSopenharmony_ci
423af6ab5fSopenharmony_ci  n1: int = 4;
433af6ab5fSopenharmony_ci  n2: int = 10000;
443af6ab5fSopenharmony_ci  static readonly expected: long = -1286749544853;
453af6ab5fSopenharmony_ci
463af6ab5fSopenharmony_ci  public run(): void {
473af6ab5fSopenharmony_ci    let result: int[] = BitopsNSieveBits.sieve(this.n1, this.n2);
483af6ab5fSopenharmony_ci    let sum: long = 0;
493af6ab5fSopenharmony_ci    for (let i: int = 0; i < result.length; ++i) {
503af6ab5fSopenharmony_ci      sum += result[i];
513af6ab5fSopenharmony_ci    }
523af6ab5fSopenharmony_ci
533af6ab5fSopenharmony_ci    assert sum == BitopsNSieveBits.expected: "Incorrect result";
543af6ab5fSopenharmony_ci  }
553af6ab5fSopenharmony_ci}
563af6ab5fSopenharmony_ci
573af6ab5fSopenharmony_cifunction main(): void {
583af6ab5fSopenharmony_ci  let a = new BitopsNSieveBits;
593af6ab5fSopenharmony_ci  a.run();
603af6ab5fSopenharmony_ci}
61