2020 CCPC Wannafly Winter Camp Day1 H

題目好像不能公開,所以我只寫個題解…
開場簽到題丟給隊友,轉而開了H題。
H的思路是選定一個數y通過gcd(y,k)來 在1到n 內唯一地確定k的值。
不難發現,通過gcd來唯一確定k
所以:
1.y必須是k的倍數,否則無法區分gcd(y,k)與k。
2.對於1~(n/kk)內的每一個素數p,y/k必須是p的倍數,否則無法區分k與pk。
最終答案等於k
(1~n/k*k內所有的質數)
質數累乘會超出longlong需要使用大數(double浮點數的精度不足)

#include <iostream>
#include <algorithm>
#include <string.h>
using namespace std;
const int maxn = 505;
int su[maxn] = {0};
struct BigInt//大數板子
{
    const static int mod = 10000;
    const static int DLEN = 4;
    int a[600], len;
    BigInt()
    {
        memset(a, 0, sizeof(a));
        len = 1;
    }
    BigInt(int v)
    {
        memset(a, 0, sizeof(a));
        len = 0;
        do
        {
            a[len++] = v % mod;
            v/=mod;
        }while(v);
    }
    BigInt operator *(const BigInt &b)const
    {
        BigInt res;
        for(int i = 0; i < len; i++)
        {
            int up = 0;
            for(int j = 0; j < b.len; j++)
            {
                int temp = a[i]*b.a[j] + res.a[i+j] + up;
                res.a[i+j] = temp % mod;
                up = temp / mod;
            }
            if(up != 0)
            {
                res.a[i+b.len] = up;
            }
        }
        res.len = len + b.len;
        while(res.a[res.len-1] == 0 && res.len > 1) res.len--;
        return res;
    }

    void output()
    {
        printf("%d", a[len-1]);
        for(int i = len-2; i >= 0; i--)
        {
            printf("%04d", a[i]);
            
        }
        printf("\n");
    }
};
void aishai()//埃篩用來篩素數
{
	su[1] = 1;
	for(int i = 2; i <= 250; i++)
	{
		if(su[i] == 0)
		{
			for(int j = 2; i * j <= 500; j++)
			{
				su[i*j] = 1;
			}
		}
	}
	return ;
}
int main()
{
	aishai();
	int t;
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);//關掉同步流,加快輸入
	cin >> t;
	while(t--)
	{
		BigInt ans(1);
		int n, k;
		cin >> n >> k;
		BigInt kda(k);
		ans = ans * kda;
		for(int i = 2; i <= n / k; i++)
		{
			if(su[i] == 0)
			{
				BigInt t(i);
				ans = ans *  t;
			}
		}
		ans.output();
	}
	return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章