HDU - 6287 口算训练 (质因子分解 + 二分)
生活随笔
收集整理的這篇文章主要介紹了
HDU - 6287 口算训练 (质因子分解 + 二分)
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題意:給一個序列,問 [l, r] 連續(xù)區(qū)間的乘積能否被 d 整除
思路:用 vector 記錄每個質(zhì)因子在序列中出現(xiàn)的位置,upper_bound(vec[p].begin(), vec[p].end, r) - lower_bound(vec[p].begin(), vec[p].end(), l) 即為質(zhì)數(shù) p 在區(qū)間 [l, r] 中出現(xiàn)的次數(shù)。(可以不篩素數(shù),直接試除
大數(shù)模板和Java都會超時
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef unsigned long long ull; typedef double db; const db eps = 1e-6; const int N = 1e5 + 7; const int inf = 0x3f3f3f3f; const int mod = 1e9 + 7;int t, n, m, l, r, d, a[N], pri[N], tot; vector<int>vec[N]; bool vis[N];void init() {tot = 0;vis[0] = vis[1] = 1;for(int i = 2; i < N; ++i) {if(!vis[i]) {pri[++tot] = i;for(int j = i + i; j < N; j += i) vis[j] = 1;}} }void solve(int x, int id) {int tmp = x;for(int i = 1; pri[i] * pri[i] <= tmp && i <= tot; ++i) {while(x % pri[i] == 0) {x /= pri[i];vec[pri[i]].emplace_back(id);}}if(x > 1) vec[x].emplace_back(id); }bool check(int l, int r, int p, int cnt) {int ans = upper_bound(vec[p].begin(), vec[p].end(), r) - lower_bound(vec[p].begin(), vec[p].end(), l);return ans >= cnt; }int main() {init();scanf("%d", &t);while(t--) {scanf("%d%d", &n, &m);for(int i = 1; i < N; ++i) vec[i].clear();for(int i = 1; i <= n; ++i) {scanf("%d", &a[i]);solve(a[i], i);}while(m--) {scanf("%d%d%d", &l, &r, &d);int dd = d;bool flag = 1;for(int i = 1; pri[i] * pri[i] <= dd && i <= tot; ++i) {if(d % pri[i] == 0) {int cnt = 0;while(d % pri[i] == 0) {d /= pri[i];++cnt;}if(!check(l, r, pri[i], cnt)) {flag = 0;break;}}}if(flag && d > 1 && !check(l, r, d, 1)) flag = 0;if(flag) printf("Yes\n");else printf("No\n");}}return 0; }總結(jié)
以上是生活随笔為你收集整理的HDU - 6287 口算训练 (质因子分解 + 二分)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: vba 定义类_类模块的定义及设计
- 下一篇: I2C总线中8bit和10bit有什么区