# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
881976 | dosts | 힘 센 거북 (IZhO11_turtle) | C++14 | 2083 ms | 10468 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
#define sp << " " <<
#define int long long
#define vi vector<int>
#define F(xxx,yyy) for (int xxx=1;xxx<=yyy;xxx++)
#define pii pair<int,int>
const int N = 1e5+1,inf = 1e18,P=2e5+1;
int MOD;
vi pri;
int sieve[P];
int add(int x,int y) {return (x+y >= MOD)?x+y-MOD:x+y;}
int mult(int x,int y) {return ((x%MOD)*(y%MOD))%MOD;}
int expo(int x,int y) {
if (!y) return 1;
int e = expo(x,y/2);
e = mult(e,e);
if (y&1) e = mult(e,x);
return e;
}
void prec() {
memset(sieve,0,sizeof sieve);
for (int i=2;i<P;i++) {
if (!sieve[i]){
sieve[i] = i;
pri.push_back(i);
for (int j=i*i;j<P;j+=i) sieve[j] = i;
}
}
}
int nck(int n,int k) {
if (n<k) return 0;
if (n<0 || k < 0) return 0;
int ans = 1;
for (auto it : pri) {
if (it > n) break;
int sm = 0;
for (int cur = it;cur<=n;cur*=it) sm+=n/cur;
for (int cur = it;cur<=k;cur*=it) sm-=k/cur;
for (int cur = it;cur<=n-k;cur*=it) sm-=(n-k)/cur;
ans = mult(ans,expo(it,sm));
}
return ans;
}
void solve() {
int n,m,k,t;
cin >> n >> m >> k >> t >> MOD;
prec();
n++;
m++;
vector<pii> pts(k+1);
F(i,k) cin >> pts[i].first >> pts[i].second;
F(i,k) pts[i].first++;
F(i,k) pts[i].second++;
sort(pts.begin()+1,pts.end());
int lim = (1<<k);
vi dp(lim,0);
for (int i=0;i<lim;i++) {
vi ys;
vector<pii> pt = {{1,1}};
for (int j=0;j<k;j++) if (i&(1<<j)) pt.push_back(pts[j+1]);
pt.push_back({n,m});
int sz = pt.size();
int p = 1;
for (int i=0;i<sz-1;i++) {
int dx = pt[i+1].first-pt[i].first;
int dy = pt[i+1].second-pt[i].second;
p = mult(p,nck(dx+dy,dx));
}
dp[i] = p;
}
int ans = 0;
for (int i=0;i<lim;i++) {
int pc = __builtin_popcountll(i);
if (pc > t && pc%2!=t%2) {
ans = add(ans,MOD-dp[i]);
}
else ans = add(ans,dp[i]);
}
cout << ans << endl;
}
signed main() {
ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#ifdef Local
freopen("in","r",stdin);
freopen("out","w",stdout);
#endif
int t = 1;
//cin >> t;
F(i,t) solve();
}
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |