hdu 6825 Set1

题目

在这里插入图片描述
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6825

思路

这个大佬说的很清楚->大佬博客

再添加以下基础在这里插入图片描述

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int mod =  998244353;
const int N=5e6+5;
int a[N],ans[N],inv_a[N],inv[N],C_tmp[N],A_tmp[N];
int quick_pow(int a,int b){
	int res=1;
	while(b){
		if(b&1) res=res*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return res;
}
int C(int m,int n){
	return a[m]*inv_a[m-n]%mod*inv_a[n]%mod;
}
int A(int m,int n){
	return a[m]*inv_a[m-n]%mod;
}
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	int t;
	cin>>t;
	a[0]=1;
	for(int i=1;i<=N;i++){
		a[i]=(a[i-1]*i)%mod;//n! 
	}
//	cout<<1<<endl;
	inv_a[N]=quick_pow(a[N],mod-2);
	for(int i=N-1;i>=0;i--){
		inv_a[i]=inv_a[i+1]*(i+1)%mod;//n!的逆元 
	}
//	cout<<2<<endl;
	for(int i=1;i<=N;i++){
		inv[i]=inv_a[i]*a[i-1]%mod;//单个n的逆元 
	}
//	cout<<3<<endl;
	C_tmp[0]=1;
	for(int i=1;i<=2500000;i++){
		C_tmp[i]=C_tmp[i-1]*C(2*i,2)%mod;
	}
//	cout<<4<<endl;
	A_tmp[1]=1;
	for(int i=3;i<=5000000;i+=2){
		A_tmp[i]=A_tmp[i-2]*(i-1)%mod;
	}
//	cout<<5<<endl;
	A_tmp[5000000-1]=quick_pow(A_tmp[5000000-1],mod-2);
	for(int i=5000000-3;i>=0;i-=2){
		A_tmp[i]=A_tmp[i+2]*(i+1)%mod;
	}
//	cout<<6<<endl;
	while(t--){
//		cout<<7<<endl;
		int n;
		cin>>n;
		for(int i=1;i<=n/2;i++) cout<<"0 ";
		for(int i=n/2+1;i<=n;i++){
			int s2=n-i,s1=i-1;
			cout<<C(s1,s2)*A(s2,s2)%mod*C_tmp[(s1-s2)/2]%mod*inv_a[(s1-s2)/2]%mod*A_tmp[n]%mod<<" \n"[i==n];
		}
	}
    return 0;
}

©️2020 CSDN 皮肤主题: 深蓝海洋 设计师:CSDN官方博客 返回首页