algorithm - Fast Integer Coordinates Inside/along a Circle Centered at origin with radius r -


consider have circle centered @ 0,0 radius r.

i'd integer points available inside circle. problem easy solve.

one may iterate on square x = -r +r , y =-r +r , see if x * x + y * y <= r * r, if so, add point result.

however, what's quickest way this? feel there should type of hack can take calculations (2r)^2 4/3 r^2

more particularly, have feeling can calculate length of inscribed square, add outer remaining components. i'm unsure of how though. math little dense. i'm refraining posting code because i'd general algorithm response, if 1 has preference, may state final benchmark used in should use jvm language.

any help?

note: similar gauss's circle problem, instead of counting number of points, want know points are.

you can values directly computing maximum y (the second coordinate of point on circle @ vertical of (x,0)) each value of x that:

for x in [-floor(r), floor(r)]     y_max = floor(sqrt(r^2 - x^2))    # pythagora's theorem     y in [-y_max, y_max]         # (x, y) !

i don't think can better (maybe can compute y_max faster won't big win) because anyway have these points in result.

edit:

this in pi*r^2 time, least can since it's number of points. can maybe save few computations doing quarter circle , getting other ones symmetry, i'm not sure it's faster, , it's longer write.


Comments

Post a Comment

Popular posts from this blog

commonjs - How to write a typescript definition file for a node module that exports a function? -

consider, toml node module can use: // toml.d.ts declare module toml { export function parse(value:string):any; } declare module "toml" { export = toml; } then: /// <reference path="../../../../../defs/toml/toml.d.ts"/> import toml = require('toml'); toml.parse(...); however, node module export single function, such 'glob' ( https://github.com/isaacs/node-glob ). the node usage of module be: var glob = require("glob") glob("*.js", {}, callback(err, files) { ... }); you'd naively expect this: // glob.d.ts declare function globs(paths:string, options:any, callback:{(err:any, files:string[]):void; ...but because typescripts 'import' semantics bit strange, seems can use 'import .. = require()' statement alias modules . trying call: /// <reference path="../../../../../defs/glob/glob.d.ts"/> import blog = require('glob'); results in: error ts2072: mod
Read more

openid - Okta: Failed to get authorization code through API call -

i'm integrating okta own idp server using okta's api. i'm implementing authorization code flow following steps below: in own server, use /api/v1/authn endpoint sessiontoken. use sessiontoken obtain authorization calling endpoint: /oauth2/v1/authorize?client_id=" + clientid + "&sessiontoken=" + sessiontoken + "&response_type=code&response_mode=query&scope=openid&redirect_uri=" + redirecturl + "&state=evanyang&nonce=" it's supposed return response status code 302 , location header containing redirect url code value. however, keep getting response status code 200 , without location header, html body saying "you using unsupported browser." , "javascript disabled on browser." according api documentation: http://developer.okta.com/docs/api/resources/oidc.html#authentication-request , sessiontoken parameter sufficient this: an okta one-time sessiontoken. allows api-based
Read more

thorough guide for profiling racket code -

did googling ("racket profiling", "racket measure performance"), didn't find , there no examples in docs. did google "profile" search on htdp - no luck. (profile (f ...)) output isn't obvious other small snippets. ideally want python's python -m cprofile usage examples. when search "python profiling", both duckduckgo , google give top result: 26.4. python profilers . (although scanned quickly, seems more of reference "thorough usage guide" examples. if had else in mind, perhaps please link that?) the equivalent racket documentation be: profile: statistical profiler . a usage example: #lang racket (require profile) (profile-thunk (thunk (function-to-profile arg0 arg1) )) here (thunk e) convenience (lambda () e) . bigger example: #lang racket (module mod racket (provide f) (define (f) (for/list ([i 10000]) i))) (require (prefix-in mod: 'mod)) (define (f) (for ([i 10000
Read more