#include #include #include #include #include using namespace std; typedef long long lint; typedef vectorvi; typedef pairpii; #define rep(i,n)for(int i=0;i<(int)(n);++i) const int N=4001; int x[N],y[N]; bool blk[N][N]; lint a[N][N]; lint b[N][N]; lint c[N][N]; int n,k; lint p; void disp(string s,lint a[N][N]){ cerr<0) tmp+=b[x+h][y-1]; } if(0){ cerr<<"f("<>n>>k>>p; rep(i,k){ cin>>x[i]>>y[i]; x[i]--,y[i]--; blk[x[i]][y[i]]=1; } memset(a,-1,sizeof(a)); queueque; rep(i,k){ que.push(pii(x[i],y[i])); a[x[i]][y[i]]=0; } int dx[6]={1,1,0,-1,-1,0}; int dy[6]={0,1,1,0,-1,-1}; while(!que.empty()){ pii t=que.front();que.pop(); int x=t.first,y=t.second; int d=a[x][y]; rep(dir,6){ int nx=x+dx[dir],ny=y+dy[dir]; if(nx<0||nx>=n||ny<0||ny>nx)continue; if(a[nx][ny]>=0)continue; a[nx][ny]=d+1; que.push(pii(nx,ny)); } } rep(j,n){ b[n-1][j]=a[n-1][j]; c[n-1][j]=a[n-1][j]; } for(int i=n-2;i>=0;--i){ rep(j,i+1){ b[i][j]=b[i+1][j]+a[i][j]; c[i][j]=c[i+1][j+1]+b[i][j]; } } if(0){ disp("a",a); disp("b",b); disp("c",c); } rep(i,n){ rep(j,i){ b[i][j+1]+=b[i][j]; } } lint ans=0; rep(i,n){ int lim=n-i; rep(j,i+1){ int hi=lim+1,lo=-1; while(hi-lo>1){ int mid=(hi+lo)/2; if(f(i,j,mid)>=p){ hi=mid; }else{ lo=mid; } } ans+=lim+1-hi; } } cout<