zoukankan      html  css  js  c++  java
  • Codeforces Round #694 (Div. 2)D. Strange Definition(hash+数学+思维)

    D. Strange Definition

    题目大意:

    给你(n)个数,(q)次询问。

    其中(x)(y)相邻需要满足(lcm(x,y) / gcd(x,y))是一个平方数。

    每经过1秒,相应的序列会变为和他相邻的集合的乘积。

    定义(d_i)是和(a_i)相邻的元素的个数。让你求最大的(d_i)

    思路:

    首先化简(lcm(x,y)/gcd(x,y)->x*y/{gcd(x,y)^2})我们发现要使其是一个平方数那么只有(x*y)是一个平方数时才满足。(x*y)是平方数只有其质因数分解相乘后所有的次幂都是偶数才满足条件,也就是说(x,y)质因数分解后的次幂是相同的,若为偶数次幂同为偶数次幂,若为奇数次幂同为奇数次幂。

    然后我们将集合的元素个数分为奇数个和偶数个,因为如果元素个数是偶数个那么他们相乘的结果仍然是偶次幂。

    如果是奇数个那么他们相乘的结果仍然有奇数次幂。

    处理时候特别需要注意(mt[1])也就是1的时候和偶次幂的相乘也仍然是平方数。

    辨别其时我们用了hash -> mt[res]的思想,因为质数相乘是不同的结果。
    情况仅分为0秒和1秒的情况,大于1秒时和1秒相同,因为集合内为偶数个和奇数个的不会再变化。

    #include<bits/stdc++.h>
    
    using namespace std;
    const int N = 3e5 + 100;
    
    typedef long long LL;
    typedef unsigned long long ull;
    int a[N];
    int prime[N];
    bool vis[N];
    int cnt = 0;
    void primes(int n) {
        vis[1] = 1;
        for(int i = 2; i <= n; i++) {
            if(!vis[i]) prime[++cnt] = i;
            for(int j = 1; j <= cnt; j++) {
                if(prime[j] > n / i)
                    break;
                vis[i * prime[j] * 1ll] = 1;
                if(i % prime[j] == 0)
                    break;
            }
        }
    }
    unordered_map<ull, int>mt;
    vector<int>v;
    void f(int x) {
    	int res = 1;
    	for (int i = 1; i <= cnt; ++i) {
    		int tot = 0;
    		while (x % prime[i] == 0) x /= prime[i], ++tot;
    		if (tot & 1) res *= prime[i];
    		if (x < 1LL * prime[i] * prime[i]) break;
    	}
    	if (x > 1) res *= x;
    	mt[res] += 1;
    	v.push_back(res);
    }
    void solve() {
    	mt.clear();v.clear();
    	int n; scanf("%d", &n);
    	for (int i = 1; i <= n; ++i) {
    		scanf("%d", &a[i]);
    		f(a[i]);
    	}
    	int t1 = 0, t2 = mt[1];
    	for (int i = 0; i < int(v.size()); ++i)t1 = max(t1, mt[v[i]]);
    	int ans = 0;
    	for (int i = 0; i < int(v.size()); ++i) {
    		if (v[i] == 1) continue;
    		if (mt[v[i]] % 2 == 0) {
    			t2 += mt[v[i]];
    			mt[v[i]] = 0;
    		}
    	}
    	t2 = max(t2, t1);
    	int q; scanf("%d", &q);
    	while (q--) {
    		int w; scanf("%d", &w);
    		printf("%d
    ", (w == 0 ? t1 : t2));
    	}
    }
    
    signed main() {
    	primes(N - 10);
        int T = 1;
    	scanf("%d",&T);
        while (T--) {
            solve();
        }
    
    }
    
  • 相关阅读:
    一个贼基础的 编码解码方式
    SQL 中循环、for循环、游标
    sql中判断是否存在 数据库、表、存储过程、函数
    sql 同步表或同步表的时候更改部分字段
    sql存储过程的建立
    POJ
    UCloud 的安全秘钥 (计蒜客初赛第五场)(待解决)
    UCloud 机房的网络搭建(计蒜客初赛第五场)
    2017 计蒜之道 初赛 第四场
    腾讯课堂的物理实验(2017计蒜客初赛第三场)
  • 原文地址:https://www.cnblogs.com/waryan/p/14288505.html
Copyright © 2011-2022 走看看