Submission #3595348
Source Code Expand
#include <bits/stdc++.h>
using namespace std;
#define NDEBUG
#include <cassert>
typedef long long ll;
typedef long double Double;
typedef unsigned long long ull;
typedef pair<int, int> ii;
typedef pair<ll, ll> llll;
typedef pair<double, double> dd;
typedef vector<int> vi;
typedef vector<vector<int>> vvi;
typedef vector<ii> vii;
typedef vector<vector<ii>> vvii;
typedef vector<ll> vll;
typedef vector<string> vs;
typedef vector<double> vd;
typedef vector<long double> vD;
#define sz(a) int((a).size())
#define pb push_back
#define FOR(var, from, to) for (int var = (from); var <= (to); ++var)
#define rep(var, n) for (int var = 0; var < (n); ++var)
#define rep1(var, n) for (int var = 1; var <= (n); ++var)
#define repC2(vari, varj, n) \
for (int vari = 0; vari < (n)-1; ++vari) \
for (int varj = vari + 1; varj < (n); ++varj)
#define ALL(c) (c).begin(), (c).end()
#define RALL(c) (c).rbegin(), (c).rend()
#define tr(i, c) for (auto i = (c).begin(); i != (c).end(); ++i)
#define found(s, e) ((s).find(e) != (s).end())
#define mset(arr, val) memset(arr, val, sizeof(arr))
#define mid(x, y) ((x) + ((y) - (x)) / 2)
#define IN(x, a, b) ((a) <= (x) && (x) <= (b))
#define cons make_pair
ll gcd(ll a, ll b) {
while (a) swap(a, b %= a);
return b;
}
const ll MOD = 1000000007LL;
ll ADD(ll x, ll y) { return (x + y) % MOD; }
ll SUB(ll x, ll y) { return (x - y + MOD) % MOD; }
ll MUL(ll x, ll y) { return x * y % MOD; }
ll POW(ll x, ll e) {
ll v = 1;
for (; e; x = MUL(x, x), e >>= 1)
if (e & 1) v = MUL(v, x);
return v;
}
ll DIV(ll x, ll y) { /*assert(y%MOD!=0);*/
return MUL(x, POW(y, MOD - 2));
}
vi accum(vi& a) {
int L = a.size();
vi aA(L + 1, 0);
rep(i, L) aA[i + 1] = aA[i] + a[i];
return aA;
}
int accumed(vi& accum, int from, int to) {
int M = accum.size() - 1;
return accum[min(M, to)] - accum[max(0, from)];
}
void see(vvi& a, FILE* fp) {
int R = a.size(), C = a[0].size();
fprintf(fp, "# see(%d, %d)...\n", R, C);
assert(R <= 300 && C <= 300);
rep(r, R) {
rep(c, C) { fputc(a[r][c] ? '#' : '.', fp); }
fputc('\n', fp);
}
}
ll solve(int H, int W, vector<string>& s) {
int L = H + W - 1;
vvi up(L, vi(W)), down(L, vi(W));
rep(r, L) rep(c, W) {
up[r][c] = IN(r - c, 0, H - 1) ? (s[r - c][c] == '#') : 0;
down[r][W-1-c] = IN(r - c, 0, H - 1) ? (s[r-c][W-1-c] == '#') : 0;
}
vvi upA(L), downA(L);
rep(r, L) {
upA[r] = accum(up[r]);
downA[r] = accum(down[r]);
}
ll ans = 0, double_counted = 0;
rep(r, L) {
repC2(i, j, W) {
if (up[r][i] && up[r][j]) {
int d = j - i;
int above = r-2*d;
if (0 <= above) {
ans += accumed(upA[above], i-d, i+1);
if (IN(i-d, 0, W-1) && up[above][i-d]) ++double_counted;
if (IN(i, 0, W-1) && up[above][i]) ++double_counted;
}
int below = r+2*d;
if (below < L) {
ans += accumed(upA[below], j, j+d+1);
if (IN(j, 0, W-1) && up[below][j]) ++double_counted;
if (IN(j+d, 0, W-1) && up[below][j+d]) ++double_counted;
}
}
}
}
rep(r, L) {
repC2(i, j, W) {
if (down[r][i] && down[r][j]) {
int d = j - i;
int above = r - 2*d;
if (0 <= above) {
ans += accumed(downA[above], j, j+d+1);
if (IN(j, 0, W-1) && down[above][j]) ++double_counted;
if (IN(j+d, 0, W-1) && down[above][j+d]) ++double_counted;
}
int below = r + 2*d;
if (below < L) {
ans += accumed(downA[below], i-d, i+1);
if (IN(i-d, 0, W-1) && down[below][i-d]) ++double_counted;
if (IN(i, 0, W-1) && down[below][i]) ++double_counted;
}
}
}
}
return ans - double_counted/2;
}
vector<string> transpose(vector<string>& m) {
int R = m.size(), C = m[0].size();
vector<string> mT(C);
rep(i, C) {
vector<char> row(R);
rep(j, R) row[j] = m[j][i];
mT[i] = string(ALL(row));
}
return mT;
}
int main() {
int H, W;
cin >> H >> W;
vector<string> s(H);
rep(i, H) cin >> s[i];
if (H < W) {
s = transpose(s);
swap(H, W);
}
cout << solve(H, W, s) << endl;
return 0;
}
Submission Info
Submission Time |
|
Task |
E - Equilateral |
User |
naoya_t |
Language |
C++14 (GCC 5.4.1) |
Score |
700 |
Code Size |
4746 Byte |
Status |
AC |
Exec Time |
309 ms |
Memory |
3452 KB |
Judge Result
Set Name |
Sample |
All |
Score / Max Score |
0 / 0 |
700 / 700 |
Status |
|
|
Set Name |
Test Cases |
Sample |
s1.txt, s2.txt |
All |
01.txt, 02.txt, 03.txt, 04.txt, 05.txt, 06.txt, 07.txt, 08.txt, 09.txt, 10.txt, 11.txt, 12.txt, 13.txt, 14.txt, 15.txt, 16.txt, 17.txt, 18.txt, 19.txt, 20.txt, 21.txt, 22.txt, 23.txt, 24.txt, 25.txt, 26.txt, 27.txt, 28.txt, 29.txt, 30.txt, 31.txt, 32.txt, 33.txt, 34.txt, 35.txt, 36.txt, 37.txt, 38.txt, s1.txt, s2.txt |
Case Name |
Status |
Exec Time |
Memory |
01.txt |
AC |
135 ms |
3452 KB |
02.txt |
AC |
133 ms |
3328 KB |
03.txt |
AC |
132 ms |
3200 KB |
04.txt |
AC |
135 ms |
3200 KB |
05.txt |
AC |
133 ms |
3200 KB |
06.txt |
AC |
135 ms |
3328 KB |
07.txt |
AC |
132 ms |
3200 KB |
08.txt |
AC |
134 ms |
3200 KB |
09.txt |
AC |
64 ms |
3328 KB |
10.txt |
AC |
63 ms |
3328 KB |
11.txt |
AC |
62 ms |
3200 KB |
12.txt |
AC |
64 ms |
3328 KB |
13.txt |
AC |
64 ms |
3328 KB |
14.txt |
AC |
63 ms |
3328 KB |
15.txt |
AC |
259 ms |
3200 KB |
16.txt |
AC |
262 ms |
3200 KB |
17.txt |
AC |
261 ms |
3200 KB |
18.txt |
AC |
261 ms |
3328 KB |
19.txt |
AC |
257 ms |
3200 KB |
20.txt |
AC |
261 ms |
3200 KB |
21.txt |
AC |
309 ms |
3328 KB |
22.txt |
AC |
306 ms |
3328 KB |
23.txt |
AC |
57 ms |
3200 KB |
24.txt |
AC |
59 ms |
3328 KB |
25.txt |
AC |
182 ms |
3328 KB |
26.txt |
AC |
179 ms |
3328 KB |
27.txt |
AC |
1 ms |
384 KB |
28.txt |
AC |
1 ms |
384 KB |
29.txt |
AC |
1 ms |
384 KB |
30.txt |
AC |
1 ms |
384 KB |
31.txt |
AC |
1 ms |
384 KB |
32.txt |
AC |
1 ms |
384 KB |
33.txt |
AC |
1 ms |
384 KB |
34.txt |
AC |
1 ms |
384 KB |
35.txt |
AC |
1 ms |
256 KB |
36.txt |
AC |
1 ms |
256 KB |
37.txt |
AC |
1 ms |
256 KB |
38.txt |
AC |
1 ms |
256 KB |
s1.txt |
AC |
1 ms |
256 KB |
s2.txt |
AC |
1 ms |
256 KB |