#include "cielib.h" #include <cstdio> #include <iostream> #include <algorithm> #include <cstring> #include <vector> #include <set> #include <map> #include <cmath> #include <list> #include <ctime> #include <sstream> #include <queue> #include <stack> #include <bitset> #include <unordered_set> #include <unordered_map> using namespace std; typedef vector<int> vi; typedef pair<int,int> pii; typedef long long ll; typedef short int sint; #define FOR(x, b, e) for(int x=(b); x<=(e); ++x) #define FORD(x, b, e) for(int x=((int)(b))-1; x>=(e); --x) #define REP(x, n) for(int x=0; x<(n); ++x) #define ALL(c) c.begin(),c.end() #define SIZE(x) ((int)((x).size())) #define PB push_back #define ST first #define ND second #define mp(x,y) make_pair(x,y) #define DEBUG 1 #define debug(x) {if (DEBUG)cerr <<#x <<" = " <<x <<endl; } #define debugv(x) {if (DEBUG) {cerr <<#x <<" = "; FOREACH(it, (x)) cerr <<*it <<", "; cout <<endl; }} #define REMAX(a,b) (a)=max((a),(b)); #define REMIN(a,b) (a)=min((a),(b)); #define wez(n) int (n); scanf("%d",&(n)); #define wez2(n,m) int (n),(m); scanf("%d %d",&(n),&(m)); const int D = 501; int strz[D]; int tr[D]; int d, r, k; // czy zmiana z kt na delta jest ostro gorsza niz strz bool gorzej(int kt, int delta) { strz[kt] += delta; // if (strz[kt] < 0 || strz[kt] > r) printf("TROLOLOOL"); bool tres = czyCieplo(strz); strz[kt] -= delta; return czyCieplo(strz); } int main() { d = podajD(); r = podajR(); REP(i, d) { strz[i] = r/2; tr[i] = (r + 1) / 2; } int runda = 0; bool isDone = false; while (!isDone) { // printf("run/da: %d, tr: %d\n", runda++, tr[0]); isDone = true; REP(wymiar, d) { // printf("\twymiar: %d\n", wymiar); strz[wymiar] = min(strz[wymiar], r - tr[wymiar]); strz[wymiar] = max(strz[wymiar], tr[wymiar]); int prev = strz[wymiar]; strz[wymiar] -= tr[wymiar]; bool r1 = czyCieplo(strz); strz[wymiar] += 2 * tr[wymiar]; if (strz[wymiar] > r) { strz[wymiar] = r; } bool r2 = czyCieplo(strz); strz[wymiar] = prev; bool r3 = czyCieplo(strz); // printf("dol : %d, gora: %d\n", dol, gora); if (r2 && !r3) { strz[wymiar] = min(strz[wymiar] + tr[wymiar], r); } if (!r2 && r3) { strz[wymiar] = strz[wymiar] - tr[wymiar]; } int &temp = tr[wymiar]; if (temp > 1) temp = (temp + 1) / 2; else temp >>= 1; if (temp > 0) isDone = false; } } REP(i, d) tr[i] = 1; REP(wymiar, d) { strz[wymiar] = max(strz[wymiar], tr[wymiar]); strz[wymiar] = min(strz[wymiar], r - tr[wymiar]); bool dol = gorzej(wymiar, -tr[wymiar]); bool gora = gorzej(wymiar, tr[wymiar]); // printf("dol : %d, gora: %d\n", dol, gora); if (dol && gora) { // printf("oba gorsze\n"); // printf("wymiar: %d, strrzl: %d, tr: %d\n", wymiar, strz[wymiar], tr); // printf("\t--- binsearch size: %d\n", 2*tr); } else if (dol) { strz[wymiar] = strz[wymiar] + tr[wymiar]; } else if (gora) { strz[wymiar] = strz[wymiar] - tr[wymiar]; } } znalazlem(strz); }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 | #include "cielib.h" #include <cstdio> #include <iostream> #include <algorithm> #include <cstring> #include <vector> #include <set> #include <map> #include <cmath> #include <list> #include <ctime> #include <sstream> #include <queue> #include <stack> #include <bitset> #include <unordered_set> #include <unordered_map> using namespace std; typedef vector<int> vi; typedef pair<int,int> pii; typedef long long ll; typedef short int sint; #define FOR(x, b, e) for(int x=(b); x<=(e); ++x) #define FORD(x, b, e) for(int x=((int)(b))-1; x>=(e); --x) #define REP(x, n) for(int x=0; x<(n); ++x) #define ALL(c) c.begin(),c.end() #define SIZE(x) ((int)((x).size())) #define PB push_back #define ST first #define ND second #define mp(x,y) make_pair(x,y) #define DEBUG 1 #define debug(x) {if (DEBUG)cerr <<#x <<" = " <<x <<endl; } #define debugv(x) {if (DEBUG) {cerr <<#x <<" = "; FOREACH(it, (x)) cerr <<*it <<", "; cout <<endl; }} #define REMAX(a,b) (a)=max((a),(b)); #define REMIN(a,b) (a)=min((a),(b)); #define wez(n) int (n); scanf("%d",&(n)); #define wez2(n,m) int (n),(m); scanf("%d %d",&(n),&(m)); const int D = 501; int strz[D]; int tr[D]; int d, r, k; // czy zmiana z kt na delta jest ostro gorsza niz strz bool gorzej(int kt, int delta) { strz[kt] += delta; // if (strz[kt] < 0 || strz[kt] > r) printf("TROLOLOOL"); bool tres = czyCieplo(strz); strz[kt] -= delta; return czyCieplo(strz); } int main() { d = podajD(); r = podajR(); REP(i, d) { strz[i] = r/2; tr[i] = (r + 1) / 2; } int runda = 0; bool isDone = false; while (!isDone) { // printf("run/da: %d, tr: %d\n", runda++, tr[0]); isDone = true; REP(wymiar, d) { // printf("\twymiar: %d\n", wymiar); strz[wymiar] = min(strz[wymiar], r - tr[wymiar]); strz[wymiar] = max(strz[wymiar], tr[wymiar]); int prev = strz[wymiar]; strz[wymiar] -= tr[wymiar]; bool r1 = czyCieplo(strz); strz[wymiar] += 2 * tr[wymiar]; if (strz[wymiar] > r) { strz[wymiar] = r; } bool r2 = czyCieplo(strz); strz[wymiar] = prev; bool r3 = czyCieplo(strz); // printf("dol : %d, gora: %d\n", dol, gora); if (r2 && !r3) { strz[wymiar] = min(strz[wymiar] + tr[wymiar], r); } if (!r2 && r3) { strz[wymiar] = strz[wymiar] - tr[wymiar]; } int &temp = tr[wymiar]; if (temp > 1) temp = (temp + 1) / 2; else temp >>= 1; if (temp > 0) isDone = false; } } REP(i, d) tr[i] = 1; REP(wymiar, d) { strz[wymiar] = max(strz[wymiar], tr[wymiar]); strz[wymiar] = min(strz[wymiar], r - tr[wymiar]); bool dol = gorzej(wymiar, -tr[wymiar]); bool gora = gorzej(wymiar, tr[wymiar]); // printf("dol : %d, gora: %d\n", dol, gora); if (dol && gora) { // printf("oba gorsze\n"); // printf("wymiar: %d, strrzl: %d, tr: %d\n", wymiar, strz[wymiar], tr); // printf("\t--- binsearch size: %d\n", 2*tr); } else if (dol) { strz[wymiar] = strz[wymiar] + tr[wymiar]; } else if (gora) { strz[wymiar] = strz[wymiar] - tr[wymiar]; } } znalazlem(strz); } |