# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
13271 | dohyun0324 | 힘 센 거북 (IZhO11_turtle) | C++98 | 102 ms | 1876 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
int x,y,n,m,k,t,z,prime[30010],top,ch[300010];
ll d[21][21],d2[21][21],dis[21][21],arr[30010],dap;
struct data{
int x,y;
bool operator<(const data&r)const{
if(x==r.x) return y<r.y;
return x<r.x;
}
}a[300010];
ll cnt(int x,int p)
{
ll c=0;
while(x){x/=p; c+=x;}
return c;
}
ll comb(int x,int y)
{
int i,j;
ll ans=1;
for(i=1;i<=top;i++){
arr[i]=cnt(x,prime[i])-cnt(y,prime[i])-cnt(x-y,prime[i]);
for(j=1;j<=arr[i];j++){ans*=prime[i]; ans%=z;}
}
return ans;
}
int main()
{
int i,j,r;
scanf("%d %d %d %d %d",&n,&m,&k,&t,&z);
for(i=1;i<=k;i++) scanf("%d %d",&a[i].x,&a[i].y); k+=2;
a[k].x=n; a[k].y=m;
for(i=2;i<=300000;i++){
if(ch[i]==1) continue;
for(j=i;j<=300000;j+=i) ch[j]=1;
prime[++top]=i;
}
sort(a+1,a+k+1);
for(i=1;i<=k;i++){
for(j=i+1;j<=k;j++){
if(a[j].y>=a[i].y) dis[i][j]=comb(a[j].x-a[i].x+a[j].y-a[i].y,a[j].y-a[i].y);
}
}
for(i=1;i<=k;i++){
for(j=i+1;j<=k;j++){
d[i][j]=dis[i][j];
for(r=i+1;r<=j;r++){
d[i][j]-=d[i][r]*dis[r][j];
if(d[i][j]<0) d[i][j]+=z;
}
}
}
t++; d2[1][0]=1;
for(i=2;i<=k;i++){
for(j=1;j<=t;j++){
for(r=1;r<=i-1;r++){
d2[i][j]+=d2[r][j-1]*d[r][i];
d2[i][j]%=z;
}
}
}
for(i=1;i<=t;i++){dap+=d2[k][i]; dap%=z;}
printf("%lld",dap);
return 0;
}
컴파일 시 표준 에러 (stderr) 메시지
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |