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