結果

問題 No.206 数の積集合を求めるクエリ
ユーザー fgwiebfaoishfgwiebfaoish
提出日時 2020-11-05 22:11:50
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 611 ms / 7,000 ms
コード長 4,431 bytes
コンパイル時間 2,549 ms
コンパイル使用メモリ 108,612 KB
実行使用メモリ 45,244 KB
最終ジャッジ日時 2023-09-29 17:34:32
合計ジャッジ時間 14,709 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 65 ms
23,856 KB
testcase_01 AC 63 ms
19,736 KB
testcase_02 AC 64 ms
21,816 KB
testcase_03 AC 64 ms
21,796 KB
testcase_04 AC 63 ms
21,788 KB
testcase_05 AC 64 ms
23,836 KB
testcase_06 AC 81 ms
24,420 KB
testcase_07 AC 81 ms
20,364 KB
testcase_08 AC 81 ms
22,388 KB
testcase_09 AC 81 ms
22,528 KB
testcase_10 AC 65 ms
21,908 KB
testcase_11 AC 65 ms
21,788 KB
testcase_12 AC 83 ms
24,572 KB
testcase_13 AC 81 ms
20,520 KB
testcase_14 AC 82 ms
22,404 KB
testcase_15 AC 82 ms
20,380 KB
testcase_16 AC 81 ms
22,516 KB
testcase_17 AC 593 ms
43,668 KB
testcase_18 AC 545 ms
35,288 KB
testcase_19 AC 586 ms
42,520 KB
testcase_20 AC 549 ms
35,792 KB
testcase_21 AC 547 ms
36,680 KB
testcase_22 AC 547 ms
36,636 KB
testcase_23 AC 587 ms
42,652 KB
testcase_24 AC 611 ms
43,524 KB
testcase_25 AC 598 ms
42,240 KB
testcase_26 AC 574 ms
43,344 KB
testcase_27 AC 555 ms
44,720 KB
testcase_28 AC 575 ms
44,840 KB
testcase_29 AC 582 ms
42,828 KB
testcase_30 AC 566 ms
45,244 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
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