Submission #405790


Source Code Expand

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cctype>
#include<cstdlib>
#include<algorithm>
#include<bitset>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<map>
#include<set>
#include<stack>
#include<cmath>
#include<sstream>
#include<fstream>
#include<iomanip>
#include<ctime>
#include<complex>
#include<functional>
#include<climits>
#include<cassert>
#include<iterator>
using namespace std;
#define MOD 1000000007

namespace math{
	long long int __gcd(long long int a, long long int b){
		if (a > b){
			swap(a, b);
		}
		while (a){
			swap(a, b);
			a %= b;
		}
		return b;
	}
	long long int lcm(long long int a, long long int b){
		long long int g = __gcd(a, b);
		a /= g;
		return a*b;
	}
	long long int ppow(long long int i, long long int j){
		long long int res = 1LL;
		while (j){
			if (j & 1LL){
				res *= i;
				res %= MOD;
			}
			i *= i;
			i %= MOD;
			j >>= 1LL;
		}
		return res;
	}
	namespace factorial{
		vector<long long int> lo;
		vector<double> l2;
		void set_long(long long int b){
			if (lo.size()){
			}
			else{
				lo.push_back(1);
			}
			for (long long int i = lo.size(); i <= b; i++){
				lo.push_back(lo.back());
				lo.back() *= i;
				if (lo.back() >= MOD){
					lo.back() %= MOD;
				}
			}
		}
		void set_log(long long int b){
			if (l2.size()){
			}
			else{
				l2.push_back(log(0.0));
			}
			for (long long int i = l2.size(); i <= b; i++){
				l2.push_back(l2.back());
				l2.back() += log((double)(i));
			}
		}
		long long int get_long(int b){
			if (lo.size() <= b){
				set_long(b);
			}
			return lo[b];
		}
		double get_log(int b){
			if (l2.size() <= b){
				set_log(b);
			}
			return l2[b];
		}
	}
	namespace combination{
		long long int simpleC(long long int a, long long int b){
			if (a < b){
				return 0;
			}
			if (a - b < b){
				b = a - b;
			}
			long long int u = 1LL;
			for (long long int j = a; j >= a - b + 1LL; j--){
				u *= j;
				if (u >= MOD){
					u %= MOD;
				}
			}
			long long int s = 1LL;
			for (long long int i = 1LL; i <= b; i++){
				s *= i;
				if (s >= MOD){
					s %= MOD;
				}
			}
			return (u*ppow(s, MOD - 2)) % MOD;
		}
		long long int C(long long int a, long long int b){
			if (a < b){
				return 0;
			}
			long long int u = math::factorial::get_long(a);
			long long int s = math::factorial::get_long(b)*math::factorial::get_long(a - b);
			u %= MOD;
			s %= MOD;
			return (u*ppow(s, MOD - 2)) % MOD;
		}
		double logC(int a, int b){
			double u = math::factorial::get_log(a);
			double s = math::factorial::get_log(b) + math::factorial::get_log(a - b);
			return u - s;
		}
		long long int H(long long int a, long long int b){
			return math::combination::C(a + b - 1LL, b);
		}
		long long int simpleH(long long int a, long long int b){
			return math::combination::simpleC(a + b - 1LL, b);
		}
	}
	namespace prime{
		vector<long long int> prime;
		vector<long long int> use;  //smallest divisor
		void init(int b){
			use.assign(b + 1, 0);
			prime.clear();
			prime.push_back(2);
			use[2] = 2;
			for (int i = 3; i < use.size(); i += 2){
				if (use[i] == 0LL){
					prime.push_back(i);
					use[i] = i;
					for (int j = i * 2; j < use.size(); j += i){
						use[j] = i;
					}
				}
			}
		}
		vector<int> factorizetion(long long int num){
			vector<int> r;
			r.clear();
			for (int i = 0; i<prime.size() && prime[i] * prime[i] <= num; i++){
				while (num%prime[i] == 0LL){
					r.push_back(prime[i]);
					num /= prime[i];
				}
			}
			if (num > 1LL){
				r.push_back(num);
			}
			return r;
		}
		int size_of_factorization(long long int num){
			int cnt = 0;
			for (int i = 0; i<prime.size() && prime[i] * prime[i] <= num; i++){
				while (num%prime[i] == 0LL){
					cnt++;
					num /= prime[i];
				}
			}
			if (num > 1LL){
				cnt++;
			}
			return cnt;
		}
		long long int number_of_div(long long int num){
			long long int way = 1LL;
			long long int cnt = 0;
			for (int i = 0; i < prime.size() && prime[i] * prime[i] <= num; i++){
				cnt = 0;
				while (num%prime[i] == 0){
					cnt++;
					num /= prime[i];
				}
				way *= (cnt + 1LL);
			}
			if (num > 1LL){
				way *= 2LL;
			}
			return way;
		}
	}
}
using namespace math;
using namespace factorial;
int main(){
	int n;
	int k1;
	scanf("%d%d", &n, &k1);
	if (k1<n){
		puts("0");
		return 0;
	}
	int s = k1/n;
	int k = k1 - s*n;
	int kk = n - k;
	long long int ans = factorial::get_long(n);
	long long int anss = factorial::get_long(k);
	anss *= factorial::get_long(kk);
	anss %= MOD;
	ans *= ppow(anss, MOD - 2);
	ans %= MOD;
	printf("%lld\n", ans);
	return 0;
}

Submission Info

Submission Time
Task B - 高橋幼稚園
User Kmcode
Language C++ (GCC 4.9.2)
Score 80
Code Size 4819 Byte
Status WA
Exec Time 26 ms
Memory 932 KB

Compile Error

./Main.cpp: In function ‘int main()’:
./Main.cpp:213:24: warning: ignoring return value of ‘int scanf(const char*, ...)’, declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d%d", &n, &k1);
                        ^

Judge Result

Set Name Sample Subtask1 Subtask2
Score / Max Score 0 / 0 80 / 80 0 / 20
Status
AC × 2
WA × 1
AC × 21
AC × 21
WA × 12
Set Name Test Cases
Sample subtask0_sample_01.txt, subtask0_sample_02.txt, subtask0_sample_03.txt
Subtask1 subtask0_sample_01.txt, subtask0_sample_02.txt, subtask1_01.txt, subtask1_02.txt, subtask1_03.txt, subtask1_04.txt, subtask1_05.txt, subtask1_06.txt, subtask1_07.txt, subtask1_08.txt, subtask1_09.txt, subtask1_10.txt, subtask1_11.txt, subtask1_12.txt, subtask1_13.txt, subtask1_14.txt, subtask1_15.txt, subtask1_16.txt, subtask1_17.txt, subtask1_18.txt, subtask1_19.txt
Subtask2 subtask0_sample_01.txt, subtask0_sample_02.txt, subtask0_sample_03.txt, subtask1_01.txt, subtask1_02.txt, subtask1_03.txt, subtask1_04.txt, subtask1_05.txt, subtask1_06.txt, subtask1_07.txt, subtask1_08.txt, subtask1_09.txt, subtask1_10.txt, subtask1_11.txt, subtask1_12.txt, subtask1_13.txt, subtask1_14.txt, subtask1_15.txt, subtask1_16.txt, subtask1_17.txt, subtask1_18.txt, subtask1_19.txt, subtask2_01.txt, subtask2_02.txt, subtask2_03.txt, subtask2_04.txt, subtask2_05.txt, subtask2_06.txt, subtask2_07.txt, subtask2_08.txt, subtask2_09.txt, subtask2_10.txt, subtask2_11.txt
Case Name Status Exec Time Memory
subtask0_sample_01.txt AC 25 ms 924 KB
subtask0_sample_02.txt AC 23 ms 792 KB
subtask0_sample_03.txt WA 25 ms 920 KB
subtask1_01.txt AC 26 ms 796 KB
subtask1_02.txt AC 25 ms 920 KB
subtask1_03.txt AC 24 ms 928 KB
subtask1_04.txt AC 23 ms 928 KB
subtask1_05.txt AC 24 ms 924 KB
subtask1_06.txt AC 24 ms 924 KB
subtask1_07.txt AC 24 ms 924 KB
subtask1_08.txt AC 24 ms 804 KB
subtask1_09.txt AC 23 ms 796 KB
subtask1_10.txt AC 23 ms 928 KB
subtask1_11.txt AC 22 ms 748 KB
subtask1_12.txt AC 24 ms 796 KB
subtask1_13.txt AC 24 ms 928 KB
subtask1_14.txt AC 22 ms 932 KB
subtask1_15.txt AC 25 ms 800 KB
subtask1_16.txt AC 24 ms 800 KB
subtask1_17.txt AC 24 ms 920 KB
subtask1_18.txt AC 24 ms 928 KB
subtask1_19.txt AC 24 ms 928 KB
subtask2_01.txt WA 24 ms 672 KB
subtask2_02.txt WA 23 ms 672 KB
subtask2_03.txt WA 23 ms 796 KB
subtask2_04.txt WA 24 ms 700 KB
subtask2_05.txt WA 24 ms 924 KB
subtask2_06.txt WA 23 ms 800 KB
subtask2_07.txt WA 24 ms 928 KB
subtask2_08.txt WA 23 ms 924 KB
subtask2_09.txt WA 23 ms 928 KB
subtask2_10.txt WA 24 ms 924 KB
subtask2_11.txt WA 23 ms 672 KB