結果

問題 No.206 数の積集合を求めるクエリ
コンテスト
ユーザー fgwiebfaoish
提出日時 2020-11-05 22:11:50
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 574 ms / 7,000 ms
コード長 4,431 bytes
コンパイル時間 4,569 ms
コンパイル使用メモリ 116,416 KB
実行使用メモリ 40,704 KB
最終ジャッジ日時 2024-07-22 11:35:31
合計ジャッジ時間 11,751 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.Collections.Generic;
using System.Collections;
using System.Collections.Specialized;
using System.Linq;
using System.Text;
using System.IO;
using System.Reflection;
using static System.Math;
using System.Numerics;
using System.Threading;
using System.Runtime.CompilerServices;
using System.Diagnostics;
//using nint=System.Int32;
static class Program{
	const int mod=(int)1e9+7;
	static void Main(){
		Sc sc=new Sc();
		var s=sc.Ia;
		var a=sc.Ia;
		var b=sc.Ia;
		var c=new int[s[2]+1];
		var d=new int[s[2]+1];
		for(int i = 0;i<a.Length;i++) {c[a[i]]++;}
		for(int i = 0;i<b.Length;i++) {d[s[2]-b[i]+1]++;}
		Console.WriteLine("{0}",string.Join("\n",Ex.Fft(c,d).Skip(s[2]+1).Take(sc.I)));
	}
}
static class Ex{
	private static long[] ga1=new long[24],ga2=new long[24];
	private const int mod=998244353;
	public static int[] Fft(int[] a,int[] b){
		int m=1,p=0;
		while(m<a.Length+b.Length){m<<=1;p++;}
		long[] aa=new long[m],ba=new long[m],c=new long[m];
		var d=new int[a.Length+b.Length-1];
		for(int i = 0;i<a.Length;i++) {aa[i]=a[i];}
		for(int i = 0;i<b.Length;i++) {ba[i]=b[i];}
		ga1[23]=15311432;
		ga2[23]=469870224;
		for(int i = 22;i>0;i--) {
			ga1[i]=ga1[i+1]*ga1[i+1]%mod;
			ga2[i]=ga2[i+1]*ga2[i+1]%mod;
		}
		Dft(aa,p);
		Dft(ba,p);
		for(int i = 0;i<m;i++) {c[i]=aa[i]*ba[i]%mod;}
		Idft(c,p);
		var y=Mi(c.Length);
		for(int i = 0;i<d.Length;i++) {d[i]=(int)(c[i]*y%mod);}
		return d;
	}
	public static void Dft(long[] a,int h){
		for(int i=0;i<a.Length;i++) {
			int j=0;
			for (int k=0;k<h;k++){j|=(i>>k&1)<<(h-1-k);}
			if(i<j){(a[i],a[j])=(a[j],a[i]);}
		}
		for(int b=1,p=1;b<a.Length;b<<=1,p++) {
			long w=1,f=ga1[p];
			for (int j = 0;j<b;j++,w=w*f%mod) {
				for (int k = 0;k<a.Length;k+=b<<1) {
					long s=a[j+k],t=a[j+k+b]*w%mod;
					a[j+k]=(s+t)%mod; 
					a[j+k+b]=(mod+s-t)%mod;
				}
			}
		}
	}
	public static void Idft(long[] a,int h){
		for(int i=0;i<a.Length;i++) {
			int j=0;
			for (int k=0;k<h;k++){j|=(i>>k&1)<<(h-1-k);}
			if(i<j){(a[i],a[j])=(a[j],a[i]);}
		}
		for(int b=1,p=1;b<a.Length;b<<=1,p++) {
			long w=1,f=ga2[p];
			for (int j = 0;j<b;j++,w=w*f%mod) {
				for (int k = 0;k<a.Length;k+=b<<1) {
					long s=a[j+k],t=a[j+k+b]*w%mod;
					a[j+k]=(s+t)%mod; 
					a[j+k+b]=(mod+s-t)%mod;
				}
			}
		}
	}
	private static int Mi(int a){
		a=(a+mod)%mod;
		int b=mod,u=1,v=0;
		while(b>0){
			int t=a/b;
			a-=t*b;(a,b)=(b,a);
			u-=t*v;(u,v)=(v,u);
		}
		u%=mod;
		if(u<0){u+=mod;}
		return u%mod;
	}
}
public class Sc{
	public int I{get{return int.Parse(Console.ReadLine());}}
	public long L{get{return long.Parse(Console.ReadLine());}}
	public double D{get{return double.Parse(Console.ReadLine());}}
	public string S{get{return Console.ReadLine();}}
	public int[] Ia{get{return Array.ConvertAll(Console.ReadLine().Split(),int.Parse);}}
	public long[] La{get{return Array.ConvertAll(Console.ReadLine().Split(),long.Parse);}}
	public double[] Da{get{return Array.ConvertAll(Console.ReadLine().Split(),double.Parse);}}
	public string[] Sa{get{return Console.ReadLine().Split();}}
	public object[] Oa{get{return Console.ReadLine().Split();}}
	public int[] Ia2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),int.Parse);}}
	public int[] Ia3(string a,string b){return Array.ConvertAll((a+Console.ReadLine()+b).Split(),int.Parse);}
	public int[] Ia3(int a){return Array.ConvertAll((Console.ReadLine()+" "+a.ToString()).Split(),int.Parse);}
	public long[] La2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),long.Parse);}}
	public long[] La3(string a,string b){return Array.ConvertAll((a+Console.ReadLine()+b).Split(),long.Parse);}
	public long[] La3(int a){return Array.ConvertAll((Console.ReadLine()+" "+a.ToString()).Split(),long.Parse);}
	public double[] Da2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),double.Parse);}}
	public double[] Da3(string a,string b){return Array.ConvertAll((a+Console.ReadLine()+b).Split(),double.Parse);}
	public T[] Arr<T>(int n,Func<T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f();}return a;}
	public T[] Arr<T>(int n,Func<int,T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(i);}return a;}
	public T[] Arr<T>(int n,Func<string[],T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(Console.ReadLine().Split());}return a;}
	public T[] Arr<T>(int n,Func<int,string[],T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(i,Console.ReadLine().Split());}return a;}
}
0